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ABSTRACT 


The  Global  Positioning  System  (GPS)  is  revolutionizing  the  science  of  navigation.  Never  before  has 
there  been  a  system  that  could  provide  real  time,  world  wide,  continuous  coverage  with  the  such  precision. 
Yet,  the  accuracy  achievable  with  GPS  alone  is  not  sufficient  to  achieve  autonomous  flight  of  an  Unmanned 
Aerial  Vehicle.  However,  when  integrated  with  an  Inertial  Navigation  System  and  other  non-inertial  sensors 
using  a  Kalman  Filter,  GPS  supplies  the  critical  positioning  information  to  permit  such  on  achievement. 
This  thesis  presents  the  selection  and  integration  of  a  GPS  receiver  using  Differential  OPS  (DOPS) 


in  support  of  a  UAV  autonomous  flight  project.  Contemporary  electronic  navigation  systems  are  surveyed, 
OPS  operation  is  reviewed,  and  a  Motorola  PVT*6  GPS  receiver  selected.  Using  the  Motorola  Proprietary 
Binary  Format  protocol,  several  software  drivers  were  written  in  C  to  interface  the  information  to  an  Intel 
80486DX  CPU  using  the  RS-232  serial  communication  standard.  Finally,  an  examination  is  made  to 
determine  the  maximum  reacquisition  time,  the  DOPS  accuracies  achievable  and  the  effects  of  pseudorange 
correction  latency  on  DGPS  accuracy. 
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I .  INTRODUCTION 


The  last  of  a  24  satellite  constellation  was  launched  into 
orbit  in  June  1993,  and  with  it,  "the  science  of  navigation  is 
poised  to  take  a  great  leap  forward"  (Carey,  1993,  p.  Bl) .  The 
Global  Positioning  System  (GPS)  revolutionized  electronic 
navigation  by  providing  worldwide,  real-time,  continuous 
coverage  with  three  dimensional  accuracies  at  least  an  order 
of  magnitude  better  than  any  other  system  in  existence.  Even 
with  these  capabilities,  it  is  not  the  panacea  some  make  it 
out  to  be.  Absolute  accuracies  (i.e.,  aocuraoies  using  a 
single  receiver)  better  than  ten  meters  within  two  standard 
deviations  cannot  be  achieved. 

Differential  GPS  (DGPS)  shows  great  promise  by  further 
improving  accuracies  to  within  one  to  three  meters.  DGPS  uses 
two  GPS  receivers  where  one  is  placed  in  a  known  location  and 
the  other  is  allowed  to  move.  By  computing  the  difference 
between  the  stationary  receiver's  actual  and  GPS  determined 
locations,  and  then  transmitting  these  differences  to  the 
other  GPS  receiver,  the  other  GPS  receiver  can  correct  its 
output  to  obtain  the  improved  position  information.  This 
improvement  would  be  sufficient  for  an  aircraft  to  achieve  an 
autonomous  landing  with  GPS  alone  were  it  not  for  the 
unreliability  of  the  GPS  signals  in  high  speed,  critical 
maneuvering  evolutions.  It  is  not  uncommon  for  a  GPS  receiver 
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to  lose  track  of  one  or  more  of  its  satellites  for  a  several 
seconds,  and  thereby,  be  unable  to  provide  positioning 
information.  While  not  a  catastrophe  for  a  slow  moving  surface 
vehicle,  an  outage  of  even  a  few  seconds  could  prove  fatal  for 
an  aircraft. 

In  an  attempt  to  achieve  the  accuracy  and  reliability 
required  for  autonomous  flight  in  a  navigation  package  light 
enough  for  a  small  unmanned  aerial  vehicle  (UAV) ,  this 
research  effort  is  part  of  a  project  where  the  GPS  information 
is  combined  with  the  outputs  of  an  Inertial  Measurement  Unit 
(IMU)  and  with  data  from  other  non-inertial  navigation  sensors 
such  as  airspeed  and  altitude  indicators  in  a  Kalman  Filter. 
A  contribution  is  made  by  selecting  the  bast  GPS  receiver, 
defining  its  method  of  employment  and  writing  the  necessary 
software  drivers  to  smoothly  interface  them  with  a  Central 
Processing  Unit  (CPU) . 

In  the  following  chapter,  contemporary  electronic 
navigation  systems  are  surveyed,  discussing  their  advantages 
and  disadvantages.  Realizing  the  obvious  benefits  of  GPS, 
Chapter  III  examines  the  User  Segment  of  GPS  emphasizing  those 
details  necessary  to  intelligently  select  and  use  the  GPS 
receiver  in  a  manner  that  best  supports  the  various  project 
requirements.  Chapter  IV  brings  it  all  together  by  detailing 
the  receiver  selection  and  capabilities  as  well  as  the 
criteria  for  its  use.  It  also  describes  the  hardware  and 
software  necessary  to  interface  the  receivers  with  the  CPU. 
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Three  tests  were  conducted  to  confirm  expected  operational 
capabilities.  They  are  described  and  the  results  from  them  are 
presented  in  Chapter  V.  The  appendices  include  receiver 
specifications,  the  ephemeris  algorithm  and  the  software 
drivers  written  in  C  code  as  well  as  other  information  to  aid 
the  user  and  the  programmer. 
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II.  SURVEY  OF  ELECTRONIC  NAVIGATION  METKODS 


The  invention  of  wireless  communications  in  the  early 
twentieth  century  ushered  in  a  new  era  in  navigation  allowing 
seafarers  to  safely  pilot  their  ships  in  the  absence  of  sight 
and  sound  for  the  first  time  in  the  history  of  mankind.  The 
first  use  of  radio  waves  required  a  vessel  to  transmit  a 
continuous  wave  at  some  frequency.  On  shore,  receivers  with  a 
rudimentary  direction  finding  capability  triangulated  the 
ship's  position  and  transmitted  it  back  to  the  ship.  The 
position  was  not  very  accurate  by  today's  standards  but  served 
its  purpose  well  in  regions  of  the  world  where  the  sky  was 
commonly  obscured  by  clouds  or  fog  when  far  from  land. 

The  next  significant  advance  came  with  the  invention  of 
highly  accurate  chronometers.  First,  the  quartz  clock,  and 
then,  the  atomic  time  standards  with  mind  boggling  stabilities 
providing  accuracies  better  than  one  second  in  30,000  years  or 
one  part  in  a  trillion  (lo'*)  (Logsdon,  1992,  p.l56).  Such 
precision  is  required  to  allow  the  precise  measurement  of 
frequency.  In  addition.  Transit  (SatNav)  and  Global 
Positioning  System  (GPS)  partially  determine  position  by 
measuring  the  amount  of  time  fo'-  the  signal  to  travel  from  the 
satellite  to  the  receiver.  Radio  waves  travelling  at  the  speed 
of  light  will  advance  one  foot  in  one  billionth  of  a  second 


4 


(Kayton,  1990,  p.33e).  It  *s  not  difficult  to  see  the 

necessity  for  such  accuracy  in  time  keeping. 

in  order  to  appreciate  the  quantum  improvement  in  the 
science  of  navigation  that  GPS  has  achieved,  it  is  worthwhile 
to  review  the  major  electronic  navigation  methods  in  use 
today.  A  brief  review  of  some  principles  and  characteristics 
of  electromagnetic  waves  is  given,  followed  by  an  individual 
description  of  each  of  the  major  navigation  systems.  While  an 
effort  has  been  made  to  succinctly  describe  the  theory  upon 
which  the  system  operates,  the  important  point  is  to  notice 
the  range,  coverage  and  accuracy  for  each  method. 

A.  PRINCIPLES  AMD  CHARACTERISTICS  OF  ELECTROMAGNETIC  WAVES 
Radio  waves  in  a  vacuum  theoretically  travel  at  the  speed 
of  light.  Wavelength  is  related  to  frequency  by  Equation  (1) : 

X  =  2.997925  x  10°  m/ sec 

where  X  is  the  wavelength  in  meters  and  /  is  the  frequency  in 
Hertz.  Figure  l  demonstrates  this  relation  for  the 
electromagnetic  spectrum.  Information  may  be  carried  on  a 
radio  wave  through  one  or  more  modifications  to  its  shape. 
These  modifications  include  changing  the  amplitude,  frequency 
or  phase  of  the  radio  wave  and  are  appropriately  called 
amplitude,  frequency  and  phase  modulation,  respectively.  In 
addition,  a  wave  may  be  intermittently  transmitted  or  pulsed. 
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This  is  called  pulsed  modulation.  Figure  2  shows  three 
representative  examples  of  modulation  techniques. 


Figure  l.  Wavelength  Related  To  Frequency  In  The 
Electromagnetic  Spectrum.  (Hobbs,  1981,  p.  217) 


Radio  waves,  like  light,  experience  such  phenomena  as 
reflection,  absorption,  refraction,  diffraction  and 
interference.  With  respect  to  the  earth's  environment,  the  two 
most  influential  features  on  the  distance  and  path  radio  waves 
travel  are  the  earth's  surface  and  the  ionosphere.  The 


Flgura  2.  Modulation  Techniques  Commonly  Used. 
(Hobbs,  1981,  p.219) 


ionosphere  has  four  layers  and  exists  at  a  height  of 
approximately  60  to  400  Km  above  the  earth's  surface.  This 
highly  charged  and  fluctuating  segment  of  the  atmosphere  not 
only  absorbs  radio  wave  energy  but  also  reflects  it  (Beck, 
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1971,  p.95).  Waves  that  are  reflected  by  the  ionosphere  are 
called  "sky  waves".  The  ionosphere  is  affected  by  diurnal  and 
seasonal  fluctuations  and  by  the  level  of  sunspot  activity 
which  in  turn  affects  the  performance  characteristics  of  sky 
waves.  Sky  wave  phenomena  is  also  greatly  influenced  by  the 
frequency,  angle  of  incidence,  and  the  height  and  density  of 
the  various  layers  of  the  ionosphere.  The  sky  wave  frequencies 
most  commonly  observed  are  in  the  300  KHz  to  30  MHz  band 
(Medium  and  High  Frequency  bands) .  Because  of  the  higher 
resistance  of  the  earth's  crust  compared  to  the  atmosphere, 
radio  waves  transmitted  parallel  to  the  earth's  surface  are 
slowed  causing  the  wave  to  bend.  The  degree  of  bending  is 
inversely  proportional  to  the  frequency.  Thus  lower 
frequencies  will  travel  over  the  horizon  in  a  direct  path. 
These  radio  waves  are  called  ground  waves.  As  it  might  be 
expected,  a  receiver  might  travel  along  the  surface  of  the 
earth  from  the  transmitter  origin  and  receive  the  ground  wave, 
enter  into  an  area  shadowed  by  the  earth's  curvature  with  no 
signal,  and  then  pick  up  the  reflected  "one-hop"  sky  wave  even 
further  out.  Figure  3  illustrates  this.  Furthermore,  when 
conditions  are  right,  one-hop  sky  waves  may  experience 
reflection  by  the  earth's  surface  and  subsequent  ionospheric 
reflection  for  a  "two-hop"  sky  wave.  Navigation  systems  such 
as  Loran-C  and  Omega  use  reflected  sky  waves  (when  they  exist) 
to  extend  their  effective  range.  (Logsdon,  1992,  pp.  95-98) 
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7lgur«  3.  Sky  Wave  And  Ground  Wave  Patterns.  Ray  1  Is 
Refracted  Through  The  Ionosphere,  Ray  2  Is  A  One-Hop 
Sky wave,  Ray  3  Is  A  Two-Hop  Sky wave  And  Ray  4  Is  A  Ground 
Wave.  (Hobbs,  1981,  p.  225) 

Finally,  several  major  navigation  systems  exploit  radio 
wave  phase  measurement  In  order  to  determine  two  dimensional 
position.  The  principle  of  operation  uses  two  synchronized 
transmitters  located  a  fixed  distance  apart.  The  line  between 
thorn  Is  called  the  baseline.  As  a  result,  at  half-wavelength 
Intervals  the  phase  difference  between  the  two  signals  Is 
zero.  In  this  type  of  system,  the  half-wavelength  segments  are 
normally  called  lanes.  As  the  distance  from  the  baseline 
Increases  the  lines  of  equal  phase  or  "Isophase  lines"  radiate 
outward  from  the  foci  of  the  two  transmitters  in  hyperbolic 
fashion  as  shown  In  Figure  4.  Thus,  methods  employing  this 
technique  are  called  hyperbolic  radio  navigation  systems.  If 
the  two  transmitters  were  to  transmit  a  pulse  Instead  of  a 
continuous  wave,  then  by  measuring  the  time  of  arrival  of  each 
synchronized  pulse,  one's  position  might  also  be  determined. 
The  points  where  the  times  of  arrival  of  transmitted  pulses 
are  the  same  describe  a  hyperbolic  locus.  This  type  of  method 


is  also  considered  a  hyperbolic  radio  navigation  system.  With 
the  exception  of  radio  beacons,  all  contemporary  terrestrial 
radio  navigation  systems  use  some  sort  of  hyperbolic 
navigation  method.  (Bowditch,  1984,  pp.  1015-1016} 


Figure  4.  Construction  Of  A  Hyperbolic  Interference 
Navigation  Pattern.  (Hobbs,  1981,  p.  227) 

B.  RADIO  BEACONS 


Radio  beacons  are  continuous  wave  transmitters  modulated 
with  a  unique  Morse  code  identification.  They  are  usually  low 
power  transmissions  in  the  low  and  mid  frequency  bands  with  a 
range  of  less  than  200  miles.  Most  cannot  be  received  beyond 


20  miles  (Hobbs,  1981,  p.  231).  They  operate  using  shipboard 
radio  direction  finding  (RDF)  equipment  where  two  or  more 
beacons  are  used  to  obtain  a  fix.  The  accuracy  is  highly 
dependent  on  the  sensitivity  of  the  RDF  equipment  as  well  as 
the  range  from  the  beacons.  As  Figure  5  illustrates,  the 
geometry  of  the  radio  beacons  relative  to  the  ship  may  also 
play  a  significant  role  in  the  accuracy.  Finally,  the  precise 
positions  of  the  beacon  transmitters  must  be  known  or  the  line 
of  bearing  is  useless.  In  other  words,  the  position 
information  cannot  be  determined  in  terms  of  latitude  or 
longitude  until  after  it  is  plotted. 

C.  LORAM-C 

Loran-C  is  the  descendant  of  the  first  attempt  at  a 
hyperbolic  navigation  system  developed  during  World  War  II  for 
both  ships  and  aircraft  (Wilkes,  1987,  pp.  31-34).  From  Figure 
6,  it  oan  be  seen  that  it  covers  almost  the  entire  northern 
hemisphere  with  only  a  very  small  area  covered  below  the 
equator  and  then  only  by  sky  waves  when  they  exist. 

e 

Loran-C  operates  in  the  low  frequency  band  at  100  KHz  with 
a  10  KHz  bandwidth.  Within  geographic  areas,  there  is  one 
master  station  and  two  or  more  secondary  stations.  All 
stations  sequentially  transmit  a  pulse  modulated  continuous 
wave  in  groups  of  eight  pulses  separated  by  one  microsecond 
with  the  exception  of  the  master  station  which  transmits  nine 
pulses.  The  ninth  pulse  is  separated  from  the  first  eight  by 
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Vlgur*  S.  Intarsaction  Of  Radio  Beacon  Signals  Showing  Error 
Spread  As  Range  From  The  Source  Increases.  (Hobbs,  1981, 
p.235) 

two  microseconds  and  Is  used  primarily  for  visual 
Identification  when  viewing  the  pulses  on  an  oscilloscope. 
This  can  be  seen  in  Figure  7.  Each  group  of  master  and 
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Flgur*  6.  Ar«as  Of  Tha  World  Covered  By  Loran-C.  In 
Particular,  Note  That  The  Coverage  Is  Almost  Completely  In 
The  Northern  Hemisphere. (Hobbs,  1981,  p.  237) 


secondary  stations  is  called  a  chain.  Currently,  there  exist 
13  chains  throughout  the  world.  (USNO,  1972,  p.  7-5) (Hobbs, 
1981,  pp.  236-239) 

Each  master  and  secondary  station  transmits  according  to 
a  predetermined  fixed  order  synchronized  using  precise  atomic 
time  standards.  Within  each  receiver,  the  group  of  eight 
pulses  are  integrated  to  form  one  pulse  of  320  microseconds. 
The  time  difference  of  arrival  is  then  measured  between  the 
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QRI 


rigur*  7.  Loran-C  Pulat  sequance  For  A  Four  Station  Chain. 
(Hobba,  1981,  p.  239) 


master  and  selected  secondary  station  pulse  from  which  may  be 
plotted  a  line  of  position  on  Loran-C  charts.  As  shown  in 
Figure  8,  Loran»C  charts  have  the  number  coded  hyperbolic 
lines  overlaid  on  charts  and  to  enable  plotting.  Chain  and 
station  identification  are  achieved  by  measuring  the  pulse 
repetition  rate  (PRR)  and  the  pulse  repetition  interval  (PRI) . 
(USNO,  1972,  p.  7-5) (Hobbs,  1981,  pp.  236-239) 

Accuracy  of  the  system  varies  from  about  700  ft  near  the 
baseline  to  about  2000  ft  near  the  extreme  range  of  the 
system.  Due  to  the  low  frequency,  baseline  distances  range 
from  1000  to  1500  miles  so  that  coverage  within  a  chain 
includes  a  great  deal  of  area.  Similar  to  the  use  of  radio 


beacons,  Loran-C  charts  and  tables  are  required  for  the  signal 
information  to  be  meaningful. 
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Flgur*  t.  A  Portion  Of  A  Loran-C  Chart  Covorad  By  The 
Northeast  U.S.  Chain.  (GRI  9960}  (Hobbs,  1981,  242) 

D.  OBCGA  NAVIGATOR  SYSTEN 

Oacoa  was  conceived  by  an  American  and  developed  by  the 
British  Admiralty  Signals  Establishment  during  World  War  II 
for  use  in  mine  sweeping  and  the  Normandy  invasion  forces  in 
1944  (Beck,  1971,  p.  65).  Today,  it  is  operated  by  the  Racal- 
Decca  Navigator  Company  in  only  a  few  but  very  heavily 
traveled  areas  in  the  world  as  shown  in  Figure  9.  It  operates 
using  the  principles  of  hyperbolic  radio  navigation  but  on  a 
much  smaller  scale  than  Loran-C. 

Each  separate  area  of  coverage  has  a  master  station  and  up 
to  three  slave  stations  arranged  in  a  star  pattern  around  the 


Flgurt  9,  Araas  Of  Th«  World  Covared  By  Dacca.  (Beck,  1971, 
p.  64) 

mastar  atatlon  and  Idantlflad  by  a  color  designator  of  green, 
rad  or  purple.  Each  station  within  mastar/alave  chain 
transmits  an  integer  multiple  of  an  unmodulated,  continuous 
wave  base  frequency  (/)  in  the  range  of  14.00  KHz  to  14.3  3 
KHz.  The  master  station  transmits  at  6/,  the  purple  station  at 
5/,  the  red  station  at  8/  and  the  green  station  at  9/.  Each 
master/slave  pair  is  phase  locked  along  the  baseline.  A  Decca 
receiver  consists  of  four  receivers,  one  for  each  of  the 
frequencies  transmitted  by  the  master  and  slave  stations. 
(Beck,  1971,  pp.  65-67) 

Within  the  receiver,  each  master/color  frequency  is 
integer  multiplied  to  get  another  unique  multiple  of  the  base 
frequency  which  is  the  first  common  multiple  of  either  of  the 
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two.  For  instance,  the  master  station  frequency  of  6/  when 
multiplied  by  four  and  the  red  station  frequency  of  8/  when 
multiplied  by  three  produce  a  new  harmonic  of  24/  for  the 
master/red  station  pair.  Each  half  wave-length  of  this  new 
harmonic  is  considered  a  "Decca  lane",  and  by  comparing  phases 
of  the  two  signals,  one  can  determine  position  within  the  lane 
as  Figure  10  illustrates. 


Lane  ambiguities  are  resolved  once  each  twenty  seconds 
when  all  four  stations  transmit  the  fundamental  frequency  / 
whose  half-wavelength  precisely  encompasses  the  common 
multiples  (i8/  for  the  master/green  pair,  24/  for  the 
master/red  pair  and  30/  for  the  master/purple  pair) .  By  using 
phase  comparison  of  the  received  signal  in  each  of  these 
"Decca  zones"  as  shown  in  Figure  ll,  the  exact  lane  can  be 
determined.  In  a  like  manner,  zone  ambiguities  are  resolved  by 
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transmitting  on  a  frequency  of  8.2/  also  every  20  seconds. 
Called  an  orange  frequency,  it  forms  a  coarse  hyperbolic 
pattern  within  which  each  seo"  difference  cycle  encompasses 
five  zones.  (Hobbs,  1981,  pp.  252-253) 


Figure  ll.  Phase  Comparison  Of  Oecca  Master /Green-Slave  Lane 
Identification  Signals.  (Hobbs,  1981,  p.  252) 

The  effective  range  of  each  chain  is  approximately  240 
miles  from  the  master  station.  At  100  miles,  the  accuracy  is 
approximately  30  yds  by  day  and  lOO  yds  by  night.  Again,  this 
system  requires  the  Decca  charts  and  receiver  in  order  to 
provide  positioning  data.  (Bowditch,  1984,  pp. 1058-1059) 

E.  CONSOL 

Consol  was  developed  in  1945  from  a  German  radio  beacon 
system  called  Sonne.  It  is  a  hyperbolic  radio  navigation 
system  with  a  very  short  baseline  such  that  a  "degenerated" 


hyperbolic  pattern  is  formed.  The  curved  portions  of  the 


hyperbolic  pattern  are  not  used,  and  so  it  acts  as  a  long 
range  radio  beacon.  It  is  used  mostly  in  Europe  and  the 
Mediterranean.  Another  system  very  similar  to  Consol,  called 
Cosolan,  is  used  in  the  United  States.  Both  have  the  advantage 
that  no  special  equipment  is  required  other  than  a  low 
frequency  receiver.  The  operator  only  needs  to  listen  to  the 
number  of  dots  or  dashes  to  know  where  within  a  particular 
sector  the  receiver  is  located.  (Beck,  1971,  113-116) 

Both  systems  use  low  frequency  radio  waves  In  the  range 
from  190  KHz  to  350  KHz.  Each  station  consists  of  three  radio 
towers  spaced  three  wavelengths  apart.  One  of  the  towers 
transmli:s  a  continuous  wave  frequency  while  the  other  two 
transmit  waves  that  undergo  180*  phase  shifts  with  respect  to 
the  CW  wave  during  a  keying  cycle.  All  the  signals  are 
modulated  with  dots  and  dashes.  This  arrangement  produces  a 
series  of  sectors  10*  to  15*  wide  as  shown  in  Figure  12.  The 
phase  shift  of  the  combined  signal  results  in  a  variable 
number  of  dots  and  dashes  being  heard  depending  on  the  radial 
position  within  the  sector.  (Beck,  1971,  113-116) 

Consol  range  can  be  up  to  1500  miles  but  not  less  than  25 
miles  due  to  the  compressed  width  between  sectors  close  to  the 
transmitters.  Keeping  in  mind  that  only  radial  lines  of 
position  may  be  obtained  from  a  Consol  station,  accuracies  of 
less  than  a  degree  may  be  obtained  at  the  maximum  range 
improving  as  the  receiver  approaches  the  station.  Other  than 
knowledge  of  the  station  location  and  the  dot  and  dash 
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7-29) 

position  correlations,  little  is  required  to  use  Consol. 
(USNO,  1972,  p.  7-5) 

F.  VOR/DKE  AND  TACAN 

VOR/DME  stands  for  VHP  Oinni-directional  Ranging/  Direction 
Measuring  Equipment  and  TACAN  stands  for  Tactical  Air 
Navigation.  It  is,  in  effect,  a  radio  beacon  system  for 
aircraft  along  heavily  travelled  airline  routes  and  between 
military  airfields  and  ships.  The  principle  of  operation  of 
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TACAN  is  the  same  as  for  VOR/DME  except  that  TACAN  uses  higher 
frequencies  and  is  more  accurate  (Beck«  1971,  pp.  26*28).  Some 
VOR  stations  are  equipped  to  respond  to  TACAN  DME  frequencies 
and  are  called  VORTAC  stations.  VOR/DME  uses  VHF  frequencies 
between  108  MHz  and  118  MHz  and  TACAN  uses  UHF  frequencies 
between  960  MHz  and  1,215  MHz  (Logsdon,  1992,  p.  36). 
Consequently,  both  systems  are  line  of  sight  only. 

VOR  employs  two  different  signals  wor}cln9  together  In 
partnership:  a  narrow-beam  rotating  "lighthouse"  transmission 
coupled  with  an  "blln)clng"  omnl-dlrectlonal  pulse.  As  the 
lighthouse  transmission  sweeps  past  magnetic  north,  the 
blln)clng  pulse  Is  transmitted.  Receiver  equipment  marlcs  the 
time  of  receipt  for  the  first  pulse  and  then  the  time  of 
receipt  for  the  lighthouse  sweep,  since  the  lighthouse  sweep 
rotates  at  a  fixed  30  rev/see.  It  Is  a  simple  matter  to 
compute  the  bearing  from  magnetic  north.  VOR  provides  only 
bearing  information.  (Logsdon,  1992,  pp.  36-38) 

For  range  from  the  VOR  station,  the  aircraft  transmits  a 
series  of  Interrogation  pulses  to  the  DME  at  the  VOR  station 
which  then  transmits  a  reply  a  fixed  time  after  receipt. 
Equipment  on  the  aircraft  computes  the  slant  range  by  dividing 
the  time  from  Interrogation  to  the  time  of  receipt  of  tho  DME 
response  by  two  and  multiplying  by  the  speed  of  light.  Figure 
13  illustrates  this  operation.  (Logsdon,  1992,  pp.  36-38) 

At  the  extreme  range  of  receipt,  VOR  can  be  off  by  as  much 
as  three  degrees.  However,  as  the  range  to  the  transmitter 
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decreases,  the  error  is  reduced.  Overall,  VOR/DME  TACAN 
positioning  accuracies  range  from  approximately  200  to  600 
feet  (Logsdon,  1992,  p.  36).  The  obvious  advantage  over  a 
radio  beacon  system  is  ack  of  need  for  RDF  equipment  on  board 
the  aircraft. 

O.  OMSGA 

The  Omega  navigation  system  was  designed  to  provide 
worldwide  coverage  and  to  be  accurate.  It  consists  of  only 
eight  stations,  each  using  a  base  frequency  of  10.2  KHz,  and 
transmitting  at  10  Kw  power.  This  makes  it  possible  to  receive 
at  least  three  and  usually  four  stations  anywhere  on  earth 
with  an  accuracy  of  approx ImateJly  one  nautical  mile.  Omega's 
principle  of  operation  is  similar  to  Dacca  in  that  it  is  a 
hyperbolic  radio  navigation  system  that  uses  the  phase 
comparison  of  two  CW  transmissions  to  obtain  a  line  of 
position.  It  is  different  in  that  any  two  stations  can  be  used 
to  obtain  a  line  of  position  and  with  special  receiving 
equipment,  range  from  the  two  stations  for  a  fix.  (Bowditch, 
1984,  pp.  1016-1036) 

Two  10.2  KHz  continuous  waves  transmitted  exactly  in  phase 
but  travelling  in  opposite  directions  produce  a  series  of 
hyperbolic  lanes  called  Omega  lanes.  Each  lane  is  eight 
nautical  miles  wide  which  corresponds  to  the  half-wavelength 
of  the  frequency.  Phase  comparison  determines  the  receiver 
position  within  the  lane  which  is  expressed  in  terms  of 


22 


VOR;  OMNI  REFERENCE  PULSE  ♦ 
ROTATING  “LIGHTHOUSE" 
SIGNAL 

(FULL  SYSTEM  PROVIDES 
HEADING  AND  SLANT  RANGE) 


ROTATINO 


AMCAAfT 


DME:  2*WAY  ACTIVE 
SPHERICAL  RANGING 


SLANT  RANGE 


rlgur*  13.  VOR/DME  Principles  Of  Operation.  (Logsdon,  1992, 
p.37) 


centicycles  (cec)  or  centilanes  (cel) .  These  are  defined  to  be 
0.01  of  the  width  of  the  lane.  Figure  14  shows  two  Onega  lanes 
produced  by  the  10.2  KHz  frequency.  (Hobbs,  1981,  pp.  259-262) 
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rlgurs  14.  Two  Half-Wav«lttngth  Omaga  Lanas  Produoad  By  Tha 
Phase  Comparison  Of  Tha  10.2  KHz  Signals.  (Hobbs,  1981,  p. 
260) 


Lana  ambiguity  Is  resolved  by  multiplexing  tha  basic 
fraguanoy  with  three  others  at  11.05  KHz,  11.33  KHz  and  13.6 
KHz  over  a  tan  second  period  according  to  tha  commutation 
pattern  shown  In  Figure  15.  Wllthln  tha  receiver,  the 
difference  frequencies  of  0.283  KHz,  1.133  KHz  and  3.4  KHz  are 
extracted  which  are  1/36,  1/9  and  1/3,  the  frequency  of  the 
10.2  KHz  signal,  respectively.  The  3.4  KHz  signal  Is  used  to 
establish  a  broader  phase  lane  that  encompasses  three  Omega 
lanes.  Likewise,  the  1.133  KHz  signal  phase  lane  encompasses 
three  of  the  3.4  KHz  phase  lanes  and  so  on.  By  performing 
successive  phase  comparisons  as  demonstrated  In  Figure  16,  the 
operator  can  easily  determine  within  which  Omega  lane  his 
receiver  Is  located.  (Hobbs,  1981,  pp.  2S9--262) 


24 


Figure  IS.  Omega  Station  Commutation  Pattern. (Hobbs,  1981, 
p.  262) 


H.  SHIFTS  INERTIAL  NAVIOATION  SYSTEM  (SIMS) 

Although  it  requires  no  external  input,  SINS  is  still  an 

extremely  useful  navigation  system.  It  is  found  primarily  on 

submarines  and  aircraft  carriers.  Inertial  navigation  is... 

"...  the  process  of  directing  the  movements  of  a  vessel 
based  on  sensed  accelerations  in  known  spatial  directions 
by  means  of  instruments  that  mechanize  the  Newtonian  laws 
of  motion,  integrating  such  accelerations  to  determine 
velocity  and  position."  (Hobbs,  1981,  p.  281) 

In  essence,  an  inertial  navigation  system  such  as  SINS  is  a 

highly  accurate  dead  reckoning  computer. 
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Vigura  16.  Thr«a  Succassiva  Omega  Phase  Comparisons  For  Lane 
Resolution  on  Surface  Ships.  (Hobbs,  1981,  p.  261) 

SINS  mounts  two  accelerometers  on  a  platform  stabilized  by 
a  system  of  three  gyros  In  such  a  way  that  It  Is  constantly 
coincident  with  a  plane  tangential  to  the  earth's  surface.  The 
two  accalaromatars  are  continually  oriented  In  a  north-south 
and  east-west  direction;  hence,  they  are  sensitive  only  to 
horizontal  north-south  and  east-west  accelerations. 
Integrating  the  accelerations  with  respect  to  time  gives  the 
component  velocities  from  which  the  ship's  true  velocity  may 
be  computed.  (Hobbs,  1981,  p. 281-282) 

Because  of  the  cumulative  affect  of  a  number  of  possible 
sources  for  errors  such  as  friction  in  the  gyro  supports  and 
disturbances  caused  by  the  dally  rotation  of  the  earth,  SINS 
accuracy  degrades  over  time  and  must  be  updated  from  external 
sources,  one  improvement  to  SINS  has  been  the  Electro-Static 
Gyro  (ESG)  which  Is  used  to  reset  or  update  the  SINS.  Simply, 
the  ESG  is  a  one  centimeter  diameter  beryllium  sphere  spinning 


at  216,000  rpm  in  a  near  perfect  vacuum.  This  rotor  is 
supported  solely  by  an  electrostatic  field.  The  sphere  is  thus 
freed  from  the  classical  gyro  bearing  friction  as  well  as  many 
of  the  associated  random  torques  that  a  mechanical  support  can 
introduce.  The  ESG  extends  the  Interim  period  between  external 
updates  by  a  factor  of  six.  (Kayton,  1990,  pp.  194-206) 

Another  Improvement  in  sins  technology  is  the  Ring  Laser 
Gyro  (RLG) .  It  differs  from  the  traditional  concept  of  a  gyro 
In  that  it  contains  no  spinning  mass.  It  is  a  closed  geometric 
laser  path  (usually  triangular  such  as  the  one  in  Figure  17) 
centered  on  an  expected  spin  axis.  Identically  phased  laser 
beams  are  continuously  generated,  which  travel  in  opposite 
directions  around  the  closed  path.  Any  rotation  about  the  spin 
axis  causes  an  apparent  phase  difference  in  the  two  beams  at 
a  measurement  point  at  one  vertex  of  the  triangle, 
proportional  to  the  speed  of  rotation.  When  combined  with  the 
output  of  three  mutually  orthogonal  accelerometers,  position 
and  velocity  are  easily  computed.  RLG's  Increase  the  precision 
over  an  ESG  monitored  SINS  by  an  order  of  magnitude.  (Hobbs, 
1981,  p. 281-282) 

Finally,  in  an  effort  to  further  improve  accuracy  and 
reduce  size  and  weight,  the  latest  development  is  the  Fiber 
Optic  Gyro.  Its  principle  of  operation  the  same  as  RLG's 
except  the  laser  path  can  be  tightly  wound  in  a  spool  and 
mirrors  are  not  required.  The  technology  is  very  young  and 
promising  and  is  still  being  perfected. 
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Figure  17.  single  Axis  Ring*'Lasar  Gyro.  By  Monitoring  The 
Interference  Pattern  Caused  By  The  Two  Counter  Directional 
Isophase  Light  Beams ^  Extremely  Small  Rotations  Can  Be 
Measured.  (Logsdon,  1992,  p.  107) 

I.  8KIPBOARD  DOPPLBR  80MAR  SYSTEM 

Another  navigation  system  that  does  not  provide  absolute 
position  as  much  as  it  provides  accurate  inputs  to  a  dead 
reckoning  system  is  the  Shipboard  Doppler  Sonar  System. 
Because  of  its  extreme  accuracy  in  measuring  very  slow 
velocities,  it  is  found  on  very  large  ships  where  ship's  speed 
over  ground  is  important.  For  example,  when  mooring  to  a  pier, 
the  ship  could  crush  the  pier  if  the  contact  velocity  were 
greater  than  a  knot  (Hobbs,  1981,  p.  285) . 

It  operates  by  measuring  the  doppler  shift  of  two  or  more 
sonar  beams  in  order  to  compute  speed  over  ground.  One  such 


system  by  Raytheon,  called  Janus,  uses  two  fixed  hull  mounted 
transducers  that  transmit  pulsed  sonar  signals  in  two  beams, 
one  forward  and  the  other  aft,  at  about  a  30°  angle  from  the 
vertical.  Figure  18  illustrates  a  similar  configuration  that 
uses  four  transducers.  The  former  of  the  two  produces  fore  and 
aft  speed  accuracies  of  0.1  kt,  depth  accuracies  to  the 
nearest  foot,  meter  or  fathom  and  distance  accuracies  to  the 
nearest  hundredth  of  a  nautical  mile.  (Kayton,  1990,  pp.  188- 
193) 

J.  TRANSIT 

Transit  Navigation  Satellite  System,  sometimes  called  just 
Transit  or  SatNav,  can  trace  its  conceptual  origins  directly 
to  the  launch  of  Sputnik  I,  It  consists  of  five  satellites  in 
polar  orbits  at  an  altitude  of  580  nm  with  periods  of 
revolution  of  107  minutes  as  shown  in  Figure  19.  The  rest  of 
the  system  Includes  ground  tracking  stations,  a  comp*iting 
center,  an  injection  station  and  receivers.  Its  geometry  is 
configured  such  that  every  satellite  comes  within  range  of 
every  position  on  earth  at  least  twice  a  day  at  12  hour 
intervals.  The  average  time  between  fix  opportunities 
approximately  100  minutes  at  the  equator  and  a  half  an  hour 
near  the  poles.  The  difference  is  due  to  the  intrinsic  nature 
of  polar  orbits.  (Bowditch,  1984,  pp.  1066-1095) 

Transit  employs  two  frequencies  in  the  UHF  band  of  150  MHz 
and  400  MHz.  It  operates  by  using  an  integrated  doppler  method 
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Flgur*  18.  A  Dual-Axis  Four-Beam  Janus  Sonar  Array  Of  The 
Raytheon  Doppler  Sonar  System.  (Hobbs,  1981,  p.  2B5) 

where  the  slant  range  from  the  receiver  to  the  satellite  is 

computed  by  fitting  the  doppler  phase  shift  curve  of  the 

satellite  to  look  up  table.  Each  satellite  transmits  a  phase 

modulated  message  that  repeats  every  two  minutes  containing 

the  satellite  time,  ephemeris  and  other  information  from  which 


Flgur*  l!l.  The  Distinctive  ‘'Bird  Cage”  Polar  Orbits  Of 
Transit  Satellite  System.  (Hobbs,  1981,  p.275) 

the  receiver  calculates  the  receivers  position  in  space  and 
computes  the  longitude  and  latitude.  Within  in  each  pass  of 
the  satellite,  the  receiver  needs  to  collect  at  least  three 
successive  messages  for  a  fix.  Four  are  preferred  and  seven 
are  optimal.  The  number  receivable  is  dictated  by  the  height 
above  the  horizon  achieved  by  the  satellite.  (Bowditch,  1984, 
pp.  1066-1095) 


The  accuracy  of  Transit  is  affected  by  ships  dead 
reckoning  accuracy  and  speed  of  advance.  For  a  stationary 
ship,  it  is  common  to  achieve  accuracies  of  35  meters  made 
possible  by  the  shorter  wavelength  of  the  UHF  signals 
(Bowditch,  1984,  pp.  1067).  It  Is  significant  to  note  that 
while  this  system  is  world  wide  and  has  improved  positioning 
accuracies  of  better  than  an  order  of  magnitude  over  systems 
discussed  so  far,  it  is  comparable  to  Omega  in  its  coverage 
and  worse  in  its  availability  to  provide  positioning 
information. 

K.  GLOBAL  POfilTlONIKQ  SYSTEM  (OPS) 

With  the  June  1993  launch  of  the  twenty-fourth  satellite, 
NAVSTAR's  Global  Positioning  System  (GPS)  became  fully 
operational  in  late  1993.  GPS  provides  real  time  three 
dimensional  positioning  and  timing  information  anywhere  in  the 
world.  It  consists  of  24  uniformly  spaced  satellites  in  six 
orbital  planes  at  an  altitude  of  20,180  Km  as  depicted  in 
Figure  20.  Like  Transit,  it  uses  two  UHF  frequencies  at  1,575 
MHz  and  1227  MHz  simultaneously  upon  which  constellation 
ephemeris  data,  atmospheric  propagation  correction  data  and 
satellite  clock  error  are  phase  modulated.  (Wells,  1987,  pp. 
4.00-4.11) 

The  principle  of  operation  for  GPS  is  to  calculate  the 
range  from  four  different  satellites  using  the  broadcast 
message,  and  to  compute  the  longitude,  latitude  and  altitude 
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Figure  20.  GPS  Constellation  With  24  Satellites  Ensures  At 
Least  Twelve  In  View  At  All  Times  Anywhere  In  The  World. 
(Hobbs,  1981,  p.  288) 


of  the  simultaneous  solution  of  the  resulting  ranges.  Because 
the  system  is  operated  by  the  Department  of  Defense  for 
military  use,  the  lower  of  the  two  frequencies  is  encrypted 
and  altered  to  deny  any  potential  enemy  real  time  positioning 
accuracy  for  medium  and  long  range  targeting.  Still, 
accuracies  using  just  the  higher  frequency  are  approximately 
100  raeters  (Logsdon,  1992,  p.64).  Using  methods  to  be 
explained  in  the  next  chapter,  this  can  be  improved  by  more 
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than  an  order  of  magnitude  and  in  non-real  time  applications 
such  as  surveying,  mind  boggling  accuracies  of  a  millimeter 
spherical  error  probable  (SEP)  can  be  achieved  using  phase 
comparison  techniques  (Logsdon,  1992,  p.64) . 

L.  SUMMARY  OF  CONTEMPORARY  ELECTRONIC  NAVIOATIOM  SYSTEMS 

It  is  easy  to  see  why,  after  surveying  the  various 
electronic  navigation  systems  available,  GPS  invites  the  most 
enthusiasm  and  excitement  from  all  users.  It  is  a  quantum 
improvement  over  anything  else  available.  There  are  other 
navigation  systems  not  discussed  here  such  as  Microwave 
Landing  Systems  (MLS) ,  the  French  Argos  system  and  the  Russian 
GPS  imitation  GLONASS.  All  suffer  from  one  or  more  limitations 
in  range,  coverage,  availability  or  accuracy.  What  is 
important  in  this  review  is  an  appreciation  of  the  various 
contemporary  electronic  methods  and  their  method  of  operation 
including  their  particular  limitations.  In  the  next  chapter, 
the  GPS  principle  of  operation  is  more  closely  examined  as 
well  as  the  various  methods  for  its  use. 
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III.  GLOBAL  POSITIONING  SYSTEM 


GPS  operational  and  technical  characteristics  must  be 

understood  in  detail,  in  order  to  follow  the  reasoning  for  GPS 

receiver  selection  and  the  concept  of  employment  in 

( 

interfacing  them  with  the  CPt.\  In  the  following  sections,  a 
comprehensive  review  of  the  User  Segment  of  GPS  is  provided. 
It  covers  the  salient  points  of  signals,  coding,  message 
content,  message  structure,  calculations,  and  various  biases. 
The  chapter  ends  by  using  this  information  in  a  description 
and  discussion  of  Differential  GPS. 

GPS  is  functionally  divided  into  three  segments:  Control 
segment.  Space  segment  and  User  segment.  The  Control  segment 
exists  to  monitor,  update  and  operate  the  whole  system.  The 
Space  segment  consists  of  the  24  satellites.  These  two 
critically  important  segments,  however,  are  of  no  concern  to 
the  user  '^nd  little  is  required  to  be  known  of  them.  The  User 
segment  is  the  most  important,  fur  understanding  its  operation 
will  greatly  influence  the  user's  choice  of  receiver  and 
subsequent  method  of  employment. 

A.  USER  SEGMENT  RECEIVER  HARDWARE  DESCRIPTION 

Currently,  receivers  cost  as  little  as  a  few  hundred 
dollars  to  as  much  as  $50,000  or  more.  Size  and  weight  also 
vary  greatly  from  receivers  as  small  as  a  package  of 
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cigarettes  weighing  approximately  the  same  to  larger  than  a 
car  battery  and  weighing  only  slightly  less.  Shape  and  type  of 
antenna  can  also  impact  the  effectiveness  of  the  receiver. 
Figure  21  illustrates  the  most  common  designs  of  antennae 
found  on  the  market. 

One  type  of  antenna  called  a  microstrip  is  relatively  thin 
(approximately  1/10"  to  1/2"  thick)  and  can  be  designed  to 
receive  one  or  both  of  the  frequencies  that  GPS  satellites 
transmit.  This  design  is  ideal  for  use  in  this  project  where 
low  drag  is  desired  on  an  aircraft.  However,  signal  reception 
from  satellites  low  on  the  horizon  usually  experience  a  much 
lower  signal-to-noise  ratio  due  to  the  lower  gain  of  the 
antenna  at  those  low  angles  caused  by  its  thin  design. 

Other  factors  that  distinguish  receivers  from  one  another 
are  the  choice  of  position  solution  algorithm  and  interface 
communication  protocols.  All  of  the  above  characteristics  are 
determined  and  limited  by  the  state  of  the  art  of  the 
technology  and  by  the  purpose  for  which  it  was  designed. 

Receivers  may  be  loosely  grouped  into  four  broad 
categories  based  on  capability  (Logsdon,  1992,  p.  66): 

1.  Number  of  channels  and  sequencing  rate. 

2.  Access  to  selective  availability  signals. 

3.  Use  of  available  performance  enhancement  techniques. 

4.  Computer  processing  capabilities. 

The  majority  of  user  receiver  sets  sold  are  continuous 
tracking  or  parallel  receivers.  They  typically  have  from  four 
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Monopole  Quadrifilar  helix 


Spiral  helices 

Figure  21.  Four  Common  Types  Of  GPS  Antennae.  (Wells,  1987, 
p.  7.02) 


to  twelve  channels.  There  are  also  single  channel  receivers 
capable  of  providing  relatively  crude  positioning.  To  obtain 
data  from  multiple  satellites,  they  employ  one  of  three 
different  sequencing  speeds:  slow,  fast  or  multiplexed  (very 
fast)  as  shown  in  Figure  22.  (Logsdon,  1992,  p.  66-67) 
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MULTIIIPLEXINO  CIHANNBIL 

signals 

1 

0  fl  0  0  0  D  D  D  fl  D  D  D 

2 

0  D  D  D  D  D  D  D  0  0  fl  D 

a 

D  D  0  D  D  D  D  D  D  D  D  D 

n 

JL  D  0  0  D  D  D  D  D  D  D  0 

^ — 1  bit  period  (20  miiliseconds) 

FAST  SEQUBNCENQ  CMANHBL 

signals 

1 

1 — 1 

2 

CD  CD 

■ 

n 

L..  ^  l=l|  1=3 

seconds  (asynchronous  with  data  rate) 

SLOW  SBQUBNCONQ  CHANNEL 

signals 

1 

2 

e 

dp  dp 

dp^ 

n 

,  , 

t - minutes  to  hours 

(asynchronous  with  data  rate) 

Flgurs  22.  Multiplexing  And  Soquencing  Channels  For  GPS 
Receivers.  (Wells,  1987,  p.  7.13) 


Selective  Availability  (SA)  is  the  primary  means  by  which 
the  United  States  Department  of  Defense  denies  the  potential 
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terrorist  or  enemy  accurate  use  of  GPS  for  targeting 
navigation  for  missiles.  GPS  messages  are  encoded  using  two 
codes.  SA  alters  and  encrypts  the  message  content  of  the  more 
accurate  Precision  (P)  cods  carried  on  both  frequencies.  The 
less  accurate  Course  and  Acquisition  (C/A)  code  message  is 
called  the  Standard  Positioning  Service  (SPS)  and  is  available 
to  anyone  on  the  higher  frequency  only.  The  P-code  information 
is  called  the  l>reoi8e  Positioning  Service  (PPS)  .  only  those 
receivers  with  a  decoding  nip  (tightly  controlled  by  the 
DOD)  can  receive  PPS  information.  (Logsdon,  1992,  p.  67) 

SPS  receivers  can  still  achieve  excellent  accuracies  by 
using  performance  enhancement  techniques.  One  such  technique 
capable  of  centimeter  accuracy  uses  interferometry.  This 
process  needs  very  expensive  antennae  and  receivers,  and  the 
data  has  to  be  post-processed.  Surveyors  are  the  most  common 
users  with  this  accuracy  requirement  and  do  not  need  the  real 
time  positioning  Information  so  denial  of  PPS  is  not  a 
significant  handicap.  Another  capability  that  enhances  SPS 
accuracy  in  near  real-time  uses  a  process  called  Differential 
GPS  (DGPS)  .  Simply,  DGPS  employs  two  receivers:  one  stationary 
at  a  known  location  and  one  mobile.  The  stationary  receiver 
computes  the  difference  between  its  known  location  and  its  GPS 
derived  location.  It  then  transmits  those  differences  to  the 
mobile  receiver  which  uses  them  to  correct  its  own  GPS  derived 
position.  Since  the  Selective  Availability  errors  are  the  same 
over  baselines  as  large  as  600  nm,  accuracies  approaching  and 
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exceeding  those  of  PPS  receivers  can  be  achieved.  DGPS  is 
explained  more  completely  later.  (Logsdon,  1992,  p.  67) 

The  last  category  of  receivers  Is  distinguished  by 
computer  processing  capability.  Mostly  esoteric,  this  Is  where 
the  manufacturer  takes  the  GPS  message  and  calculates  the 
output  message  accounting  for  all  the  nuances  of  the  solution 
algorithm.  These  Include  such  Items  as  whether  or  not  to  use 
just  four  satellites  or  more  when  they  are  available,  whether 
or  not  to  account  for  Ionospheric  and  tropospheric  delays 
(sometimes  modeling  these  delays  can  Introduce  more  errors 
than  they  eliminate) ,  when  to  shift  to  other  satellites  that 
will  reduce  the  geometrical  error,  etc.  In  addition,  the 
receiver  can  be  enhanced  through  the  use  of  a  Kalman  Filter  to 
help  smooth  the  position  and  time  data  as  well  as  reduce 
receiver  processor  time.  (Logsdon,  1992,  p.  68) 

B.  QPB  8XQMAL8  AMD  CODXMO 

GPS  satellites  continuously  transmit  two  frequencies:  LI 
at  1575.42  Mhz  and  L2  at  1227.60  Mhz.  The  respective 
wavelengths  for  Ll  and  L2  are  approximately  20  cm  and  25  cm. 
Each  of  these  la  modulated  with  the  GPS  message  at  50  bps. 
Since  each  satellite  uses  the  same  frequencies  to  send  its  own 
unique  message,  the  receiver  must  be  able  to  distinguish  one 
signal  from  the  other.  Furthermore,  the  transmitting  power  of 
each  satellite  so  Is  low  that  by  the  time  the  signal  reaches 
the  receiver,  its  signal-to-noise  ratio  is  less  than  the 
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ambient  noise  level.  GPS  overcomes  these  obstacles  by  using 
spread  spectrum  and  code  division-multiple  access  techniques. 
As  mentioned  earlier,  the  messages  are  encoded  using  two  codes 
ca.lled  C/A  code  and  P-code.  The  chipping  rate  of  the  C/A  code 
and  P-code  is  1.023  Mbps  with  a  virtual  wavelength  of  290  m 
and  10.23  Mbps  with  a  virtual  wavelength  of  29  m, 
respectively.  Each  satellite  uses  its  own  unique  code  with 
which  the  receiver  may  auto-correlate  a  duplicate  and  thus 
phase  lock  on  to  it.  The  combination  of  the  transmission, 
coding  and  message  chipping  rates  provides  for  a  very  robust 
signal  with  which  it  is  very  difficult  to  Interfere.  (Wells, 
1987,  pp.  4.5  -  4.10) 

The  C/A  code  uses  a  1023  bit  Gold  code  with  a  period  of 
one  mllllseoond.  The  Gold  code  was  chosen  because  of  its 
length  and  because  of  the  property  that  Gold  codes  are  nearly 
orthogonal.  It  is  also  very  easy  to  replicate  and  with  its 
short  period,  easy  to  achieve  correlation  with  itself.  A  G  ' 
signal  can  only  be  received  if  the  code  assignment  is  known 
ahead  of  time.  The  Gold  code  assigned  to  each  satellite  is 
published  in  the  Interface  Control  Document  ICD-‘GPS-200 
published  by  Rockwell  Space  Systems  Division.  (Wells,  1987,  p. 
6.07) 

The  P-code  uses  a  much  longer  code  with 
235,469,592,765,000  bits  and  has  a  period  of  approximately  267 
days.  It  is  divided  into  32  seven  day  segments,  and  a  segment 
is  assigned  to  each  satellite.  The  DOD  further  encrypts  the  P- 
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code  to  produce  the  Y-code  to  deny  its  use  to  unauthorized  PPS 
users.  The  encryption  and  length  of  the  P-code  make  it 
extremely  difficult  to  correlate  to  without  keying  information 
and  decryption  hardware.  (Wells,  1987,  p.  6.07} 

The  GPS  message  is  phase  modulated  onto  its  respective 
codes  which  in  turn  are  phase  modulated  onto  their  respective 
carrier  frequencies.  As  mentioned  earlier,  the  Ll  frequency  is 
modulated  with  both  the  C/A  code  and  the  P-code.  This  is 
accomplished  by  phase  modulating  its  quadrature  component  with 
the  P-code.  This  feature  primarily  allows  the  military  to 
employ  PPS  redundancy  in  certain  applications.  The  table  in 
Figure  23  compares  the  important  points  of  the  code  and  data 
features.  (Wells,  1987,  pp.  €.3  -  6.8) 

0.  GPS  MBSSAGB  COMTBNT 

A  GPS  receiver  determines  its  position  by  computing  the 
position  of  four  or  more  satellites  within  view  and  then 
measuring  its  range  from  each.  By  solving  four  simultaneous 
equations,  the  receiver  calculates  its  coordinates  in  an  Earth 
Centered  Earth  Fixed  (ECEF)  cartesian  coordinate  system.  The 
system  chosen  is  defined  to  be  consi8T:ent  with  the  World 
Geodetic  Survey  of  1984  (WGS-34)  which  places  the  origin  at 
the  intersection  of  the  Z  axis  (through  the  earth's  center  of 
mass  and  parallel  to  the  north  pole  of  the  Conventional 
Terrestrial  Pole  (CTP))  and  the  J^-axis  (the  intersection  of 
the  WGS-84  Reference  Muridian  Plane  parallel  to  the  zero 
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Parameter 

C/A  Code 

PCode 

Data 

CHIPPINU  RATE 

SO  bits/suc 

SPATUL  LENGTH 
(PER  BIT) 

mmm 

6000km 
(3235  n.  m.) 

REPETITION 

INTERVAL 

7  days 

Not  applicable 

CODE  TYPE 

Gold  code 

257  day  pseudo 
random  code 

Nut  applicable 

TOTAL  NO.  OF 
CODES 

37  unique 
Gold  codes 

37  seven-day  sections 

Nut  applicable 

Easy  to 
acquire 

Slightly  more  accuiatc. 
Resistent  to  jamming 
and  spoofing. 

Rejection  of  multipath. 

Provides  handover  from 

C/A  to  P  cniio,  uphcmcris 
data  and  clock  correction. 

rigurt  23.  GPS  cod«  comparison.  (Applisd  Rsssaroh 
Laboratorlss,  1993,  p.  2.2) 


msridian  plans,  and  ths  plans  of  ths  CTP'a  squator) .  Ttis  V* 
axis  completes  the  ECEF  right  hand  orthogonal  coordinate 
system  measured  in  the  equatorial  plane  90”  to  the  east  as 
shown  in  Figure  24  (Applied  Research  Laboratories,  1993,  p. 
2.1).  Having  defined  the  reference  frame  for  GPS  position,  it 
is  necessary  to  understand  a  few  more  concepts  in  describing 
position  within  an  ECEF  system. 

The  GPS  receiver  translates  the  ECEF  coordinates  into 
geodetic  coordinates  for  the  user.  The  equ^^tions  to  translate 
coordinates  from  ECEF  to  geodetic  coordinates  and  vice-versa 
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rigur*  24.  GPS  Earth  Cantarsd  Earth  Flxad  Coordinata  System 
As  Defined  By  The  WGS  84.  (BIH  Stands  For  Bureau 
International  De  L'Hsure)  (Applied  Research  Laboratories, 
1993,  p.  2.1) 

are  provided  in  Appendix  A.  By  using  the  reference  ellipsoid 
model  of  the  world  defined  by  WGS  84,  GPS  is  capable  of 
relating  pin  point  positions  between  continents.  Before  the 
World  Geodetic  Surveys  of  1972  and  1984,  inter-continent 
surveying  was  virtually  impossible.  Cartographers  could  only 
establish  local  datums  to  which  all  surveyed  points  were 
reforenced.  These  datums  were  quite  accurate  for  their 


respective  geographical  regions  at  best.  One  of  the 
accomplishments  of  WGS  84  was  in  relating  all  the  different 
local  datums  to  the  ECEF  coordinate  system.  By  modelling  the 
datum  as  a  point  of  tangency  of  an  ellipsoid  as  shown  In 
Figure  25,  the  ellipsoid  origin  could  easily  be  referenced  to 
the  origin  of  the  ECEF  coordinate  system.  With  the  advent  of 
WGS  84  and  GPS,  the  Importance  of  the  local  datums  has  been 
diminishing.  Most  GPS  receivers  contain  a  library  of  the  major 
local  ellipsoid  references  to  the  WGS  ECEF  coordinate  system 
so  charts  and  maps  that  use  different  local  datums  can  be 
employed  with  GPS. 

Once  the  reference  frame  Is  established  and  understood, 
satellite  position  can  be  described  using  the  six  classical 
Kepler Ian  orbital  elements.  They  are  briefly  described  In 
Figure  26  but  they  are  not  sufficient  for  the  most  accurate 
description  of  the  satellite  orbit.  This  Is  because  the 
predicted  position  of  a  satellite  is  the  result  of  least** 
squares  curve  fitting  which  Is  not  valid  for  an  entire  orbit. 
Thus,  It  is  necessary  to  Include  several  other  parameters  to 
more  accurately  model  satellite  positions  in  space  and  time. 
All  the  parameters  are  listed  In  Appendix  B.  All  told,  there 
are  17  specific  parameters  Included  In  the  GPS  message. 

D.  GPS  CALCULATIONS 

The  equations  required  to  determine  the  ECEF  coordinates 
of  the  satellites  from  the  ephemeris  are  presented  in  Appendix 
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rlgur«  25.  Two  Dlmenaional  Illustration  Of  A  Theoretical 
Ellipsoid  Tangent  To  The  Geoid  At  An  Arbitrary  Datum  P. 
(Leick,  1990,  p.  189) 


B.  If  the  receivers  and  satellites  could  be  perfectly 
synchronized  with  GPS  time,  three  parameters  would  be  unknown: 
X,  Y,  Z.  It  would  then  be  possible  to  determine  position  using 
only  three  satellites.  However,  all  receivers  use  inexpensive, 
lightweight  and  comparatively  inaccurate  quax'tz  clocks  and 
each  has  a  bias  from  GPS  time  for  which  there  must  be  an 
accounting.  Therefore,  for  three  dimensional  position,  there 
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THE  SIX  KEPLERIAN  ORBITAL  ELEMENTS 


THE  SIX  KEPLERIAN  ORBITAL  ELEMENTS 

a.  ft.  l.aa>.T 


Figure  26.  Six  Classical  Keplarian  Orbital  Parameters. 
(Logsdon,  1992,  p.  143} 
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are  four  unknowns  and  four  simultaneous  equations  as  shown  In 
Equation  (2): 


Pk  = 

Pk  = 

pj  = 
= 


-  a;)*  +  -  i;)*  +  +  c  At, 

+  |jr!  _  +  (^2  .  2,)*  +  c  At, 

^{x^  -  x^)^  *  (y^  -  +  (z^  -  z,)^  +  c  A t, 

^{x*  -  +  (z^  -  z,)^  +  c  At, 


(2) 


where  P/  is  the  pseudorange  for  the  »cth  epoch  In  time,  #  is 
the  satellite  number,  c  is  the  speed  of  light  and  At, is  the 


receiver  clock  bias  from  GPS  time. 


Parenthetically,  it  should  be  noted  that  although  GPS  time 
is  synchronized  within  one  microsecond  of  UTC,-  it  does  not 
recognize  UTC  leap  seconds.  The  last  simultaneous  epoch  of  GPS 
and  UTC  time  was  In  January  1980.  Consequently,  GPS  time 
differs  from  UTC  by  the  integer  number  of  leap  seconds  since 
then. 

£.  OPS  MESSAGE  STRUCTURE 

The  entire  GPS  message  format  is  illustrated  in  Figures  27 
through  29.  It  is  defined  by  ICD-GPS~200  written  by  Rockwell 
International  who  was  the  prime  contractor  for  GPS.  The  basic 
message  unit  is  called  a  Frame.  It  is  subdivided  into  five 
subframes  which  are  in  turn  subdivided  into  ten  words  each.  A 
word  has  30  bits.  For  positioning,  a  receiver  requires  only 
subframes  one  through  three.  Those  contain  the  necessary 
ephemeris  and  clock  corrections.  However,  the  designers  saw  a 
need  to  provide  the  user  information  about  all  the  satellites 
in  the  entire  GPS  constellation.  They  used  subframes  four  and 
five  to  do  this.  Since  this  "Almanac"  of  extra  information 
exceeded  the  size  of  those  two  subframes  alone,  it  is 
transmitted  using  subframes  four  and  five  of  25  successive 
message  frames.  Each  pair  of  subframes  four  and  five  of  the  25 
successive  frames  is  referred  to  as  a  page.  The  entire  block 
of  25  pages  and  subframes  one  through  three  is  considered  a 
Master  Frame.  Figures  27  and  28  illustrate  this  best.  At  a 
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data  rate  of  50  bps,  a  frame  needs  3  0  seconds  and  a  master 
frame  requires  12.5  minutes  for  complete  transmission. 


BASIC  MESSAGE  UNIT  IS  ONE  FRAME  (ISOO  BITS  LONG) 


1  FRAME  -  5  SUBFRAMES 


1  SUBFRAME  -  10  WORDS 


1  WORD  -  30  BITS 


ONE  I  MASTER  FRAME  I  INCLUDES 
ALL  25  PAGES  OF  SUBFRAMES  A  Sc  5 
-  37,500  BITS  TAKING  12.5  MINUTES 


Subframes  4  and  5  hf 
have  25  PAGES  | 


Figure  27.  GPS  Message  Format  Breakdown.  (Rockwell,  1987,  p. 


7. 13) 


F.  GPS  BIAS 

There  are  numerous  factors  that  affect  GPS  performance 
accuracy.  Selective  Availability  is  primary  among  these.  While 
the  exact  nature  of  Selective  Availability  is  classified,  it 
is  sufficient  to  know  that  it  denies  single  receiver  position 
accuracy  of  better  than  100  meters  95%  of  the  time.  In 
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addition,  the  Ionosphere  and  Troposphere  affect  the  tine  of 
receipt  by  slowing  down  the  signal.  If  the  receiver  could 
receive  two  frequencies,  these  atmospheric  errors  could  be 
determined  and  eliminated.  For  the  receivers  that  are  capable 
of  only  single  frequency  reception  (most  receivers) , 
compensation  for  these  errors  Is  achieved  by  using 
mathematical  models  which  adequately  reduce  their  affects  for 
most  positioning  purposes. 


The  GPS  Data  Message  Content 


1  2  3 


Klagi  (L2  code  &  data;  week  #; 
satellite  accuracy  and  health) 

Age  of  data 

Satellite  clock  correction  coefficients 


Almanac  for  satellites  1>24  (pages  1*24) 
Health  of  satellites  1*24  (page  25) 


Almanac  for  satellites  25-32  (pages  2, 3. 4, 5, 7, 8,  9, 10) 
Ionospheric  model,  and  UTC  data  (page  18) 

Aotispoof  flag  -  32  satellites  (page  25) 

Satellite  configuration  •  32  satellites  (page  25) 

Health  of  satellites  25-32  (page  25) 


Figure  28.  GPS  Message  Content  Organization  Within  The 
Message  Format.  (Applied  Research  Laboratories,  1993,  p. 

2.2) 

Signal  reflections,  commonly  referred  to  as  maltipath,  are 


another  source  of  position  errors.  Multipath  interferes  with 
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the  fundamental  concept  of  measuring  the  time  from 
transmission  to  the  time  of  receipt  in  order  to  determine 
range  from  the  satellites.  Fortunately,  with  adequate 
filtering  and  by  paying  careful  attention  In  choosing  the 
antenna,  almost  all  the  effects  of  multipath  may  be 
eliminated. 


The  GPS  Data  Message 

50  Bits  Per  Second  Data  Rate 


*  Fufmit  of  FmiVM  Twcniy-Fiw  Uianin.  I 

**  123  MIouIm  Bifori  Thi  finiiie  Repeita.  j 

Figure  29  GPS  Message  Organization.  (Applied  Research 
Laboratories,  1993,  p.  2.2) 

Finally,  satellite  geometry  relative  to  the  receiver  plays 
a  role  second  only  to  Selective  Availability  In  determining 
the  precision  achieved.  The  measure  of  "goodness"  of  satellite 
geometry  Is  called  Geometric  Dilution  Of  Precision  or  GOOP. 
Optimal  GDOP  (the  lower  the  GDOP  number  the  better)  is  when 
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three  of  the  satellites  are  on  the  horizon  equally  spaced  120° 
apart  and  the  fourth  is  directly  over  head  as  shown  in  Figure 
30.  Since  the  constellation  is  designed  such  that  there  are 
usually  12  satellites  in  view  at  all  times,  a  low  GDOP  of  two 
or  three  can  regularly  be  achieved. 


/ 

Positioning 

accuracy 


=  DOP  •  ( 

Geometry 
(Dilution  of 
Precision) 


Measurement 

accuracy 


POOR  GDOP 
satellites  bunched 
together 


GOOD  GDOP 
(ideal  case) 

•  onr  satellite  overhead 

•  3  on  horizon, 

120**  apart  in  azimuth 


Figure  30.  Geometric  Relation  Between  GPS  Satellites  Affects 
The  Position  Accuracy.  (Wells,  1987,  p.  4.22} 


a.  DIFFERENTIAL  GPS 


SPS  does  not  guarantee  accuracies  less  than  300  m.  It  is 
possible  to  achieve  accuracies  two  orders  of  magnitude  better 


through  a  process  called  Differential  GPS  but  it  requires  at 
least  two  receivers!  The  concept  is  simple  and  generates  the 
greatest  amount  of  promise  in  commercial  applications  of  GPS. 
By  placing  one  of  the  two  receivers  in  a  pre-surveyed  location 
whera  its  geodetic  coordinates  are  known  as  precisely  as 
possible,  it  can  compare  its  true  position  to  the  position  GPS 
provides  and  calculate  the  differences  between  latitude, 
longitude  and  altitude.  If  the  resulting  error  signs  are 
changed,  the  errors  become  corrections  that  may  added  to  the 
second  receiver's  GPS  position  for  a  greatly  improved  position 
accuracy,  over  baselines  of  up  to  600  nm,  the  errors  are 
approximately  the  same. 

In  practice,  there  are  a  few  restrictions  in  employing 
DGPS.  The  most  important  is  that  the  position  results  from 
both  receivers  must  be  derived  using  the  same  satellites  or 
else  the  improvements  in  position  are  diminished  by  the 
introduction  of  new  errors.  The  other  restriction  is  that  the 
corrections  should  be  applied  at  the  earliest  point  in  the 
algorithm  used  to  determine  position  coordinates  in  order  to 
minimi^je  the  propagation  of  errors.  The  reasons  for  these 
restrictions  will  become  clear  after  a  brief  description  of 
the  algorithm  used  to  compute  DGPS  correotlons. 

The  stationary  GPS  receiver  is  receiver  is  required  to 
send  corrections  to  the  mobile  receiver.  It  would  be  simple  to 
subtract  the  true  longitude,  latitude  and  altitude  from  the 
geodetic  position  obtained  from  the  GPS  satellites.  However, 
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the  errors  introduced  from  SA,  satellite  geometry  and 
atmospheric  interference  affect  the  psuedoranges  from  each 
satellite  in  varying  degrees,  one  satellite  psuedorange  will 
always  have  a  greater  error  than  the  next.  Since  the  GPS 
position  is  calculated  using  a  least  squares  reduction  of  the 
four  simultaneous  equations  listed  in  equation  two,  ^1.1  the 
individual  psuedorange  errors  will  be  "blended"  thereby 
permitting  the  greatest  error  to  act  on  all  the  psuedoranges 
and  the  final  position  solution. 

In  order  to  avoid  this  problem,  the  stationary  receiver 
computes  the  expected  position  of  and  range  from  each 
satellite  (based  on  the  time  of  travel  of  the  individual 
signal)  and  compares  it  to  the  position  and  range  that  the  GPS 
messuge  asserts  in  its  ephemeris.  The  difference  between  the 
two  becomes  a  psuedorange  correction.  At  least  four  of  these 
corrections  are  transmitted  to  the  mobile  receiver  and  applied 
to  the  corresponding  psuedoranges  obtained  by  that  receiver. 
Then  the  now  psuedoranges  are  used  in  equation  two  to 
calculate  the  longitude,  latitude  and  altitude  of  the  mobile 
receiver.  Since  the  individual  corrections  are  applied  before 
equation  two,  the  final  position  is  more  accurate. 

In  practice,  GPS  makes  all  of  its  calculations  using  the 
Cartesian  Coordinate  ECEF  system  so  the  stationary  receiver 
must  translate  its  true  position  from  geodetic  coordinates  to 
ECEF  coordinates.  In  a  similar  manner,  the  mobile  receiver 
must  translate  its  corrected  position  from  ECEF  coordinates  to 
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geodetic  coordinates.  To  summarize,  the  following  steps  are 
involved  in  improving  the  GPS  position  using  DGPSj 

1.  Convert  the  true  geodetic  position  of  the  stationary  GPS 
receiver  to  ECEF  coordinates. 

2.  Compute  the  satellite  positions  from  the  ephemeris. 

3.  Compute  the  expected  satellite  positions  by  multiplying 
the  time  from  transmission  to  time  of  receipt  by  the  speed 
of  light. 

4.  Compare  the  received  pseudoranges  to  the  expected 
pseudoranges  and  calculate  the  pseudorange  corrections  for 
each  satellite. 

5.  Transmit  the  pseudorange  corrections  to  the  mobile  GPS 
receiver  and  add  to  the  received  pseudoranges. 

6.  Compute  the  mobile  receiver's  position  using  the 
corrected  pseudoranges. 

7.  Convert  the  mobile  receiver's  position  from  ECEF 
coordinates  to  geodetic  coordinates. 

H.  SUMMJUiy 

The  User  segment  of  GPS  draws  together  the  results  of 
several  decades  of  study  and  experimentation.  Hardware  design 
and  miniaturization  are  pushing  the  state  of  the  art.  Signal 
construction  and  frequency  selection  optimize  its  immunity  to 
noise  and  interference.  Orbital  kinematics  referenced  to  an 
ECEF  coordinate  system  simplify  and  standardize  the 
information  contained  in  the  GPS  message  structure.  Even  so, 
natural  and  man-made  biases  corrupt  the  integrity  and  accuracy 
of  GPS.  Through  DGPS  users  can  overcome  the  degradations 
caused  by  these  biases  and  achieve  real  time  accuracies 
unheard  of  heretofore.  Understanding  these  principles  upon 
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which  the  User  segment  is  founded  will  enable  the  user  to 
better  select  a  receiver  that  suits  his  purpose  and  integrate 
it  into  whatever  application  he  so  desires. 
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IV.  OPB  RECEIVER  SELECTION  AND  INTERFACE  WITH  CPU 

A.  MISSION  REQUIREMENTS  FOR  OPS 

Bafore  any  aarloua  naaroh  for  a  GPS  receiver  can  be  made, 
it  la  naoesaary  to  define  the  extent  to  which  it  will  be  used, 
Theae  criteria  come  from  the  mlaaion  atatemant.  This  project's 
goals  are  to  integrate  an  Inertial  Measurement  Unit  (IMU)  with 
GPS  and  other  non-lNS  sensors  through  a  smoothing  and  control 
filter  to  achieve  autonomous  flight  of  an  Unmanned  Aerial 
Vehicle  using  off  the  shelf  hardware  and  technology.  As  a 
result,  several  requirements  are  defined  for  the  GPS 
receivers. 

To  begin  with,  the  position  output  should  be  frequent 
enough  to  provide  regular  corrections  to  the  iMU.  Also,  the 
GPS  receiver  must  actually  update  its  position  solution  as 
often  as  possible  to  ensure  accuracy  in  the  information  of  its 
outputs.  The  receiver  that  is  to  be  in  the  aircraft  should  be 
small  and  lightweight.  It  should  be  easily  interfaced  to  the 
CPU,  and  it  should  provide  the  necessary  information  to 
conduct  phase  measurements  of  the  incoming  GPS  signal.  (This 
last  feature  should  provide  the  positioning  accuracy  to  auto¬ 
land  the  aircraft.)  Finally,  it  must  be  readily  available  on 
the  open  market.  Intrinsic  in  the  mission  statement  is  that  it 
must  be  affordable  given  the  project  budget.  With  theae 
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guidelines,  an  informed  selection  can  be  made  in  choosing  the 

best  GPS  receivers  for  the  job. 

B.  OPS  RECEIVER  SELECTION 

Having  reviewed  the  principles  of  GPS  operation  and 
variations  in  receiver  design  features,  an  extensive  survey  of 
216  different  GPS  receivers  from  54  companies  was  conducted. 
In  the  January  1993  issue  of  GPS  World  magazine,  many  of  the 
critical  facets  of  each  receiver  were  listed  for  comparison 
(Chan,  1993,  pp.  52-63).  In  order  to  stay  within  budget,  only 
those  reoeivera  less  than  $5,000.00  were  considered.  (The 
highest  priced  receiver  was  $55,000.00  with  many  between 
$10,000.00  to  $20,000.00.)  Among  the  remaining  units,  those 
that  were  capable  of  tracking  six  or  more  satellites 
simultaneously  were  selected  for  comparison.  This  left  the 
list  at  seventeen  receivers. 

Of  the  remaining  seventeen  receivers,  all  were  capable  of 
the  same  approximate  accuracy  and  information  output  rata  ac 
once  per  second.  The  one  notable  exception  was  Navstar's  XR5 
with  an  output  rate  at  1/10  second,  but  ultimately,  it  was  not 
chosen  because,  at  $1000.00,  it  cost  three  times  as  much  as 
the  Motorola  PVT-6.  Since  all  the  receiver  units  performed  at 
virtually  the  same  level,  the  deciding  factors  came  down  to 
cost,  weight,  size  and  power  consumption.  Table  I  shows  the 
top  eight  receivers  with  these  features  listed. 
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Tabl«  I.  COMPARISON  OF  TOP  EIGHT  CANDIDATES  FOR  THIS 
PROJECT'S  GPS  RECEIVERS.  (CHAN,  1993,  PP.  52-65) 


GPS 

Maker 

RCVR 

NAME 

Max» 

chan 

Size  (in^) 

Price 

Furuno 

OB-92/ 

ON-72 

a 

■ 

8.74 

2 

$750 

Magellan 

OPS 

Brain 

s 

B 

6.53 

B 

$445 

Magnavox 

Turbo  Sng 

6 

3.2 

11.44 

<1.5 

$3750 

Motorola 

PVT-6 

6 

B 

4.64 

1.3 

$369 

Navatar 

XRS 

Kernel 

12 

5 

8.28 

2.5 

$995 

Novatel 

Perform 

OEM 

10 

6.2 

15.44 

5 

$3195 

Sony 

m 

6 

8 

15.65 

2.61 

$1000 

Trimble 

SV6 

6 

2.9 

2.85 

1.B5 

$995 

The  Trimble  price  quote(i  includes  a  starter  evaluation 
kit.  The  Motorola  receiver,  when  all  options  were  installed, 
was  actually  more  expensive.  By  the  criteria  used,  Trimble  had 
the  better  unit  would  have  been  chosen  over  Motorola  but  they 
were  very  slow  in  answering  inquiries  for  characteristics  and 
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pricing  information.  When  the  decision  was  made,  they  were  not 
considered.  However,  two  factors  make  the  choice  of  Motorola's 
PVT*'€  the  best  selection.  One,  they  also  make  a  portable 
differential  GPS  base  station  called  a  LGT-1000.  It  is 
conceivable,  in  the  future,  this  project  will  choose  to 
expeziment  with  the  specialized  equipment.  And  two,  Trimble 
has  since  stated  they  do  not  intend  to  pursue  the  small, 
single  receiver  user  market.  Instead,  they  will  specialize  in 
industrial  scale  receivers  for  use  at  airports  and  in 
surveying. 

C.  OPS  RECBXVER  DESCRIPTION 
1 .  Hardware 

The  Motorola  PVT-6  is  an  original  equipment  for 
manufacturer  (OEM)  module.  Although  evaluation  software  is 
provided  in  a  starter  kit,  it  is  not  useful  when  the  CPU  is 
multi-tasking  with  several  other  processes.  PVT-6  stands  for 
Position,  Velocity  and  Time  with  six  channels.  Figure  31  shows 
the  dimensions  of  the  unit  without  the  antenna.  The 
specifications  are  listed  in  Appendix  C.  The  PVT-6  is  capable 
of  tracking  six  satellites  simultaneously.  It  uses  only  the 
C/A  coded  message  from  the  LI  carrier.  The  code  tracking  is 
carrier  aided.  It  can  be  powered  by  a  12  Vdc  unregulated  or  5 
Vdc  regulated  source.  (Motorola,  1992,  p.  3-2) 

The  PVT-6  is  supplied  with  a  low  profile  antenna  of  a 
microstrip  design  called  a  strip  patch  antenna.  The  antenna 
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Figure  31.  Physical  Dimensions  Of  The  PVT-6 . (Motorola,  1992, 
p.  3-3) 

gain  pattern  is  provided  in  Appendix  D.  The  Ll  frequency  is 
collected  by  the  antenna,  passed  through  a  narrow  band  filter, 
amplified  and  sent  to  the  receiver  unit  via  an  RG-58  coaxial 
cable.  It  provides  the  antenna  5  Vdc  @  25  mA  from  the  receiver 
and  returns  the  amplified  Ll  signal  from  the  antenna.  It 
should  be  noted  the  power  loss  along  the  cable  should  not 
exceed  6  dB  at  a  frequency  of  1575.42  MHz  which  equates  to  a 
maximum  .length  of  six  meters  for  RG-58  coaxial  cable.  RG-58  is 
the  recommended  cable  because  it  allows  a  sufficient  length  of 


cable  while  achieving  lighter  weight  due  to  the  small  gauge 
without  the  corresponding  losses  due  to  increased  resistance 
that  accompany  even  smaller  gauge  coaxial  cables  at  that 
frequency.  Antenna  placement  should  position  the  patch  plane 
to  be  as  level  as  possible  and  have  an  unobstructed  view  of 
the  horizon.  This  ensures  a  direct  line  of  sight  to  all  the 
visible  satellites.  Table  II  and  Figure  32  provide  the  power 
and  data  connector  pin  assignments.  (Motorola,  1992,  p.  D-1) 


Table  II  SERIAL  PORT  AND  POWER  PIN  ASSIGNMENTS.  (MOTOROLA, 
1992,  p.  3“2) _ 


Pin0 

Signal  Name 

Description 

1 

+12V/+5V  BATT 

(Optional)  -fS  Vdc  regulated  or  Vde 

unregulated  for  running  Real  Time  Clock  and 
retention  of  satellite  ephemeris  information 
stored  in  Static  RAM  memory. 

2 

+5V  !4AiN 

(Optional)  ‘*■5  Vdc  regulated  for  power 
requirements  of  the  entire  GPS  receiver (+12  Vdc 
signal  not  used). 

3 

ESBSBSEBEII 

4 

V  ' 

Flash  memory  (BPROM)  programming  voltage. 

B 

'fl2V  MAIN 

+12Vde  unregulated  for  power  requirements  of  the 
entire  ops  receiver. 

6 

1  PPS’ 

(Option  A)  1  pulse  per  second  output. 

7 

IPPS  RTN 

(Option  A)  1  pulse  per  second  return. 

a 

RS232  TXD 

Serial  RS232  data  output. 

9 

RS232  RXD 

Serial  KS232  data  input. 

10 

RS232  RTN 

Signal  return  for  RS232  signals. 

'Not  uied  (uaed  only  by  Motorola  for  updating  of  software). 

"1  Hx  pulse  train  with  rising  edge  coincident  with  UTC/GPS  1  second 
tick. 


As  Figure  33  shows,  the  LI  analog  signal  enters  the 
receiver  and  is  down  converted  to  an  intermediate  frequency 
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Figure  32.  Side  View  Of  The  PVT-S  Showing  The  Pin  Numbering 
Order.  (Motorola,  1992,  p.  3-2) 


(IF) .  The  IF  is  passed  to  the  6-channel  code  and  carrier 
correlator  section  where  it  is  digitized  and  split  for  input 
into  six  separate  channels  for  correlation,  filtering,  carrier 
tracking,  code  tracking  and  signal  detection.  The  output  is 
synchronously  routed  to  the  MPU  (PVT-6  CPU)  where  the 
information  is  processed  and  formatted.  Also  within  the  MPU  is 
the  RS-232  interface  for  full-duplex,  asynchronous  data 
communication  with  the  host  equipment.  (Motorola,  1992,  p.  2- 
2) 

The  PVT-6  is  configured  to  operate  at  9600  b" ud  with 
eight  bit  words,  one  start  and  stop  bit  and  no  parity.  Its 
communications  port  is  considered  Data  Communications 
Equipment  with  respect  to  the  RS-232  serial  communications 
standard.  However,  close  examination  of  the  available  pins 
for  serial  communications  reveals  the  presence  of  only  the  TX 
(pin  8)  ,  Rx  (pin  9)  and  Gnd  (pin  10)  pins.  This  dictates  there 
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rlgura  33.  Functional  Block  Diagram  Of  Tha  PVT-e.  (Motorola ; 
1992,  p.  2-1) 


will  be  no  modem  line  configuration  (DTK  and  DSR)  and  at  least 
no  hardware  flow  control  (CTS  and  RTS).  In  fact,  there  ie  no 
software  flow  control  (XON  and  XOFF)  either  implying  the  Data 
Terminal  Equipment  (DTE-  host  equipment)  must  monitor  the  port 
continuously  or  miss  valuable  Information.  Likewise,  the  PVT-6 
must  do  the  same.  In  practice,  it  receives  all  serial  input 
into  a  two  kilobyte  buffer  from  which  it  reads  750  bytes  per 
second.  The  buffer  is  large  enough  that  it  never  overflows. 
(Motorola,  1992,  pp.  3-2  -  4-3) 
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Finally,  the  PVT-6  is  provided  with  static  random 
access  memory  (SRAM)  for  retention  of  the  satellite  ephemerls 
data.  To  prevent  the  loss  of  this  information,  an  external 
power  source  such  as  a  battery  may  be  hooked  up  to  the 
+12V/+5V  BATT  connection  (pin  1)  when  the  GPS  unit  is  powered 
off.  In  addition,  nonvolatile  EEPROM  is  used  for  storage  of 
custom  operating  parameters,  almanac  information  and  other 
information  such  as  the  receiver  serial  number  and  option 
list.  (Motorola,  1992,  p.  2-2) 

2.  Motorola  Binary  Message  Format 

The  PVT-6  receiver  is  capable  of  providing  its 
information  in  any  one  of  three  formats:  Loran,  NMEA-0183  and 
Motorola  Proprietary  Binary  (Motorola,  1992,  pp.  4-2  -  4-6). 
Both  Loran  and  MMEA  (National  Marine  Electronics  Association) 
would  severely  restrict  this  project's  flexibility  and 
accuracy,  and  so  are  not  used.  Motorola  Binary  format  is 
logically  organized  and  provides  tremendous  flexibility  to  the 
user  to  apply  GPS  to  what  ever  purpose. 

Motorola's  Binary  format  has  47  different  messages  it 
sends  and  receives  ranging  in  length  from  7  to  294  bytes.  Each 
message  begins  with  (double  ASCII  0x40)  which  is  followed 
by  a  two  letter,  case  sensitive  message  ID.  Then,  0  to  287 
bytes  of  data  follow,  a  byte  size  check  sum  and  an  ASCII 
carriage  return-line  feed  (OxOD  OxOA) .  The  check  sum  is  the 
exc].usive-or  of  all  the  preceding  bytes  in  the  message.  The 
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data  that  follows  the  message  ID  can  be  thought  of  as  scaled 
floating  point  data  stored  in  integer  format.  This  avoids  the 
floating  point  storage  convention  as  well  as  saves  space  (only 
one  byte  is  needed  instead  of  four  or  more) .  The  user  must  be 
very  careful  to  observe  the  units  of  measurement  provided  in 
the  Motorola  GPS  Technical  Reference  Manual.  (Motorola,  1992, 
p.  4-2) 

For  every  command  to  the  GPS  receiver,  there  is  a 
response  usually  using  the  same  message  ID.  The  receiver  has 
no  way  to  determine  if  the  information  it  has  received  from 
the  CPU  is  correct  othar  than  by  verifying  the  several 
intrinsic  properties  of  the  message  itself: 

•  The  message  begins  with 

•  The  message  ends  with  CR-LF. 

•  The  length  of  the  message  matches  that  expected  for  the 
message  ID. 

•  The  checksum  is  valid. 

•  The  message  data  is  within  the  expected  maximum  and 
minimum  ranges. 

If  the  command  input  to  the  PVT-6  changing  a  parameter  is 
invalid  for  any  of  the  above  reasons,  it  responds  as  expected, 
but  the  data  values  will  not  reflect  the  ordered  change.  It  is 
prudent,  therefore,  to  compare  every  input  message  with  the 
subsequent  PVT-6  response.  (Motorola,  1992,  pp.  4-2  -  4-3) 

Input  commands  may  be  of  the  type  that  change  a 
particular  configuration  parameter  of  the  receiver.  Any 
command  with  a  message  ID  that  begins  with  a  upper  case  "A" 
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falls  into  this  category.  Input  commands  may  also  be  of  the 
type  that  enable  or  disable  the  output  of  data  or  status 
messages.  They  are  the  commands  that  the  CPU  will  use  for 
measuring  position,  velocity,  time,  pseudorange  and  satellite 
ephemeris  data.  These  output  status  messages  include  12  that 
are  the  work  horses  of  the  PVT-6.  Table  III  shows  only  nine 
because  three  of  the  commands  are  the  reciprocals  used  by  the 
CPU  to  talk  to  the  PVT-6.  The  nine  shown  are  listed  in  order 
of  precedence  for  transmission  from  the  GPS  receiver  to  the 
CPU.  (Motorola,  1992,  p.  4-3) 

Table  III  DATA  MESSAGE  OUTPUT  RATES  SHOWN  IN  ORDER  OF 
PRECEDENCE  OF  TRANSMISSION  ."ROM  THE  PVT-6  TO  THE  CPU. 
(MOTOROLA,  1992,  P.  4-4) 
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Whan  Raquaatad 

3 .  PVT-C  Operation 

The  PVT-6  has  some  limitations  that,  if  understood, 
may  be  overcome.  Data  is  output  once  per  second  from  the 
receiver.  This  information  is  calculated  from  the  ephemeris 


67 


and  time  measurements  of  the  previous  second  time  tic  called 
the  measurement  epoch.  In  order  to  provide  the  most  accurate 
position  given  the  latency  of  data,  the  position,  velocity  and 
time  represent  the  best  estimate  of  the  data  advanced  one 
second.  Figure  34  graphically  illustrates  this  process.  The 
first  byte  of  information  for  a  particular  measurement  epoch 
is  transmitted  anywhere  from  zero  to  50  ms  after  the 
subsequent  measurement  epoch.  Only  the  time  for  the  Satellite 
Range  Data  message  and  the  Pseudorange  Correction  message  is 
not  propagated  forward.  Careful  consideration  addressed  to 
this  latency  issue  in  the  design  of  any  Kalman  Filter  or 
smoothing  filter  that  integrates  GPS  position  with  the  IMU  and 
other  sensors  will  improve  performance.  (Motorola,  1992,  pp. 
4-13  -  4-17) 

To  assist  in  synchronisation  of  the  GPS  data  with 
other  sensor  data,  a  one  pulse  per  second  (IPPS)  signal  is 
output  using  pins  six  and  seven  of  the  serial  port.  The  pulse 
is  approximately  200  milliseconds  long  with  its  leading  edge 
corresponding  precisely  to  the  measurement  epoch.  The  PVT-6 
has  the  option  of  additionally  offsetting  the  IPPS  up  to 
0.999999999  seconds.  If  there  is  a  requirement  for  the  IPPS  to 
be  precisely  synchronized  /ith  the  UTC  or  GPS  measurement 
epoch,  the  user  may  also  offset  the  receiver  measurement  epoch 
as  well  as  account  for  antenna  cable  delay.  For  this  project, 
these  latter  conveniences  will  not  be  required.  This  is  shown 
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rlgurt  34.  Positlon/Channal/Status  Message  Data  Output 
Latency.  (Motorola,  1992,  p.4~is) 

in  the  timing  diagram  In  Figure  35.  (Motorola,  1992,  pp.  4''16 


-  4-17) 


Finally,  the  PVT-6  suffers  from  periodic  outages  like 


most  other  receivers  regardless  of  quality.  These  outages  may 
be  caused  by  temporary  loss  of  satellite  track  do  to  obscuring 
by  a  building,  tree  or  wing  tip  if  the  aircraft  were  in  a 
banking  turn.  Sudden  accelerations  may  also  cause  loss  of 
track.  In  this  event,  the  time  for  the  receiver  to  reacquire 
the  track  and  output  data  is  important.  As  shown  in  Table  IV, 
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I SECON 


output  Data.  (Motorola,  1992,  p.  4>17) 

this  tlma  nay  taka  as  long  as  IS  aaconds.  If  the  aircraft  were 
attempting  to  land  using  GPS  alone,  this  could  be  fatal. 
(Motorola,  1992,  p.  4-7) 

Table  IV  NOMINAL  REACQUISITION  TIMES  FOR  THE  PVT-6. 

(MOTOROLA,  1992,  P.  4-7) 
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D.  IHTERFACINO  THE  GPS  RECEIVERS 

The  two  GPS  receivers  to  be  used  for  DGPS  will  be  linked 
to  the  ground  CPU  and  the  UAV  CPU  through  RS-232  serial  ports. 
The  big  picture  of  the  project  showing  where  GPS  fits  in  is 
shown  in  Figure  36.  The  ground  control  station  will  be 
interfaced  through  a  33  MHz  Intel  804S6DX  personal  computer. 
It  will  be  the  primary  data  flow  conduit  for  the  uplink  and 
downlink  as  well  as  provide  formatting  and  mass  storage.  The 
additional  link  in  the  ground  station  shown  in  Figure  36  is 
for  redundancy  in  the  event  the  main  link  fails  and  emergency 


control  of  the  UAV  is  required. 


Figure  36  Archytas  Project  Control/Communications  Diagram 
Showing  Where  GPS  Fits  Into  The  Big  Picture. 
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The  UAV  CPU  uses  off  the  shelf  hariware  that  economizes  on 
size,  weight  and  price.  This  is  done  by  using  six  computer 
cards  and  a  passive  backplane.  The  CPU  is  a  3  3  MHz  Intel 
80486DX  with  256  Kb  double  cache  memory.  Non-volatile  memory 
for  the  operating  system  and  programs  is  provided  by  a  2.88  Mb 
RAM/ROM  Disk  card  instead  of  a  relatively  heavy  hard  drive.  In 
addition,  servo  control  and  sensing  are  managed  using  an 
Analog-to-Digital  (A/C) /Digital-to-Analog  (D/A)  card  and  a  ten 
channel  Timer/Counter  card.  Serial  I/O  is  handled  through  the 
two  standard  serial  communication  ports  on  the  CPU  card  (COMl 
and  COM2)  and  an  intelligent  8  port  RS-232  serial  I/O  card. 
The  GPS  receiver  communicates  to  the  CPU  through  port  1  on 
this  card  (not  to  be  confused  with  COMl  on  the  CPU  card) .  In 
order  save  space,  all  the  cards  are  connected  to  a  common  bus 
using  a  passive  backplane.  A  simplified  drawing  of  the  final 
arrangement  is  shown  in  Figure  37. 

1.  PCL-744  8-Fort  Intelligent  R8-237  interfaoe  Card 

The  UAV  CPU  needs  multiple  serial  ports  for  other  data 
inputs  from  the  link  modem,  the  IMU  and  the  non-INS  sensors. 
These  multiple  serial  inputs  are  provided  through  Advantech's 
PCL-744  8-Port  Intelligent  RS-232  Interface  Card. 

The  PCL-744  is  an  intelligent  card  in  that  it  has  its 
own  processor  and  64  Kb  RAM  buffer.  It  requires  only  one 
interrupt  number  fc  ;  eight  ports  although  each  may  configured 
independently.  This  configuration  includes  determining  the 
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Flgura  37  simplified  Illustration  Of  UAV  CPU  Hardware. 


size  of  the  transmit  and  receive  ring  buffers  and  what  type  of 
flow  control  is  used.  These  facts  enable  the  UAV  CPU  to 
perform  multi-tasking  by  allowing  it  to  access  the  ports  and 
process  data  without  losing  any  information  from  the  other 
ports  arriving  simultaneously.  The  PCL-744  takes  care  of  that 
by  multiplexing  between  its  eight  ports  and  placing  the 
incoming  data  in  buffers  that  can  be  accessed  by  the  UAV  CPU. 

The  PCL-744  comes  with  a  library  of  40  user  functions 
written  in  several  high  level  languages  including  Microsoft  c 
and  Turbo  C  by  Borland.  This  project  uses  the  Turbo  C  routines 
for  all  serial  communications.  These  functions  allow  the  user 
to  use  the  card  in  anyway  desired  for  serial  I/O.  It  is  also 
necessary  to  ensure  the  TSR  drivers  are  loaded  prior  to  PCL- 
744  operation.  They  are  called  STD-DRV.EXE  and  744-DRV.EXE  and 
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should  reside  in  the  AUTOEXEC.BAT  file  to  ensure  their 

loading. 

The  description  of  the  software  drivers  for  the  GPS 
receivers  that  follows  is  the  same  in  functional  design  for 
both  the  ground  and  UAV  units,  but  there  are  slight 
differences  in  coding.  This  is  due  to  the  fact  that  the  UAV 
receiver  must  communicate  through  the  PCL~744,  and  the  ground 
receiver  communicates  directly  with  the  CPU  through  a  standard 
DOS  COM  port.  The  port  interface  for  the  ground  CPU  uses  a 
standard  C  library  Interrupt  routine,  and  the  port  interface 
for  the  UAV  that  communicates  through  the  PCL-744  uses  the  PC- 
ComLIB  serial  communication  library  routines  provided  with  the 
card. 

2.  OPS  Software  Driver 

In  communicating  with  the  PVT-6,  all  transmissions  are 
without  any  flow  control.  It  is  necessary,  therefore,  to 
provide  a  storage  buffer  where  the  contents  of  the 
transmission  may  be  stored  for  holding  until  the  CPU  is  ready 
to  handle  the  data.  The  information  sent  to  the  PVT-6  is  put 
into  a  two  kilobyte  buffer  in  the  receiver  itself.  All  message 
handling  is  done  by  the  receiver  including  parsing, 
verification  and  interpretation,  and  is  transparent  to  the 
user.  The  user's  only  responsibility  is  to  ensure  that  the 
message  sent  to  the  receiver  is  properly  formatted.  However, 
when  receiving  messages  from  the  PVT-6,  parsing,  verification 


74 


and  interpretation  must  be  done  by  the  user.  This  allows 
flexibility  for  the  user  and  is  one  of  the  minor  reasons  for 
selecting  a  Motorola  GPS  receiver. 

This  section  describes  the  software  driver  required  to 
accept  information  from  the  GPS  receivers.  The  software  in 
Appendix  E  is  written  in  modular  format  to  allow  flexibility 
and  ease  of  handling.  It  is  structured  so  it  can  be  expanded 
and  built  upon.  The  following  paragraphs  describe  the  flow 
chart  shown  in  Figure  38  which  illustrate  the  concept  of  the 
GPS  software  driver  in  Appendix  E  and  the  structure  of  the 
software  driver. 

The  PCL-744  requires  the  eight  ports  to  be  initially 
configured  using  the  vendors  software.  This  configuration  is 
stored  in  the  PCL-744  driver  and  the  card  is  thus  reconfigured 
upon  startup  every  time  power  is  applied.  This  does  not  make 
the  card  ready  for  operation,  though.  Each  port  must  be 
individually  opened  and  set  to  the  communication  parameters 
required  by  its  user.  For  instance,  the  PVT-6  uses  9600  BAUD 
with  an  eight  bit  word,  no  parity  and  one  start/stop  bit.  In 
addition,  RS-232  modem  control  (i.e.,  DTR  and  RTS),  hardware 
control  (i.e.,  DTS  and  RTS)  and  software  flow  control  (i,e., 
XON/XOFF)  are  set  to  off.  This  is  all  accomplished  using  the 
function  "Open_ports'' .  At  the  conclusion  of  using  the  ports 
they  should  be  formally  closed  and  this  is  accomplished  by 
using  the  function  "Close  ports” . 
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Flgura  38  Functional  Diagram  Of  The  GPS  Software  Drivers  For 
Receiving  Information  From  The  PVT-6, 


The  main  functions  that  handle  incoming  information 
are  called  "Master^^gps”  for  the  ground  CPU  and  "Slave_gps''  for 
the  UAV  CPU.  The  ''Master__gps''  function  first  gets  the  buffer 
length  of  the  message  waiting  in  the  Rx  buffer  for  port  number 
one  on  the  PCL-744.  It  then  creates  a  buffer  of  exactly  that 
size  and  extracts  the  information  from  the  Rx  buffer  and  puts 
it  into  the  temporary  buffer  just  created.  Next,  it  creates 
the  input  buffer  which  will  contain  the  temporary  buffer 
appended  to  the  excess  buffer  contents.  The  excess  buffer  has 
the  partial  message  left  over  from  the  previous  call  to 
''Master_gps" .  With  the  input  buffer  created  and  the  contents 
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of  the  excess  and  temporary  buffers  copied  into  it, the  excess 
and  temporary  buffer  memory  is  freed. 

The  pointer  to  the  input  buffer  is  then  passed  to  the 
function  "Parse^inbuf f "  where  the  first  message  is  parsed  from 
the  input  buffer  and  returned.  If  a  complete  message  is  not 
returned,  the  function  places  the  partial  message  in  the 
static  excess  buffer,  frees  the  memory  for  the  input  buffer 
and  exits  the  function  "Master_gps'‘.  if  a  complete  message  is 
returned,  it  is  passed  to  the  function  ''M3g__verif ication” 
where  it  is  checked  for  the  correct  length,  check  sum,  ASCII 
prefix  and  hexadecimal  suffix  xOD  xOA  (ASCII  carriage 
return  and  line  feed) .  When  it  has  been  verified,  the  message 
is  placed  in  a  buffer  to  be  passed  out  to  the  calling  routine 
for  transmission  to  the  UAV  GPS  if  it  is  a  pseudorange 
correction.  Otherwise,  it  is  written  to  disk  for  post¬ 
processing. 

In  the  UAV  GPS  driver,  the  ”Slave_gps"  function 
operates  exactly  the  same  way  except  after  a  message  has  been 
verified,  it  is  put  into  a  structure  and  scaled  according  to 
the  format  in  the  Motorola  protocol  for  subsequent  use  by  the 
Kalman  filter.  One  of  the  more  complicated  structures  for  the 
ephemeris  message  has  been  written  as  an  example  by  Motorola 
and  incorporated  in  the  header  file  GPSCONV.H.  All  message 
output  from  the  UAV  GPS  is  passed  out  to  the  calling  routine 
to  be  down  linked  to  the  ground  in  binary  format  (to  reduce 
downlink  time)  for  post-processing. 
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Two  other  functions  are  mentioned  but  not  yet 
completely  written  which  play  an  important  role.  They  are  the 
''Initiali2e__master_gps''  and  •'lnitialize_slave_gps”  functions. 
They  will  pre-conf igure  the  receivers  for  such  things  as  true 
geodetic  coordinates  of  the  ground  antenna «  time,  date,  output 
rates,  desired  files  to  be  collected  etc.  They  should  be  made 
interactive  to  allow  the  user  to  change  and  set  any 
parameters.  The  code  framework  for  these  functions  is  in  the 
file  GPSPORTS.H  in  Appendix  E. 

Finally,  the  source  file  GPSAN.C  with  GPSCRNCH.H  was 
used  to  post-process  data  collected  for  this  thesis.  It  was 
written  to  take  the  Position/Channel/Status  message  from  an 
ASCII  file,  convert  it  to  binary  and  then  place  the  message  in 
a  structure  where  the  scaled  and  formatted  data  could  be 
analyzed.  Then  GPSCRNCH.H  "crunched"  the  numbers  and  computed 
the  maximum,  mean  and  standard  deviation  for  the  various  time 
delays. 
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V.  VERIFICATION  OF  OPERATIONAL  CAPABILITIES 


GPS  is  a  complicated  new  tool  for  navigators.  Each 
receiver  manufacturer  employs  unique  and  innovative  techniques 
to  improve  its  accuracy  and  flexibility.  Motorola  provides 
specifications  and  operational  capabilities  for  its  receiver 
and  users  have  published  various  accuracies  achieved  using 
DGPS.  It  is  prudent  to  verify  those  features  important  to  this 
project's  success.  Those  two  capabilities  are  examined  below. 
They  are: 

1.  Verification  of  the  manufacturer's  reacquisltlon  time 
specifications. 

2 .  Determination  of  the  effects  of  pseudorange  correction 
delay  in  getting  from  the  stationary  receiver  to  the  mobile 
receiver. 

In  the  following  sections,  these  concerns  are  investigated. 
For  each  case,  the  experimental  setup  is  described  and  the 
results  are  presented. 

A.  REACQUISITION  TIME  VERIFICATION 

This  project  Intends  to  use  DGPS  with  one  of  the  two 
antennae  on  a  UAV.  During  some  flight  maneuvers,  it  is 
expected  that  the  GPS  antenna  will  be  masked  from  one  or  more 
of  the  satellites  it  is  using  to  calculate  its  position.  It  is 
important  for  the  designers  of  the  Kalman  Filter  to  understand 
the  normal  reacquisition  times  given  such  a  circumstance.  In 
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Table  IV  in  Chapter  IV,  Motorola  presents  its  results  for  the 
reacquisition  time  of  the  PVT-6  receiver  given  the  antenna  was 
masked  for  15,  30,  45  and  60  seconds.  The  method  they  used  to 
collect  this  data  is  not  well  documented,  and  it  is  not 
indicated  whether  or  not  the  receiver  was  in  three  dimensional 
fix  mode  or  two  dimensional  fix  mode.  The  difference  could  be 
significant  because  the  former  requires  four  satellites  for  a 
fix  and  the  latter  only  three. 

The  reacquisition  time  was  investigated  using  this 
project's  PVT-6  receiver.  Table  V  shows  the  results  next  to 
the  Motorola  claims.  Data  was  collected  using  the  software 
supplied  by  Motorola  with  an  evaluation  kit.  The  antenna  was 
completely  masked  360*  around  above  and  below  using  four  to 
six  layers  of  tin  foil.  The  mask  was  left  in  place  for  the 
specified  amount  of  time  and  then  removed  for  at  least  30 
seconds.  Data  files  were  imported  into  a  spread  sheet  where 
the  time  marked  data  was  analyzed. 

The  results  do  not  exactly  agree  with  those  that  Motorola 
publishes,  but  they  are  close.  Significantly,  data  obtained 
during  daylight  hours  were  within  one  standard  deviation  of 
Motorola's  claimed  reacquisition  times  whereas  data  collected 
after  daylight  hours  were  not.  No  effort  was  made  to  confirm 
this  suggested  phenomenon  nor  was  any  effort  made  to  correlate 
it  to  ionospheric  thickness  or  tropospheric  conditions  such  as 
rain  or  cloudiness.  However,  until  the  IMU  is  integrated  with 
GPS  using  a  Kalman  Filter  and  the  effects  of  a  two  to  three 
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second  delay  in  reacquisition  time  can  be  determined, 

investigation  of  this  discovery  is  left  for  further  study. 

Table  V  COMPARISON  OF  MOTOROLA  AND  EMPIRICAL  REACQUISITION 
TIMES. 


Time  Obscured 
(sec) 

Motorola 
Reacquisition 
Time  (sec) 

Empirical 
Reacquisition 
Time  (sec) 

Standard 

Deviation 

(sec) 

IS 

<12 

13.75* 

2.2 

30 

<13 

17.08** 

2.0 

45 

<14 

16.67** 

3.2 

60 

<15 

13.67* 

2.1 

*  TJiea*  obB0rv*tionM  were  teAen  from  1230  to  1500, 
**  TAeee  oba^rvtionM  were  takan  from  2200  to  2400, 


B.  DOPS  ACCURACY  AND  PSUEDORANGE  CORRECTION  LATENCY 

The  accuracy  that  DGPS  can  achieve  is  well  documented.  The 
PVT-6  is  designed  and  built  to  permit  DGPS.  In  fact,  the 
Motorola  Binary  format  includes  master/slave  DGPS  commands  and 
automatically  incorporates  the  data  in  its  calculations.  It  is 
worthwhile  to  validate  the  expected  results  especially  in 
light  of  the  fact  the  communications  software  to  be  used  is 
not  written  by  the  manufacturers  and  may  adversely  affect  the 
accuracy. 

In  designing  the  test  for  DGPS,  the  length  of  che  baseline 
between  the  two  receivers  was  an  issue  that  was  researched. 
The  results  indicated  the  baseline  range  between  the  two 
receivers  did  not  significantly  affect  DGPS  accuracy  up  to 
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about  50  km.  Beyond  50  km,  out  to  about  1000  km,  the  effects 
slowly  increased  the  error  (Kremer,  1990,  pp.  307-312).  This 
project's  maximum  baseline  is  not  expected  to  exceed  15  km. 
For  all  intents  and  purposes,  a  zero  meter  base  line  is  as 
good  as  a  1000  meter  baseline  and  certainly  is  easier  to 
test.  Therefore,  this  test  simply  mounted  the  two  antennae 
side  by  side  at  a  pre-surveyed  location  (See  Appendix  F)  . 
Data  was  collected  by  sending  pseudorange  corrections  from  the 
master  antenna  to  the  slave  antenna  using  the  software  driver 
in  Appendix  E.  The  PVT-6  puts  out  a  position  message  at  a  rate 
of  once  per  second.  These  positions  were  differenced  using  the 
ECEF  equations  in  Appendix  A,  and  the  mean,  maximum  and 
standard  deviation  were  calculated. 

Table  VI  shows  the  DGPS  accuracy  with  various  delays  in 
transmitting  the  correction  to  the  slave  receiver  and  Figure 
39  graphically  illustrates  the  same  results.  The  delay  was 
artificially  introduced  through  the  function  call  to 
''Pr__delay''  in  the  header  file  GPSDELAY.H.  The  delays  were  in 
integer  amounts  of  seconds  as  shown  below. 

Interestingly,  Motorola  technicians  quote  adverse  effects 
of  pseudorange  corrections  received  after  90  seconds.  The 
results  above  suggest  an  immediate  degrading  effect  of  any 
amount  of  delay  and  that  the  error  grows  with  delay  time,  it 
is  hypothesized  that  the  decrease  in  error  at  the  15  second 
delay  point  is  a  spurious  result  caused  by  any  number  of 
factors.  It  is  left  for  further  study  to  confirm  this. 
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Table  VI  OBSERVATION  OF  PSEUDORANGE  CORRECTION  INPUT  DELAY 
ON  THE  EFFECTIVE  ACCURACY  OF  DIFFERENTIAL  GPS. 


Delay 

(sec) 


Average 

Maximum 

Standard 

Offset 

Offset 

Deviation 

(m) 

(m) 

(m) 

1.17 

3.76 

1.00 

2.58 

5.16 

0.87 

2.75 

3.84 

1.00 

2.72 

3.86 

0.96 

0.96 

3.28 

0.92 

2.41 

5.53 

0.97 

1.85 

3.85 

1.39 

2.45 

8.03 

2.02 

4.07 

8.83 

2.14 

4.51 

10.80 

2.19 

4.94 

13,85 

3.36 

1.96 

5.91 

1.48 

6.09 

21.21 

4.59 

21.40 

42.16 

10.90 

6 


Showing  The  Effects  Of  Time  On  The  Accuracy  Of  DGPS 
Corrections. 
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VI.  CONCLUSION 


This  thesis  reviews  the  major  contemporary  electronic 
navigation  systems  available  and  demonstrates  that  GPS  offers 
the  best  accuracies  and  world  wide  coverage  on  a  realtime 
basis.  The  user  segment  of  GPS  Is  explained,  and  Its 
capabilities  and  deficiencies  discussed.  Differential  GPS  Is 
presented  and  shown  to  be  able  to  achieve  enhanced  accuracies 
sufficient  to  significantly  contribute  to  this  project. 

An  extensive  review  of  available  GPS  receivers  was 
conducted  and  Motorola's  PVT-6  was  found  to  combine  all  the 
benefits  of  GPS  in  a  small,  lightweight,  low  power  and 
inexpensive  package.  Motorola's  proprietary  binary  Interface 
protocol  allows  the  user  great  flexibility  with  little 
sacrifice  of  convenience.  In  addition,  the  standardized  RS*>232 
serial  communications  port  enables  easy  access  to  a  variety  of 
different  systems.  All  things  considered.  It  Is  the  best 
choice  over  all  other  GPS  receivers  for  this  project. 

The  PVT-6's  reacquisition  times  compare  favorably  with  the 
manufacturer's  claims.  When  two  of  these  receivers  were  used 
In  a  DGPS  mode  with  a  zero  baseline,  the  results  easily 
confirmed  nominal  results  of  one  to  three  meters  desired 
accuracy.  The  most  encouraging  result  Is  the  relative  Immunity 
of  DGPS  accuracy  to  short  delays  (one  to  seven  seconds)  of  the 
pseudorange  correction  Inputs. 
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The  succ<’ss  of  this  project's  ultimate  goal;  to  perform  an 
autonomous  landing,  will  depend  on  improving  processing  speeds 
and  accuracies.  A  few  areas  where  further  investigation  into 
GPS  could  contribute  to  this  end  are  listed  below. 

•  Conduct  dynamic  tests  to  confirm  expected  nominal 
accuracies  using  DGPS. 

•  Conduct  additional  data  collection  to  determine  if  the 
accuracy  of  DGPS  or  reacquisition  time  is  affected  by  the 
time  of  day  or  by  the  weather. 

•  Collect  more  data  to  determine  if  psuedorange  correction 
latency  steadily  degrades  or  has  anomalies  (wherein  the 
distance  error  actually  improves  with  time)  as  Figure  39 
suggests  at  the  15  second  mark. 

•  Determine  the  effects  of  the  aircraft's  RF  environment  on 
the  GPS  performance. 

•  Investigate  whether  programming  in  assembly  or  some  other 
language  might  improve  processing  speeds. 

•  Determine  if  the  PVT-6's  capability  to  compute  two 
dimensional  position  by  constraining  altitude  (relying 
solely  on  barometric  altitude  sensor  data  for  altitude 
inputs)  might  Improve  robustness  of  GPS  satellite 
tracking. 

•  Write  a  serial  port  driver  that  is  independent  of  the  PCL** 
744  multiple  serial  port  I/O  card  in  order  to  provide 
flexibility  in  using  the  PVT-6  and  possibly  reduce  the 
total  hardware  weight  through  the  elimination  of  the  PCL- 
744  itself. 

•  Write  a  mission  initialization  interface  routine  to  setup 
up  the  PVT-6  receivers. 

•  Write  a  post-processing  routine  to  interpret  and  present 
the  recorded  results  in  a  desired  format. 

•  Write  a  real  time  graphical  interface  to  plot  the 
coordinates  on  a  map.  Also  include  either  a  separate 
window  or  color  coded  routine  to  be  able  to  see  current 
altitude. 
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GPS  has  revolutionized  electronic  navigation.  The  airline 
industry  and  other  nations  decry  the  U.S.  DOD's  introduction 
of  Selective  Availability.  They  are  actively  lobbying  the  U.S. 
government  to  force  the  DOD  to  turn  it  off  and  leave  it  off. 
It  is  an  uphill  battle  but  not  a  hopeless  one.  In  time,  as 
civilian  applications  become  common  place  in  our  society,  the 
DOD  may  acquiesce  and  turn  it  off.  If  that  happens,  GPS 
accuracies  achieved  in  this  project  will  undoubtedly  improve 
and  may  warrant  reconsideration  in  the  way  they  are  used  in 
the  Kalman  Filter  that  integrates  data  from  all  the  sources. 
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APPENDIX  A:  COORDINATE  TRANSFORMATIONS 

A.  GEODETIC  COORDINATE  TRANSFORMATION  TO  ECBF  COORDINATES 
Given  the  geodetic  coordinates  longitude  (X) ,  latitude  (tp) 
and  altitude  (h) ,  the  Earth  Centered  Earth  Fixed  coordinates 
(X,y,Z)  are  (Leick,  1990,  p.  184-lfiS) : 

X  •  (  N  +  A  )  cos<J>  cosX 
y  •  (  N  +  jh  )  cosd)  sinX 
Z  -  [N(l-e2)  *  h  ]  sin4> 

where  h  is  the  height  above  the  ellipsoid  and 

e^  -  2r  - 

f  -  ^  ^ 

a 

a  -  semi -major  axis 
b  -  semi-minor  axis 

N  -  - 

yTi  -  sin0 


B.  ECEF  COORDINATE  TRANSFORMATION  TO  GEODETIC  COORDINATES 
Given  the  ECEF  coordinates  (X, F, Z),  the  Geodetic 
coordinates  (X,0,h)  are  (Leick,  1990,  p.  184-185); 


88 


(p  -  arctan 


1 


1  - 


X 


=  arctan^ 


COB(p 


_  _ ^ _ 


-  N 
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APPENDIX  B:  EPHEMERI8  ALGORITHM 

The  following  parameters  are  provided  continuously  for 
each  satellite  of  the  GPS  constellation.  Their  exact  placement 
and  description  within  the  ephemeris  message  is  defined  by  the 
Interface  Control  Document,  ICD-GPS-200  as  published  by 
Rockwell  International  Space  Systems  Division.  Any  serious 
manipulation  of  the  ephemeris  algorithm  using  GPS  data 
requires  careful  consideration  of  this  document.  There  are 
several  important,  yet  subtle,  facts  about  the  time  tagging 
and  the  corresponding  validity  of  the  Issue  of  Data  messages. 
Mean  Anomaly  at  Reference  Time 
Mean  Motion  Difference  from  Computed  Value 
Eccentricity 

Square  Root  of  the  Semi-Major  Axis 

Longitude  of  Ascending  Node  of  Orbit  Plane  at 
Weekly  Epoch 

Inclination  Angle  at  Reference  Time 
Argument  of  Perigee 
Rate  of  Right  Ascension 
Rate  of  Inclination  Angle 

Amplitude  of  the  Cosine  Harmonic  Correction 
Term  to  the  Argument  of  Latitude 

Amplitude  of  the  Sine  Harmonic  Correction  Term 
to  the  Argument  of  Latitude 


Mo 

An 

e 

(A)*'* 

no 

io 

0) 

OMEGADOT 
I  DOT 
Cue 
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C„  Amplitude  of  the  Cosine  Harmonic  Correction 

Term  to  the  Orbit  Radius 

Cn,  Amplitude  of  the  Sine  Harmonic  Correction  Term 

to  the  Orbit  Radius 

Cjc  Amplitude  of  the  Cosine  Harmonic  Correction 

Term  to  the  Angle  of  Inclination 

Cj,  Amplitude  of  the  Cosine  Harmonic  Correction 

Term  to  the  Angle  of  Inclination 

t^  Reference  Time  Ephemeris 

lODE  Issue  of  Data  (Ephemeris) 

The  following  tables  list  the  number  of  bits,  the  location 
within  words  three  through  ten  (omitting  the  parity  bits) ,  the 
scale  factor  of  the  LSB,  the  range  and  the  units  for  each  of 
subframes  two  (Table  VII)  and  three  (Table  VIII) .  Parameters 
indicated  with  an  aster iB]c  are  to  be  two's  complimented  with 
the  sign  bit  in  the  MSB.  Unless  otherwise  listed  in  the  value 
range  column,  effective  range  is  the  maximum  range  attainable 
with  the  indicated  bit  allocation  and  scalr  factor. 

The  following  elements  of  the  ephemeris  algorithm  are 
listed  in  order  of  expected  solution.  All  the  required 
information  is  given  or  may  be  found  using  the  above  from  the 
GPS  ephemeris  message. 

2  WGS  84  value  of  the  earth's 

=  3.986005  X  10^*  meters  universal  gravitational 

sec^  parameter 
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Table  VII  BIT  ALLOCATION  AND  SCALING  FACTORS  FOR  WORDS  THREE 
THROUGH  TEN  OF  SUBFRAME  TWO.  {ROCKWELL,  1987,  PP ,  60-81) 


lODE 


No.  of 
bits 

Location 
of  bits 

8 

1-8 

16* 

9-24 

16* 

25-40 

mi 

41-72 

16* 

73-88 

32 

89-120 

16* 

121-136 

32 

137-168 

16 

169-184 

8 

185-192 

units 


see  ICD- 
GPS-200 


meters 


semi¬ 

circles/ 

sec 


semi¬ 

circles 


radians 


radians 


meters*'* 


604,784  seconds 


other 


^Para/netera  indicated  ar«  aignad  intagara.  All  othara  ara  unaignad 


see  ICD- 
GPS-200 


.  WGS  84  value  of  the  earth's 

0  =  7  2921151467  x  10’®  rotation  rate 

sec 


A  = 


JL 

sec 


tv  =  t  -  t. 


Semi-major  axis 
Computed  mean  motion 


Time  from  ephemeris  reference 
epoch 


In  the  above  equation,  t  is  GPS  system  time  at  time  of 
transmi  nn  (i.e,,  GPS  cotreoted  for  transit  timn 
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Table  VIII  BIT  ALLOCATION  AND  SCALING  FACTORS  FOR  WORDS 
THREE  THROUGH  TEN  OF  SUBFRAME  THREE.  (ROCKWELL,  1987,  PP. 
60-81) 


Parameter 

No.  of 
bits 

Bit 

location 

LSB  scale 
factor 

Value 

range 

Units 

Cic 

16* 

1-16 

2-29 

radians 

no 

32* 

17-48 

Si 

semi¬ 

circles 

16* 

49-64 

radians 

32* 

65-96 

m 

semi¬ 

circles 

CfC 

16* 

97-112 

2* 

meters 

Q 

32* 

113-144 

2-31 

semi¬ 

circles 

OMEGADOT 

24* 

145.-168 

2-43 

semi¬ 

circles/ 

sec 

lODE 

8 

169-176 

see  ICD- 
GPS-200 

IDOT 

14* 

177-190 

2-*» 

semi¬ 

circles/ 

sec 

other 

2 

191-192 

see  ICD- 
GPS-200 

^Parameter*  indicafcl  ara  aigntd  int»g»rM.  All  oth»rM  ara  unaign»d. 


range/speed  of  light.,  .  Furthermore,  shall  be  the  actual 
total  tine  difference  between  the  time  t  and  the  epoch  time 
t„,  and  must  account  for  beginning  or  end  of  wee)c  crossovers. 
That  is,  if  t^  is  greater  than  302,400  seconds,  subtract 
604,800  seconds  from  t,,.  If  t|(  is  less  than  302,400  seconds, 
add  604,800  seconds  to  t^. 


n  =  Hq  +  An 


Corrected  mean  motion 


Mk  “  Mo  +  nt^ 


Mean  anomaly 
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The  next  equation  solves  for  the  eccentric  anomaly  (EJ  by 
iteration.  Mu  and  e  are  known  from  above. 

Mk 

Let  E,  » 

'  1  -  e 

Then  M,  ■  -  e  sin  E  If  «  Mj,#  then  E2  ■  Ek 

otherwise  repeat  and 
Ii]  -  III  compare  again. 

^1  •  ^  -  e  coL  E  (Tattersfield,  1984,  p.  36) 

Ej  -  E,  +  A, 

Having  found  Ek,  the  rest  of  the  solution  is  simply  solving 
equations. 


Vk  ■  arctad 


VI  -  sinE> 
cosEk  -  e 


True  anomaly 


®  -  Vk  +  CJ 


Argument  of  latitude 


Second  harmonic  perturbations: 


fiUk*Cu,  sin2<Pk  +  Cue  cos24>k 

cos2<l>k  +  C„  sin2®k 
fiik=Cie  cos2<C>k  Ci,  sin2<l>k 

Ik  “  A  (1  -  e  cosEk)  +  fifk 
ik  =  io  +  6ik  +  (IDOT)  tk 


Argument  of  lat  correction 
Radius  correction 
Correction  to  inclination 

Corrected  argument  of  latitude 

Corrected  radius 

Corrected  inclination 
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Positions  in  orbital  plane 


Xk  *  Ijt  COSUjj 

yi  =  sinuj^ 

n  /rt  ^  t.  Corrected  ?  jnfjitude  of 

tijt  ■  Oq  +  (ti  -  n,)  tjj  -  H,t<„  ascending  node 


Xjt  =  Xk  cosO^  -  yi  cosijj  sinfljj 
Xjj  *  xi  sinflit  -  Vk  cosily  cosO,^ 
»  yi  sinik 


Earth-Centered , Earth-Fixed 
fECEF)  coordinates 
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APPENDIX  C:  PVT-6  SPECIFICATIONS 


Dynwnict  •  Vtlodly  1000  Knots  (SIS  m/tso) 

•  AooilartUon  4g 

•  Jsrk  Sm/s* 


Antotvw  •  Activ#  mierostrip  patch  Antanna  Modulo 

•  Poworod  by  Rooeivur  Moduli  (S  Vde  9  25  mA) 

Antenna  to  Rooolvor  •  Single  coaxial  oafalo  (6dB  max  lou  at  LI ;  167S.42  MHx) 

Intereonrwelion  •  Typteal  length  with  ROSS  coaxial  oablo;  20  loet  (6  maters) 


Sorial  Output  •  RS232C  Interface 

Accuracy  •  Less  than  28  meters,  SEP  (without  SA) 

*  DoO  may  invoka  Salaotlvo  Availability  (SA),  potontially  degrading 
accuracy  to  100  meters  (2d  RMS) 

Altitude  •  80,000  feel  (16  Km)  (maximum) 

Operating  Temperature  •  Aotiva  Antenna  -40*0  to  ^100'C  (•40*P  to  #21 2*P) 

•  Receiver  Module  -30*0  to  •ra0*O  (-22^  to  176*F) 

HumidHy  •  96%  non  «ondenting  4-30*0  to  -460*0  (•fSS*P  to  -r-l  40‘F) 

Physkal  Dimensions  •  ReosiverPCB  3.94  X  2.76  X  0.86  in  (100  X  70  X  16.6  mm) 

•  Antenna  Module  4.01  (dia)  X  0.89  In  (102  (dia)  X  22.6  mm) 

Weight  •  Receiver  and  Houelng  4.5  oz  (128  g) 

•  Antenna  Module  4,8  oz  (136.2  g) 

SwHehed  Power  •  9  to  18  Vdo  or 

•  8tp28Vde 

•Keep-AHvo*  BATT  Power  *  4.76  •  1 6  Vdo:  0.3  mA 

Power  Consumption  (typioal)  •  1.3  W  9  6  Vdo  Input 

•  1.6  W  9  12  Vdo  Input 

Rooeiver  Intereonneotlon  •  AMP  tO  pin,  #104369-4  on  Reooivsr  Module 

-  RS2320 

-  Power  Input 

•  08X  oonneetor  on  Receiver  and  Antenna 

MTBF  (Mean  Time  Before  Psilura)  •  $6,000  hours  (ogtlmaied) 
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APPENDIX  D:  PVT-6  ANTENNA  GAIN 


APPENDIX  E:  SOFTWARE  DRIVER  SOURCE  CODE 


A.  GPSAN.C 


/•  GPSAN.C 

Reads  ascii  position  data  "©©Ba"  from  a  file,  converts  it 
it  to  binary  and  places  the  scaled  and  formatted  messages 
in  an  array  of  structures. 

Input  is  the  name  of  the  ascii  file  without  the  extension. 

Output  is  a  file  with  the  analyzed  data  from  GPSCRNCH.H. 

V 

jfinclude  < string. h> 

#include  <stdlib.h> 
jflnclude  <stdio.h> 
jfinclude  "c  :\borlandc\twit6\gpstruct.  h " 

#lnclude  "c:\borlandc\twite\gpsconv. h* 

/^include  "e:\borlandc\twite\gpscrnch.h" 

^Include  "c:\borlandc\twite\gpsfun.h" 

void  main(void) 

{ 

#defina  COL  130 
#define  NEWCOL  66 

char  namel40I,  path(30)-''c!\\borlandc\\etwork\\"; 
int  i,j,k,m,n,s, master, slave,  ROW; 

ONEBYTE  first,  second,  third,  inbuffiCOLl; 

ONEBYTE  ••mbuff,  •*sbuff,  •*mbuffbin,  **sbuffbin; 
struct  T  POS  CHAN  STATUS  •norst,  ’sarst,  pos  chan; 

FILE  •fpna'ime;  “  ” 

printfr\nEnter  the  input  filename  (without  the  exten8ion(?7.ext):  "I; 

8canf|"%s",name): 

strcatipath,  name); 

strcatipath,  ".txt"); 

prlntf("%i",poth); 

fpname  •  fopen(path,"r"|; 

prlntf("\nEntar  the  number  of  rows  (Individual  messages):  "); 

8CBnf("%d",&ROW); 

printf("\n"); 
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k  =  0; 
m  =  0; 
s  -  0; 

r  Create  mbuff  and  sbuff. 

•/ 

mbuff  =  calloc(  ROW,  sizeof(  ONEBYTE  *  )|; 
sbuff  -  calloc(  ROW,  sizeofi  ONEBYTE  *  1); 
for  (  i-0;  i<ROW;  i+  +) 

{ 

mbuffli]  -  calloci  COL,  sizeofi  ONEBYTE  )); 
sbuffli]  -  calloci  COL,  sizeofi  ONEBYTE  )); 

); 


while  ID 

{ 

/*  Read  in  file  to  ascii  arrisys  mbuff  and  sbuff. 

•/ 

if  I  k  >  -  70  ) 
break; 

first  fgetclfpnameli 
if  I  first 
{ 

second  «  fgetclfpnamej; 
if  Isecond  >  «  '@'1 
{ 

third  «  fgetclfpname); 

/*  Store  inbuff  array  in  either  mbuff  or  sbuff. 

*/ 

if  Ithird  -  - 

{ 

fraadi  sbuffisl,  sizeoflONEBYTE),  COL,  fpname  I; 
k  -  0; 
s+  +; 

) 

else 

{ 

mbufflmllO)  -  third; 
for  |i-1;  i<COL;  i  +  +  ) 

mbufflmlli]  -  fgetclfpname); 
k  -  0; 
m  +  + ; 

); 

); 

); 

k+  + ; 

};  /*  end  while  •/ 
fcloselfpname); 

master  ■>  rn; 
slave  -  s: 
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/•  Create  mbuffbin  and  sbuffbin, 

*/ 

mbuffbin  ■  calloci  master,  sizeofi  ONEBYTE  *  )); 
for  (  i  =  0;  i<  master;  i  +  +) 

mbuffbinli)  =  calloci  NEWCOL,  sizeofi  ONEBYTE  )); 
sbuffbin  «  calloci  slave,  sizeofi  ONEBYTE  *  )); 
for  I  i  =  0;  i<  slave;  i  +  +  ) 

sbuffbinli]  -  calloci  NEWCOL,  sizeofi  ONEBYTE  )); 

r  Convert  hex  buffer  to  binary  buffer. 

•/ 

for  |i  »  0;  i  <  master;  i  +  +  ) 

{ 

mbuffblnlillOl  »  mbuffdKO]; 
mbuffbinlillll  -  mbuffliHll; 
k  >  2; 

forlj-  2;j  <  C0L;j+-2l 

{ 

mbuffbin(i][k]  «  |Hex|mbuff(i][jll  <<  4)  1  Hex|mbuff[il[j  + Hi; 
k  +  -  1; 

}; 

): 


for  II  ■>  0;  i  <  slave;  i  +  +  ) 

{ 

sbuffbinlillOI  -  sbufflillOl; 
sbuffbinlillll  -  sbuff[i][1]; 
k  »  2; 

for  I  j-  2;  j  <  COL;  j+  -21 

{ 

sbuffbin!l][kj  -  |Hex|sbuff(l][Jl)  <  <  4)  |  Hex|sbuff(i][j  + 1  J); 
k  +  -  1; 

); 

}: 


/•  Free  mbuff  and  sbuff  memory. 

*/ 

fur  li-0;  i<R0W;  i+  +  » 

{ 

freelmbufflil); 

freelsbufflil); 

}; 

freelmbuff); 

freelsbuff); 

/•  Temporarily  store  mbuffbin  in  a  file  "mtemp.bin"  and  free  memory. 

*/ 

fpname  -  fopeni  "mtemp.bin",  "wb"  ); 
for  |i  -0;  i<  master;  i  +  +  ) 

fwritelmbuffbinli],  sizeoflONEBYTE),  NEWCOL,  fpnamel; 
far  |i  -0;  i<  master;  i  +  +  ) 
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free(mbuffbinlil); 

free(mbuffbin); 

fclose(fpname); 

/*  Temporarily  store  sbuffbin  in  a  file  "stemp.bin*  and  free  memory, 

*/ 

fpname  «  fopani  "stemp.bin",  "wb"  ); 
for  (i«0;  i<slave;  i-t-  +) 

fwrite(sbuffbin(il,  sizeof(ONEBYTE),  NEWCOL,  fpname); 
for  (i-0;  i<slave;  i+  +) 
free(sbuffbinli]); 
free(sbuffbin); 
fclose(fpname); 

/*  Create  an  array  of  structures. 

*/ 

marst  ■  calloci  master,  sizeofi  struct  T_POS_CHAN_STATUS  )); 
sarst  ■  calloci  slave  ,  sizeofl  struct  T_POS_CHAN_&rATUS  )): 

fpname  «  fopeni  "mtemp.bin",  "rb"); 
for  (I  -  0;  I  <  master;  i  +  + ) 

{ 

freadlinbuff,  sizeoflONEBVTE),  NEWCOL,  fpname); 
Convert^Pos_Chan_Stati(s(  inbuff ); 
marstli]  *  pos  chan; 

}: 

fclosel  fpname  ); 

remove(*mtemp.bin*); 

fpname  -  fopeni  "stemp.bin",  "rb"); 

for  I  i  «  0;  I  <  slave;  i  +  + ) 

{ 

freadlinbuff,  sizeoflONEBVTE),  NEWCOL,  fpname); 
Convert_Pos_Chan_Status|  inbuff ); 
sarstli]  -  pos  chan; 

}; 

fclosel  fpname  ); 
remov6rstemp.bin"); 

Ave_Max_Stdev|marst,  sarst,  master); 

};  /*  end  main  */ 


101 


B.  GPSCONV.H 


/•  GPSCONV.H 

Convert_Ephemeris|) 
Con'/ert..Pos_Chan_Status() 
ConvertlxDOPO 
Convert__Sat_RanQe_Data() 
Convert  Pos  Chan  Status  ExtO 

*1 

fifndef  GPSCONV  H 
#define  _GPSCONV_H 

^include  <math.h> 

#ifnd8f  _GPSTRUCT_H 

finclude  "c:\borlandc\twite\gpstruct.h'' 

#endif 

. . . . . 


Convert_Ephemeris() 

This  procedure  converts  the  raw  binary  ephameris  data 
into  a  representation  usabie  for  post-processing 
applications.  The  equations  coded  in  this  procedure 
are  taken  from  the  ICD-GPS-200  (the  GPS  NAVSTAR  ICO). 

Input  is  an  array  with  the  binary  ephemeris  message. 
Although  it  is  unsealed,  it  is  in  the  array 
format  of  the  ICD-GPS-200  format. 

Output  is  a  structure  with  the  scaled  and  formatted 
ephemeris  message. 


void  Convert  Ephemeris!  ONEBYTE  hex  ephemerlsdOl  ) 

{ 

/*  Sign  extends  a  neg  hex  byte  into  its  TWOBYTE  (integer)  equivalent 

*1 

^define  MACRO_SET_UPPER(word_to..setl  \ 

if  T  word_to_set  &  0x0080  )  v/ord_to_set  |  ■  OxFFOO  ; 

/*  This  MACRO  takes  two  hex  bytes  and  merges  them  into  a  single 
TWOBYTE  (integer)  value. 

•/ 

Jddefine  MACRO  MAKE  EPH  WORO(e  wd,  w1,  b1,  w2.  b2)  \ 

e  wd  (1  unsigned"  TWOBYTE  )  hex_ephemeri8lw1 1(b1 1  <  <  8  I  |  \ 
(  UNSIGNED  TWOBYTE  )  hex  euhemerislw21[b2)  ); 
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UNSIGNED_TWOBYTE  word!,  word2;  /*  Temporary  storage  for  words  built.  */ 
struct  T_FLOAT_EPHEMERIS  ’eph  ; 

eph  =»  mallocisizeofi  struct  T_FLOAT_EPHEMERlSI  I; 

/*  Time  of  ephemeris 

•/ 

MACRO  MAKE  EPH  WORD(word1,  IB,  0,  15,  1  ); 
eph->toe  -  (doubleT  ( ((UNSIGNED_FOURBYTE)  wordi)  *  16  )  ; 

/*  Mean  anomaly  at  reference  time  (  mO  ). 

V 

MACRO  MAKE  EPH  WORDiwordI,  9,  2,  10,  0  ); 

MACRO  MAKE  EPH  WORD(word2,  10,  1.  10,  2  I; 

aph'>mO  -  (  n  FOURBYTE  )  wordi  <<  16  H  word2  )  •  PITWOM31: 

/*  Mean  motion  difference  from  computed  value  (  delta^n  ). 

V 

MACRO  MAKE  EPH  WORDiwordI,  9,  0,  9,  1  ); 
eph->delta_n  -  (  TWOBYTE  )  wordi  •  P1TWOM43; 

/*  Eccentricity  I  e  ). 

•/ 

MACRO  MAKE  EPH  WORDiwordI,  11,  2,  12,  0); 
MACRO”MAKE“EPH”wORD|word2,  12,  1,  12,  2  ); 

eph->e"-  I  1 1  UNSIGNED_FOURBYTE  )  wordi  <<  16  )  |  word2  )  *  TWOM33; 

/*  Square  root  of  the  semi^major  axis  I  aqrt.a  I. 

•/ 

MACRO  MAKE  EPH  WORDiwordI,  13.  2.  14,0); 

MACRO  MAKE  EPH  WORD|word2.  14,  1,  14,  2  ); 

eph->8qrt_a  -“l  I  I  UNSIGNED.FOURBYTE  )  wordi  <<  16  )  |  word2  )  *  TWOM19; 

/*  Longitude  of  ascending  mode  of  orbit  plane  at 
weekly  epoch  I  omega  0  ). 

♦/ 

MACRO  MAKE  EPH  WORDiwordI,  16,  2,  17,  0  ); 

MACRO  MAKE  EPH  WORD|word2,  17,  1,  17,  2  ); 

oph->omega_0"-  I  Fl  FOURBYTE  )  wordi  <<  16  )  |  word2  )  *  PITWOM31; 

/*  Inclination  angle  at  reference  time  I  iO  ). 

*/ 

MACRO  MAKE  EPH  WORDiwordI,  18,  2.  19,0); 
MACRO“MAKE"EPH"WORD|word2,  19,  1,  19,  2  ); 
eph->iO" -  I  I  TfOURBYTE  )  wordi  <  <  16  )  |  word2  )  *  PITWOM31; 

/•  Argument  of  perigee  I  w  ). 

•/ 

MACRO  MAKE  EPH  WORDiwordI,  20,  2,  21,  0  ); 

MACRO'mAKE  eph  WORD|word2,  21,  1,  21,  2  ); 

ephow"-  I  I  TfOURBYTE  I  wordi  <  <  16)1  word2  )  *  PITWOM31; 
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/•  Rate  of  right  ascension  (  omega  dot ). 

*! 

wordi  =  (  UNSIGNED_TWOBYTE  )  hex_ephemerisl221I01; 

MACRO  SET  UPPER(wordl) ;  /*  Set  upper  bits  if  negative. 

MACRO"maKE  EPH  WORD(word2,  22.  1,  22,  2  ); 

eph->omega_dot  -  (  <  I  FOURBYTE  )  wordi  <  <  16  )  |  word2  )  *  PITWOM43; 

r  Rate  of  inclination  angle  ( i  dot ). 

•/ 

MACRO_MAKEJPH_WORD( wordi,  23.  1.  23,  2  ) ; 
wordi  -  wordi  >  >  2  ; 

if  (  wordi  &  0x2000  )  /*  Set  upper  bits  if  negative.  */ 

wordi  I  >  OxCOOO; 

eph->l_dot  -  (  TWOBYTE  )  wordi  *  PITWOM43; 

/*  Amplitude  of  the  cosine  harmonic  correction  term  to  the 
argument  of  latitude  (cun). 

•/ 

MACRO  MAKE  EPH  WORDIwordI,  11,0,  11.  1  ); 
eph->cuc  -  I  fwOBYTE  )  wordi  *  TWOM29; 

/*  Amplitude  of  the  sine  harmonic  correction  term  to  the 
argument  of  latitude  (ous). 

•/ 

MACRO  MAKE  EPH  WORDIwordI,  13,  0,  13,  1  ); 
eph->cus  -  (  fwOB^YTE  )  wordi  *  TWOM29; 

/*  Amplitude  of  the  cosine  harmonic  correction  term  to  the 
orbit  radius  (crc). 

*/ 

MACRO  MAKE  EPH  WORD(word1,  20,  0,  20,  1  ); 
eph->crc  -  I  TWOBYTE  )  wordi  *  TWOMO; 

I*  Amplitude  of  the  sine  liarmonic  correction  term  to  the 
orbit  radius  (crs). 

•/ 

MACRO  MAKE  EPH  WORDIwordI,  B,  1,  8,  2  ); 
oph->cr8  -  (  TWOBYTE  )  wordi  *  TWOM5: 

r  Amplitude  of  the  cosine  harmonic  correction  term  to  the 
angle  of  iitclination  (cic). 

*i 

MACRO  MAKE  EPH  WORDIwordI,  16,0,  16.  1  ); 
eph‘>cic  -  I  TWOBYTE  )  wordi  *  TWOM29; 

/•  Amplitude  of  the  sine  harmonic  correction  term  to  the 
angle  of  inclination  Ids). 

•/ 

MACRO  MAKE  EPH  WORDIwordI,  18,  0,  18,  1  ); 
eph->cis  -  I  TWOBYTE  )  wordi  *  TWOM29; 
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/•  Estimated  group  delay  differential  ( tgd  1. 

*/ 

wordl  -  (  UNSIGNED_TWOBYTE  )  hax_ephemeris(41(21; 

MACRO  SET  UPPERIwordI  1  ;  /*  Set  upper  bits  if  negative.  */ 
eph->tgd  -  TtwOBYTE  )  wordl  *  TWOM31; 

/*  Clock  data  reference  time  ( toe  ). 

*/ 

MACRO_MAKEJPH^WORD(word1,  5,  1,  6,  2  ); 
aph>>toc  «  wordl  *  16.0; 

/•  Polynomial  coefficient  (  af2  ). 

•/ 

wordl  -  (  UNSIGNED  TWOBYTE  )  hex  ephemerlt[6]I0]; 

MACRO  SET  UPPERIwordI ) ;  /*  Set  upper  bits  if  negative.  */ 
eph->af2  -  TtwOBYTE  1  wordl  •  TWOM56; 

r  Polynomial  coefficient  (  afl  /. 

*/ 

MACRO  MAKE  EPH  WORDIwordI,  6,  1,  6.  2  ); 
eph->afl  -  (  TWOBYTE  )  wordl  •  TWOM43; 

/*  Polynomial  coefficient  I  afO  ). 

*/ 

wordl  -  (  UNSIGNED_TWOBYTE  )  hex^ephemerlsITKOl; 

MACRO,SET_.UPPER(word1) ;  /•  Sat  upper  bits  If  negative.  */ 

MACRO.MAKE_EPH_WORDIwofd2,  7,  1,  7,  2  )  ; 

word2  -  word2  &  OxFFFC  ;  /•  clear  2  bits  */ 

eph->afO  -III  FOURBYTE  )  wordl  <<1611  word2  )  •  TWOM3126; 

/*  Semi-major  axis  (  a  ). 

•/ 

eph->a  -  eph->sqrt_a  *  eph->eqrt_e; 

/*  Corrected  meirn  motion  (  n  |, 

•I 

eph->n  -  SQRMU  / 1  eph->a  *  eph->8qrt_a  )  +  eph->delta_n; 

t*  Square  root  of  eccentricitv*2  subtracted  from  1  t  esqrt  I. 

*/ 

eph->e_8qrt  -  sqrt  I  1 .0  •  eph->e  *  eph->8  ); 

r  Omega  dot  -  WE  (  wkl  i. 

•/ 

eph->wk1  -  eph->omegB_dot  -  WE; 

/*  OmegaO  •  WE  *  toe  (  wkO  ). 

•/ 

eph->wkO  •«  eph’>omega_0  -  WE  *  eph->toe; 

/•  Sine  w  I  sin  w  ). 
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*/ 

eph'>8in_w  -  sin  (  eph->w  ); 

r  Cosine  w  (  cos  w  ). 

•/ 

eph->C08_w  -  cos  (  eph->w  ); 

r  Square  root  of  semi-major  axis  *  eccentricity  *  FCONST. 

*/ 

eph- > sqrta^ejconst  ■  eph->8qrt_a  ♦  eph->e  *  FCONST; 

/*  Issue  of  data,  ephomeris. 

•/ 

eph->iode  -  hex_ephemeriB[23][0]; 

return  *eph; 

}; 

h* . . . * . * . 


Convert,  Po8_Chan_Statu80 

This  procedure  parses  and  decodes  a  posItion/status/data 
record  received  from  the  Motorola  PVT-S.  The  flow  of  this 
procedure  is  as  follows  : 

•  pass  in  the  BINARY  buffer  containing  66  bytes 

•  read  the  BINARY  data  from  the  buffer  and  convert  directly 
into  the  final  data  structure. 

Input  is  a  buffer  of  type  ONEBYTE  with  a  single  binary 
position/status/data  message. 

Output  is  the  structure  with  the  scaled  and  formatted 
position/status/data  message. 


struct  T  POS  CHAN  STATUS  Convert  Pos  Chan  StatuslONEBYTE  pos[]) 
{  ■  "  " 

struct  T  POS  CHAN  STATUS  pos  chan; 

UNSIGNED  ONEBYTE" I,  j  ; 

UNSIQNED"0NEBYTE  tempchar ; 

UNSIGNED  FOURBYTE  tempu4bvte  ; 

FOURBYTE  tempsAbyte  ; 

UNSIGNED_F0URBYTE  nsecs  ; 
double  degrees,  minutes,  seconds  ; 

/*  Read  and  scale  the  of  the  data. 

*/ 

pos_chan. month  >  pos|2]  ; 
pos_chan.day  -  posOI ; 

tempchar  ~  posl4]  ; 
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pos  chan.vear  =  (  tempchar  <  <  8  )  +  pos(5)  ; 

pos_chan. hours  -  posISl ; 
pos_chan.rrinutes  -  pos[7] ; 

tempchar  ■  posI8] ;  /•  integer  seconds  */ 
tempu4bvto  «  posid]  ; 

tempu4byte  «  ( tempu4byte  <<8)4^  pos[10] ; 
tempu4byte  >  ( tempu4byte  <<8)4  posd  1] ; 
tempu4byte  >  |  tempu4byte  <<8)4  posd  2] ; 
po8_chan. seconds  >  (double)  tempchar  4  \ 

(  ( (double)  tempu4byte  )  /  1,0E4  9  ); 

temps4byte  -  posd  3] : 

temp84byte  ( temp84byte  <<8)4  posd  4] ; 
temps4byte  ■  ( temps4byte  <<8)4  posd  5] ; 
tsmps4byte  •  ( temps4byte  <<8)4  posd  6] ; 
degrees  -  (double)  temps4byte  *  MSECS.TO.DEGREES  ; 

poa.ohan.latitude.degrees  -  (TWOBYTE)  degrees  ; 

if  (  degrees  <  0  ) 

degrees  >  fsbs  (  degrees  ) ; 

minutes  -  (  degrees  ■  (TWOBYTE)  degrees  )  ♦  60.0  ; 

pos_chan.latitude.mlnute8  ■  (TWOBYTE)  (  minutes  ) ; 

poi_chan.latitude.second8  -  (  minutes  •  (TWOBYTE)  minutes  )  *  60.0  ; 

temp84byte  ■  posd  7] ; 

temps4byte  •  ( temps4byte  <<8)4  posd  8] : 
temp84byte  -  ( temp84byte  <<8)4  posd  9) ; 
temps4byte  >  ( temps4byte  <<8)4  pos[20] ; 
degrees  -  (double)  temps4byte  *  MSECS.TO.DEGREES  ; 

pos.chan. longitude. degrees  ■  (TWOBYTE)  degrees  ; 

if  (  degrees  <  0  ) 

degrees  ■  fabs  (  degrees  ) ; 

minutes  -  (  degrees  ■  (TWOBYTE)  degrees  )  *  60.0  ; 

po8_chan. longitude. minutes  -  (TWOBYTE)  (  minutes  ) ; 

pos.chan. longitude. seconds  ■  (  minutes  •  (TWOBYTE)  minutes  )  *  60.0  ; 

temps4byte  ■  pos[21] ; 

temps4byte  -  ( temp84byte  <<8)4  po8(22) ; 
temps4byto  -  ( temps4byte  <<8)4  po8(23) ; 
temp84byte  -  ( temps4byte  <<8)4  pos[24) ; 
pos.chan.datum_heioht  -  (double)  temp84byte  / 100.0  ; 

temps4bYte  -  posl25]  ; 

temp84byte  «  ( temp84byte  <<8)4  po8(261  ; 
temp84byte  «  ( temp84bvte  <<8)4  po8(27)  ; 
temps4byte  -  ( tamp84bvte  <<8)4  po8(28)  ; 
pos.chan. msl.height  »  (double)  temp84byte  / 100.0  ; 
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tempchar  =  pos[29|  ; 

pos_cnan.  velocity  =  (double)  (  ( tempchar  <  <  8  )  +  posOOl  )  /  100,0  ; 


tempchar  -  pos[31 1  ; 

pos_chan. heading  s  (double)  ( ( tempchar  <  <  8  )  +  pos[32] )  / 10.0  ; 
tempchar  -  pos[33l ; 

poi_chan.current_dop  (double)  ( ( tempchar  <<  8  )  -t-  pos[34] )  / 10.0; 

pua_chan.dop_type  ■  posOB] ; 
pos_chan.vislble_sats  -  pos[36] ; 
pos^chan.satsjracked  *  pos(371 ; 

J-0; 

for  (i  «  0;  i  <  NUM  CHANNELS;  l-f  •»•) 

{ 

poa^chan.channaKll.svId  -  pos[38  -t-  j); 
poa^chan.channalllj.moda  poa[39  -f  j); 
pos_chan.channel[l],atrangth  »  po8[40  +  jl; 
poa  chan.channeKD.flags  ■  po8[41  +  j); 

J+-4) 

) 

p08_chan,rcvr_8tatUB  ->  pos|62]; 

return  po8  chan; 

}) 

#endlf  /*  .QPSCONV^H  •/ 
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C.  GPSCRMCH.H 


/*  GPSCRNCH.H 

Used  to  crunch  the  position  data  to  determine  the  mean, 
maximum,  and  standard  deviation  of  the  position  error 
obtained  by  DGPS.  Geodetic  coordinates  are  converted  to 
ECEP  coordinates  and  the  distance  formula  is  used  to  compute 
the  distance  error. 

Input  is  an  array  of  structures  that  contain  the  scaled 
and  formatted  poaition/chan/statue  message  data. 

The  user  is  asked  to  supply  the  output  file  name 
without  the  extension. 

Output  is  a  file  with  the  distance  error  for  each  data 
point,  and  the  mean,  max  and  standard  deviation  of 
all  the  data  points. 


#ifndef  GPSCRNCH  H 
#deflne  _GPSCRNCH_H 

#inoludo  <math.h> 
finclude  <stdio.h> 
finclude  <  string. h> 

#lfndef  _GPSTRUCT.H 

#lnclude  *c!\borlandc\twlte\gp8truct.h" 

#endif 

FILE  ‘fp; 

void  Ave  Max  Stdevistruct  T  POS  CHAN  STATUS  ‘mast,  \ 

struct  T  POS  CHAN  STATUS  *8lav,  int  master) 

( 

Int  I,  count  -  0,  start; 

char  nam8l401,  outpathOO]  ■"ciWborlandcWtwiteW"; 

double  mlat,  mlon,  malt,  slat,  slon,  salt,  mN,  sN,  mean,  dlstsum-0; 

double  mx,  my,  mz,  sx,  ey,  sz,  dist[200],  max»-1.0,  stdev,  distsum2>«0i 

prIntfrEnter  the  output  filename  without  the  extension  (777.txt);  "I; 

scanfi*%s",nBme); 

strcatloutpath,  name); 

strcatloutpath,  ".out"); 

fp  -  fopen(outpath,"w''); 

fprlntflfp,''\n%8\n",outpath); 

prlntf(”Cnter  the  starting  data  point  number;  "); 
scBnfr%d",&start); 
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for  (i  =  start;  i  <  master;  i  +  +  ) 

{ 

/*  Test  to  ensure  the  seconds  match  */ 

If  ( (ONEBYTE)  mastdl. seconds  >•  «  (ONEBYTE)  slav[i|. seconds  ) 

{ 

/•  Convert  to  mx,  my,  ms,  sx,  sy,  sz. 

•/ 

mlat  (mastliMatitude. seconds  *  (1/60)  +  \ 

mastdl. latitude. minutes)  *  (1/60)  ■<' mastdl. latitude. degrees; 
mion  (mastdl. longitude.seconds  *  (1/60)  -f  \ 

mastdl. longitude. minutes)  *  (1/60)  -f  mastdl. longitude. degrees; 
malt  mastdl. datum^helght; 

slat  ■  (alevdl. latitude. seconds  *  (1/60)  +  \ 

tlavdl. latitude. minutes)  *  (1/60)  4 slavdl.latltude. degrees; 
slon  >  (slavdl. longitude.seconds  *  (1/60)  4  \ 

slavdl. longitude. minutes)  *  (1/60)  4  alavdl.longitude. degrees; 
salt  •«  slavdl. datum_height; 

mN  »  SEMIMAJOR  AXIS  /  sqrt(  1  •  E  SQUARED  *  sinlmlat)  ): 
sN  -  SEMIMAJOR.AXIS  /  aqrt(  1  •  E.SQUARED  *  sin(slat) ); 

mx  «  (  mN  4  malt )  •  coslmlat)  *  cos(mlon); 
my  «  (  mN  4  malt )  *  cos(mlst)  *  sln(mlon); 
mz  ■«  (  mN  •  (  1  •  E_SQUARED  )  4  malt)  *  sln(mlat); 

sx  «  (  sN  4  salt )  *  cos(slat)  *  oosislon); 
sy  •  (  sN  4  salt )  *  cosislat)  *  sinlslon); 
sz  -  (  sN  *  (  1  •  E.SQUARGD  )  4  salt)  *  sin(slat); 

distdl  ■  8qrt(  (mx-sx)*(mx-sx)  4  (my8y)*(my8y)  4  (mz*,sz)*(mz*8z) ); 

If  (  distdl  >  1000  ) 

{ 

fprlntflfp."\nMtime  -  %2d;%10.7f  <->  %2d:%10.7f  -  Stime  \ 
Dist  #%3d  »  XXXXX",m8stdl. minutes,  mastli). seconds,  \ 
slavdl. minutes,  slavdl. seconds,  I); 
continue; 

}; 


/*  Compute  average,  max  and  standard  deviation. 

*1 

count  4  4 ; 

if  (  dlst(il  >  max  )  max  -  distdl; 
distsum  4  «  distdl; 
diBtsum2  4  o  distil]  *  distdl; 

/*  Print  out  the  results  to  a  file. 

*/ 

fprlntf(fp,*\nMtlmo  -  %2d;%10,7f  <  •>  %2d:%10.7f  «  Stime  \ 
Dist  ^%3d  >  %6  2^, mastdl. minutes,  mastdl. seconds,  \ 
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slavli], minutes,  slav[i|. seconds,  i,  distdl); 


else 

{ 

fprintf(fp,"\nMtime  -  %2d:%10.7f  <->  %2d:%10.7f  -  Stime  \ 
Oist  #%3d  «•  XXXXX*,mast[ll. minutes,  mastli]. seconds,  \ 
slavti], minutes,  slavli). seconds,  i); 

); 

}i 


/*  Compute  the  mean  and  standard  deviation. 

V 

mean  >  distsum  /  count; 

stdev  ■  iqrti  (di8tsum2/oount)  •  (mean  *  mean)); 

fprintfifp,  ''\n\nMaan  -  %B.2f\nMax  -  %B.2f\nStdBv  -  %5.2f\n",\ 
mean,  max,  stdev); 

fcloselfp); 

);  /*  end  Ave_Max_Stdev  •/ 

#endlf  /*  _QPSCRNCH,H  */ 
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D 


6PSDEFIN.H 


/•  GPSDEFIN.H 

Store  most  of  the  definitions  here. 

•/ 

#ifndef  GPSDEFIN  H 
fdefine". GPSDEFIN.H 

/****•**••  type  definitions  ♦***•♦**•/ 

typedef  char  ONEBYTE;  /*  data  type  allowing  for  the  use  of  8  bits.  *1 
typedef  short  TWOBYTE;  I*  data  type  allowing  for  the  use  of  1 6  bits.  */ 
typedef  long  FOURBYTE;  /*  data  type  allowing  for  the  use  of  32  bits.  */ 

typedef  unsigned  char  UNSIGNED.ONEBYTE;  /*  unsigned  8  bits  */ 
typedef  unsigned  short  UNSIGNED  TVyOBYTE;  /*  un&iigned  16  bits  */ 

typedef  unsigned  long  UNSIGNED>OURBYTE;  I*  unsigned  32  bits  */ 


/****••**•  constant  definitions  ♦♦♦••/ 

j|«definaPI  3.1416926535898  /*  GPS  value  of  PI  V 


#deflna  WE  7.2921 161 467E-5  /*  WGS84  earth's  rot  rate  (rads/seci  */ 

#definaSQRMU  1 .9964981 84321 74E7  /*  Square  root  of  MU.  */ 

^define  FCONST  •4.442807633E-1 0  /«  FCONST  is  used  in  relativistic 

clock  correction.  Defined  on 
page  73  of  ICD-GPS-200.  */ 

Idefine  INVERSE  FLATTENING  298.267223563  /*  WGS-84  unitless  */ 

#define  FLATTENING  1 /INVERSE  FLATTENING 
#define  SEMIMAJOR  AXIS  63781 37.0‘  /*  WGS-84  meters  V 

#dafine  E.  SQUARED  2*FLATTENING-FLATTENING*FLATTENING 

fl'define  TWOM5  0.03125000000000  /*  2*-5  */ 

#dafineTWOM11  (TWOM5  /  64, Ui 
Idefine  TWOM19  (TWOM1 1  /  256.01 
^define  TWOM20  (TWOMI9/  2.01 
#define  TWOM21  (TWOM20  /  2,01 

#define  TWOM23  (TWOM21  /  4.01 

#define  TWOM24  ITWOM23  /  2.01 

#deftne  TWOM27  {TWOM24  /  8.01 

Idefine  TWOM29  (TWOM23  /  64,0) 

Idefine  TWOM30  (TWOM29  /  2.0) 

Idefine  TWOM31  (TWOM30  /  2.0) 

Idefine  TWOM33  (TWOM31  /  4,0) 

Idefine  TWOM38  (TWOM33  /  32.0) 

Idefine  TWOM43  (TWOM38  /  32.0) 
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(((define  TWOM50  (TW0M43  /  128,0) 
(((define  TWOM55  (TWOM50  /  32.0) 
#define  TWOM3126  (0.25  *  TWOM31) 


#define  TWOP1 1  2048,0  /*  2"11  */ 

#define  TWOP1 2  (TW0P1 1  *  2.0) 

(((define  TWOP1 4  (TW0P1 2  *  4.0) 

(((define  TWOPl  6  (TWOPl  4  *  4.0) 

(((define  TW0P1 2  INTEGER  4096  /•  2*12  •/ 

(((define  TWOPl  6JNTEGER  6B536  /*  2*ie  */ 

(((define  PITW0M31  (PI  *  TW0M31) 

(((define  PITWOM43  (Pi  *  TWOM43) 

/*********  Twite  definitions  **•**•••/ 

(((define  FALSE  0 

#defins  TRUE  1 

#define  RX  0  /*  receive  */ 

#define  TX  1  /*  send  */ 

#deflne  MGPS.  PORT  3  /*  master  gps  receiver  port  */ 

(Kdefine  SGPS'PORT  4  /*  slave  gps  receiver  port  •/ 

(Kdefine  NUM  CHANNELS  6 

(((define  MSECS  JO.DEGREES  ((1.0/  1 000.0  I  /  3600.0  ) 

#define  EPH.NUM_WORDS  24  /•  at  3  bytes/word  */ 

(Kendlf  /*  ^GPSDEFIN  H  */ 
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E.  aPSDELAY.H 


/*  GPSDELAY.H 

Used  to  introduce  an  artificial  delay  in  transmitting 
psuedorange  corrections  from  the  master  to  the  slave. 

PR_delay()  accepts  a  pointer  of  type  PANDL  to  the  buffer 
containing  the  psuedorange  correction  message 
(®<S>Ce...)  and  an  integer  number  of  seconds  it 
is  to  be  delayed.  The  output  is  to  the  PCL-744 
serial  I/O  card. 

V 

#lfndef  GPSDELAY  H 
#deflne  _GPSDELAY_H 

#lnclude  <stdlo.h> 
linclude  <  string. h> 

^Include  <stdl!b.h> 

#lnclude  <time.h> 
linclude  <do8.h> 

^Include  "c:\borlandc\twite\head-c.h* 

#ifndef  _GPSTRUCT_H 

#lnclude  "c:\borlandc\twlte\gpstruct.h" 

#endlf 

#deflne  MAXDELAY  300 

static  PANDL  storlMAXDELAYl;  /*  Ring  buffer  •/ 

static  Int  num_waitlng  ••  0;  /*  Number  of  messages  In  queue  * 

static  Int  next,8tor  -  0;  /*  Next  storage  location  to  be  used  */ 

static  int  next_xmit  ■  0;  /*  Next  message  to  be  transmitted  */ 

void  PR  delay!  PANDL  *Ce  msg,  int  delay  I 

{ 

PANDL  ’ptr; 

/•  Copy  the  contents  of  the  buffer  with  the  Ce_m5g  to  a 
new  buffer  for  holding  during  the  delay. 

•/ 

ptr  -  malloci  sizeoflPANDL)  |; 

ptr->len  ■  Ce_msg->len; 

ptr- > ptr  -  calloci  ptr->len,  sizeoflONEBYTE)  ); 

memcpyl  ptr-'* ptr,  Ce_msg->ptr,  Ce_msg->len  ); 

/*  Store  the  structure  of  type  PANDL  in  the  ring  b"ffer  "storll". 

*/ 

If  (num_waiting  <  MAXDELAY  I 
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{ 

stor(next_stor]  =  *ptr; 
nunn_  waiting  +  +  ; 
next_stor  +  +  ; 

if  (next  stor  «  »  MAXDELAY) 
next  stor--  MAXDELAY; 

} 

else 

{ 

printf(*There  are  %d  \”®(SiCe..A”  messages  in  the  ring  buffer.\n",\ 
MAXDELAY); 

printfCOne  correction  message  was  lost.*); 

froe(ptr->ptr); 

free(ptr); 

return; 

} 

/*  Transmit  the  oldest  message  if  there  ara  more  messages 
than  the  delay  length.  The  delay  is  premised  upon  the 
fact  that  correction  messages  come  approximately  one  second 
apart. 

•/ 

while  (  (num  waiting  -  delay)  >  0  ) 

{ 

sio_putb(  SGPS_PORT,  storlnext^xmitl.ptr,  stor[next„xmitl,len  ); 

free(stor[next_xmitl.ptr); 

next  xmit  +  -f  ; 

if  (next  xmit  -  »  MAXDELAY) 
next.xmlt--  MAXDELAY; 
num  waiting-  ; 

): 


free(ptr); 

return; 


}: 


#endi*  r  _GPSDELAY_H  */ 
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p 


GPSFUN.H 


/•  GPSFUN.H 

A  collection  of  functions  to  manipulate  GPS  messages. 
HexO 

Check_sumO 

Bin_to_asciiO 

Mso_verification() 

ParseJnbuffO 
Master_aps(l 
Slave  gpsO 

*1 

X^ifndef  GPSFUN  H 
#define  _GPSFUN_H 

^include  <stdlib.h> 

^include  < string. h> 

/^include  <stdio.h> 

^include  ''c:\borlandc\twite\head-c.h" 

Jdlfndef  _GPSTRUCT_H 

#include  ‘*c:\borlandcUwite\gp8truct.h” 

#endif 

FILE  •m_all_flle;  /*  master  file  for  all  data  •/ 

FILE  ♦m_p08jlle;  /•  master  file  for  pos  data  •/ 

/*  Global  matrix  for  converting  ephemeris 

•/ 

UNSIGNED_ONEBYTE  hex_ephemeri8l24K31 ; 

/*  (wordllbytel  :  24  words  (3  bytes  per  word)  corresponding  to 
words  3*10  of  subframes  1'3  of  the  satellite  nav  message. 

•/ 


/ 


This  procedure  converts  an  ASCII  value  to  its  hex 
equivalent. 

Input  is  a  single  ONEBYTE. 

Output  is  a  single  ONEBYTE. 
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ONEBYTE  HexiONEBYTE  ASCIIvaiue) 

{ 

/•  0  is  0x30(48) ->  9  is  0x39(57) 

A  is  0x41(65) ->  Pis  0x46(70) 
a  is  0x61(97)  •  >  f  is  0x66(102) 

•/ 

if  ((ASCIIvaiue)  <  0x40) 
return((ASCIIvalue)  -  0x30): 

else 

if  {(ASCIIvaiue)  <  0x60) 
return(IASCIivalue)  ■  0x37); 
else 

return((ASCIIvalue)  •  0x57); 

}; 


/ 


Check_8um() 

Calculates  the  checitsum  of  a  massage  In  a  buffer  of  type 
DNEBYTE.  Since  the  message  is  assumed  to  be  in  the  Motorola 
Binary  Proprietary  format,  all  but  the  last  three  bytes 
are  used  in  computing  the  checicaum. 

Input  is  the  pointer  of  type  PANDL. 

Output  is  the  checksum  of  type  ONEBYTE. 


ONEBYTE  Check  8um(  PANDL  *msg  ) 

{ 

ONEBYTE  chk^sum  -  0; 
int  i; 

I*  Compute  the  Exclusive*Or  of  all  but  the  last  three  bytes 
V 

for  ( i  ■  0;  I  <  (m£g">l8n  -  3);  i  +  +  ) 
chk_sum*-  msg->ptr[il; 
return  ohk  sum; 

}: 


/ 


Binjo.asciiO 

Converts  a  message  from  binary  to  its  hexidecimal 
ASCII  equivalent.  All  bytes  of  the  input  message  are 
converted  to  ASCII  except  the  first  n  values  as  specified 
in  the  integer  input  argument.  (This  was  incorporated  to 
permit  the  use  of  this  function  with  the  Motorola  GPS 
proprietary  message  format  where  the  first  four  bytes  are 
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already  in  ASCII.) 

Input  is  a  pointer  of  type  PANDL  and  the  interger  value 
of  the  number  of  bytes  to  skip  before  conversion. 
Output  is  a  pointer  of  type  PANDL  to  the  new  buffer 
with  the  ASCII  converted  massage. 


PANDL  *Bin  to  ascii!  PANDL  *bin  msg,  int  skip  ) 

int  i,  j.  index,  tmplen  *  0  ; 

ONEBYTE  nibbie(21,  •ptrjool; 

PANDL  •ascii^msg; 

tmplen  =  bin_mso->len  *  2  -  skip; 
ptrjooi  ■  bin^msg-^ptr; 

/*  Create  new  memory  for  the  converted  message. 

*/ 

a8cii_m8a  »  malioo!  slzeofi  PANDL  )); 
aBCii_msa>>ptr  -  calioc!  tmplen.  sizeoffONEBYTE) ); 
a8Cii_msa->len  «  tmplen; 

/*  Copy  any  initial  values  to  the  new  buffer  that  don’t  need  conversion. 

*1 

memcpy(ascii_m8g*>ptr, ptrjool, skip); 

for  ( I  ■  skip;  I  <  tmplen;  i  +  +  I 

{ 

/*  Split  the  byte  into  its  upper  and  lower  nibbles. 

•/ 

nibble[0]  ■  (  ptrjooKil  >  >  4  )  &  OxOP;  /*  upper  •/ 
nibble!!]  -  !  ptrjooKi)  &  OxOF  );  /•  lower  */ 

for  ( j  -  0;  j  <  2;  j  +  +  ) 

{ 

/*  Convert  the  nibbles  to  ASCII. 

•/ 

if  I  nibble!])  <  10  ) 

nibble!]]  +  -  0x30;  /•  nibble +  0x30 -0x30:0x39  */ 

else 

nibble!]]  +  -  0x37;  /•  nibble +  0x37 -0x41  ;0x46  •/ 

): 

/*  Compute  the  index  and  insert  into  the  ascii  msg  buffer 

*/ 

index  -  2  *  i  •  skip; 
ascii_msg->ptr!  index  ]  -  nibblelOl; 
ascii  msg‘>ptr!  index  +  1]  -  nibble!!]; 

): 

return  ascii_msg; 
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Msg.verificationO 


Specifically  used  to  verify  the  Integrity  of  the 

Motorola  GPS  proprietary  message  format. 

Input  are  a  pointer  of  type  PANDL  and  an  integer  designating 
whether  the  message  was  received  from  the  GPS  receiver  (Rx) 
or  is  to  be  transmitted  to  the  GPS  receiver  (Tx). 

Output  is  the  integer  flag  indicating  whether  or  not  the 
message  passed  verification.  True  means  it  passed. 


int  Msg  verificationi  PANDL  *msg,  int  msg  direction  ) 

{ 

int  I,  len  ■>  *1 ,  array.col; 
int  status  ■  TRUE;  ~ 


ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 

ONEBYTE 


bytaO  «  m8g*>ptr[0];  /*  Msg  prefix  @  byte  *! 

bytel  “  mag‘>ptr[li;  /*  Msg  prefix  ®  byte*/ 

byte2  -  m8g->ptr(21;  /*  Msg  id  first  byte*/ 

byteS  -  msg‘>ptr(3];  /*  Msg  id  second  byte  */ 

byteCS  -  m8g->ptr[  m80'>len-3  1;  /*  Check  Sum  byte  */ 
byteCR  •  mig->ptrl  mso->len-2  1;  /*  Carr  Retn  byte  */ 
byteLF  «  mag*>ptr(  rnsg*>len>1  );  /*  Line  feed  byte  */ 


/*  Determine  if  message  verification  is  for  Tx  or  Rx  because  the 
number  of  byte  in  the  message  is  different  for  a  Tx  message 
than  for  an  Rx  measage. 

•/ 

if  (  msg^direction  ) 

array^co!  »  0;  /*  Tx  message:  to  GPS  receiver  */ 

else 

array_col  •  3;  I*  Rx  message;  fm  GPS  receiver  */ 

/*  Message  ID  existence  verification  and  expected  length. 

*/ 

if  (  byte2  -  -  'A'  ) 

for  ( i  0;  I  <  30;  i  +  +  ) 

{ 

if  (  byte3  -  -  A  arraylillarray  col] ) 

{ 

if  (  msg  direction  ) 

len  »  A_arrav(i][1);  /*  Msg  size  to  PVT-6  */ 

else 

len  «  A_arravli][2];  /*  Msg  size  fm  PVT-6  */ 
break; 

); 
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}: 

if  I  byte2  -  «  'B'  ) 

for  ( i  -  0;  i  <  1 1 ;  i  +  +  i 

{ 

If  (  byte3  ■>  -  B  arrayli] [array  col) ) 

( 

if  I  msg  direction  ) 

len  -  B.arraylild  1.  r  Msg  size  to  PVT-6  */ 

else 

len  -  B.array[il[2];  /*  Msg  size  fm  PVT-6  */ 
break; 

}; 

}; 

If  (  byte2  -  -  'C'  ) 

for  ( I  -  0;  I  <  11;  1+  +  ) 

{ 

if  (  byteS  ■  C  array  [I]  (array  cot] ) 

{ 

if  (  msg_direction  ) 

len  -  C.arraydlll);  /*  Mag  size  to  PVT-6  */ 

else 

len  -  C„array(l][2];  /•  Mag  size  fm  PVT-6  •/ 
break; 

}i 

)l 

if  (  byte2  I  -  'A'  &&  byte2  I  -  'B'  &&  byte2  I  -  'C*  I 

{ 

status  «  FALSE; 
return  status; 

}; 

/*  Message  length  verification.  ("CJ"  message  is  a  special  case. 
See  GPSINIT.H.) 

*1 

if  (len  I  -  msg->len  &&  byte2  I  -  X’  &&  bytaS  I  -  T  ) 

{ 

status  -  FALSE; 
return  status; 

}i 


r  Message  check  sum  verification. 

•/ 

if  (  byteCS  I  «  Check  sum(  msg  ) ) 

{ 

status  -  FALSE; 
return  status: 

); 

r  Message  begins  with  @@  prefix. 

•/ 

if  (  byteO  I  «  '(g)'  1 1  bytel  !-'(§)') 
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{ 

status  =  FALSE; 
return  status; 


}; 


I*  Message  ends  with  CR-LF. 

•/ 

if  (  byteCR  I  -  OxOD  1 1  byteLF  I  -  OxOA  ) 

{ 

stitus  ••  FALSE; 
return  status; 

}) 

return  status; 


/ 


ParseJnbuffO 

This  function  accepts  a  pointer  of  type  PANDL  to  the 
contents  extracted  from  the  serial  port  Rx  buffer.  As 
this  function  repeatedly  parses  each  succeeding  massage 
from  this  buffer,  the  buffer's  length  decreases  until 
the  only  remaining  information  is  a  partial  message. 

The  pointer  is  a  static  pointer  and  wili  not  lose  its 
information  upon  compietion  of  this  function  so  it  will 
be  rejoined  with  the  rest  of  itself  when  the  next  grab 
from  the  Rx  buffer  occurs. 

Inputs  are  a  pointer  of  type  PANDL  and  a  pointer  to  a  flag 
indicating  when  the  only  remaining  contents  constitute 

a  partial  message. 

Output  is  a  pointer  to  a  new  buffer  of  type  PANDL  containing 
the  next  complete  single  message  in  the  buffer,  Also, 
the  inbuff  pointer  passed  into  Parsejnbuff  has  been 
altered. 


PANDL  sparse  InbuffI  PANDL  "start  of  next  message,  int  "partial  mag  flag  ) 

{ 

Int  I  ■  0, 1  “  0; 

int  msg_8tartJound_before  -  FALSE; 

PANDL  "parsed  msg; 

ONEBYTE  •ptrjool; 

parssd_msg  ■  malloci  sizeofi  PANDL  ) ); 

ptrjool  -  start_of_next_message->ptr; 

/*  CASE  #1 :  A  complete  message  has  been  found.  The  pointer 
tool  points  to  the  beginning  of  a  msg  (indicated  by  "©(g)'') 
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and  "(s)(§)"  has  been  found  once  before.  Move  the  place  keeper 
pointer  "start_of_next_msa->ptr"  and  adjust  the  length 
"8tart_of_next_m8g->len".  Copy  the  "found"  message  to  a 
new  buffer  and  return  its  location  in  a  pointer  of  type 
PANDL. 

*/ 

for  ( i  -  0;  i  <  start  of  next  message- >len:  !•«-+) 

{ 

If  (  ptr  tooim  -  -  '(S>'  &&  ptr  tool[l  +  11  -  -  I 

{ 

If  (  msg  start  found  before  I 
parsed.mig- >  len  »  I  -  j; 

parsed~msg->ptr  *  calloc(parsed_m8g->len,  sizeoflONEBYTE)); 
memcpy(par8ed_msg->ptr,  ptrjool,  parsed_m8g->len); 

start_of_next_mas8age->ptr  +  ■»  I; 
itart^of^next^mesiage-  >  len  *  -  1; 

*partiaLmsgJlag  ■  FALSE: 
return  parsed  msgi 

}) 

msg  start  found  before  TRUE; 

J  -  ij 

}i 

)} 

r  The  Stan  of  a  second  message  was  not  found.  The  contents 
are  either  exactly  a  complete  message  or  they  constitute  a 
penial  message.  Thus,  the  Parsejnbuff  function  Is 
given  enother  chance  otherwise  the  "panial-msgjiag"  is 
returned  and  the  panlal  message  contents  will  be  prepended 
to  the  next  buffer  brought  Into  Parse  Inbuff. 

•/ 

r  CASE  #2;  There  is  exactly  one  message  in  the  buffer. 

•/ 

If  (  ptrjooKOl  -  -  '®'  &&  \ 
ptr_toolI1 1  ■  -  &&  \ 

ptrjooK  8tart_of_next_mes8age->len  ■  2  1  ■  -  OxOD  &&  \ 
ptr  tool!  start  of  next  message- >  len  •  1  )  -  OxOA  ) 

{  ■ 

parsed_m8g->len  ■  I  •  j; 

par8ed_msg->ptr  «  calloc(par8edjnsg->len,  sizeoflONEBYTE)); 
memcpy(paried_m8g->ptr,  ptrjool,  par8ed_m8g->len); 

8tan_of_next_me88age->ptr  -  NULL; 

8tart_of_next_me8sage->len  -  0; 

"partial  msg  flag  *  FALSE; 

} 
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/*  LA;.E  There  is  only  a  par'iial  rviesragb  left  in  the  buffer. 

*/ 

•Isa 

{ 

^partial  msg  flag  -  TRUE; 

)J 

return  parsed  msg; 

}; 


Maater^gpsO 

For  use  with  a  Motorola  GPS  reoaiver  configured  in  the 
master  station  mode.  It  builds  two  binary  files  on  disk; 

One  records  all  output  messages  and  the  other  records  Just 
position  massages.  The  function  returns  a  pointer  of  type 
PANOL  that  points  to  a  buffer  with  psuedorange  correction 
messages  ready  for  transmission  to  the  slave  receiver. 

Input  Is  nothing. 

Output  is  a  pointer  of  type  PANOL. 


PANDL  *  Master  gps(void) 
{ 


static 

static 

static 


PANDL  mxsbuff; 

ONEBYTE  mxtOOO);  /*  buffer  for  partial  msg  */ 
ONEBVTE  matart  -  TRUE  ; 

ONEBYTE  •tall; 

PANDL  *inbuff,  ‘tmpbuff,  •paraejool,  *000; 

int  *partlal.m8g_flag; 


partial_m8g_flag  -  malloci  8izeof(  ONEBYTE  ) ); 
•partlal_m8g_flag  -  TRUE; 


inbuff  -  malloci  sizeofi  PANDL  )  ); 
tmpbuff  -  malloci  sizeofi  PANDL  )  I; 
parsejool  >  malloci  sizeofi  PANDL  ) ); 
ben  malloci  sizeofi  PANDL  )  ); 


/*  Determine  amount  of  info  in  Rx  buffer  and  import  to  tmpbuff. 

•/ 

tmpbuff->len  ■  iONEBYTE)  siojqueuei  MGPS_PORT  ); 
tmpbuff- >ptr  -  calloci  tmpbuff- >len,  sizeof iONEBYTE)  ); 
8lo_readi  MGPS_PORT,  tmpbuff- >ptr,  tmpbuff- >len  ); 

/*  Initialize  the  excess  buffer  that  is  to  contain  the  partial  msg. 

•/ 

if  i  mstart ) 

{  /*  This  "if"  statement  Is  invoked  */ 
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mxsbuff.len  «  0;  /*  only  the  first  time  the  program  •/ 

mstart  =  FALSE;  /*  is  run.  Note  the  scope  of  mstart.*/ 
mxsbuff.ptr  «  mxs; 

}; 


/*  Append  tmpbuff  messages  to  the  partial  message  that  is  in  mxsbuff. 

*/ 

inbuff- >len  -  tmpbuff- >len  +  mxsbuff.len: 
inbuff- >ptr  >  calloci  inbuff- >len,  sizeoffONEBYTE) ); 
memcpyi  inbuff- >ptr,  mxsbuff. ptr,  mxsbuff.len  ); 
tail  -  Inbuff- > ptr  +  mxsbuff.len; 
memcpyi  tail,  tmpbuff- >  ptr,  tmpbuff- >len  ); 

/*  Free  tmpbuff  memory  and  reset  partial  message  buffer  length  to  zero. 

•/ 

free!  tmpbuff- >  ptr ); 
free!  tmpbuff ); 
mxsbuff.ptr  «  mxs; 
mxsbuff.len  •«  0; 

/*  Initialize  "ben*  for  use  in  collecting  pauedorange  correction  messages. 

*/ 

ben  ox  inbuff; 
ben->len  >■  0; 

while  (  TRUE  ) 

{ 

Get  a  full  message  from  Inbuff. 

*/ 

parsejool  *>  Parse JnbuffI  Inbuff,  partlaLmsgJIag  ): 
if  (  ^partial  msg  flag  ) 

{ 

free!  parsejool- >  ptr ); 
break; 

)1 


/*  Verify  the  received  message  validity  and  indicate  when 
the  message  fails  verification. 

*1 

if  (  Msg  verification!  parse  tool,  RX  )  l>  TRUE  ) 

{ 

printfC Message  from  master  receiver  fails  verificatlonhn”); 

free!  parsejool- > ptr  I; 

continue; 

)i 


/*  Store  only  psuedorange  corrections  in  old  inbuff  buffer  for 
return  to  the  calling  routine. 

•/ 

if  (parsejool- >ptrl21  »  «»  ’C'  &&  parsejool- >ptr( 3 1  -  “  'o'  ) 

memcpy(ben->ptr  -f-  ben->len,  parsejool- > ptr,  parsejool- >len); 
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/*  Write  all  messages  to  this  file. 

*/ 

m_all_file  «  fopeni  "mall.txt",  "ab"  ); 

fwritalparse  jool-  >  ptr,sizeof(ONEBYTE),parse_tool-  >  len,m_all_file); 
fcloael  m_all_file ); 

/*  Write  only  position  messages  to  this  file. 

*/ 

If  (parse  tool->ptr[2]  ->»  'B'  &&  parse  tool->ptr[3|  -  -  'a'  I 

{ 

m_po8_flle  -  fopeni  "mpoa.txt",  "ab"  I; 

fwrite(parse_tool->ptr,ai2eof(ONEBYTE),pars8_tool->len,m^po8_file); 
fcloael  m  pos  file  ); 

h 


free!  parsejool->ptr ); 

);  /*  end  while  */ 

/*  Store  any  partial  message  in  the  excess  buffer, 

•/ 

if  I  Inbuff*  Men  >  0  ) 

{ 

memcpylmxibuff.ptr,  Inbuff- >ptr.  inbuff* >  ten); 
mxsbuff.len  ■  inbuff*  >len; 

)i 


frael  inbuff  )) 
free!  paraa.tool  )i 
free!  partial^magjlag  ); 

return  ben; 

}; 


/ 


Slave_gps() 

For  use  with  a  Motorola  GPS  receiver  in  the  slave 
station  mode.  It  converts  five  data  messages  into 
scaled  and  formatted  data; 

Position/Status/Data  message 
xDOP  Table  Status  message 
Ephemerls  Data  Output  message 
Satellite  Range  Data  Output  message 
Position/Status/Data  Extension  message 
In  addition,  these  and  all  other  output  messages  are 
returned  to  the  calling  routine  for  transmission  to 
the  master  receiver  for  post-processing. 

Input  is  a  pointer  of  type  PANDL  with  any  input  mo.ssages. 
Output  is  a  pointer  of  type  PANDL  with  all  output  messages. 
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/ 


PANDL  'Slave  gps(  PANDL  'input) 

{ 

static  PANDL  sxsbuff; 

static  ONEBYTE  sxsOOOl;  /*  buffer  for  partiai  msg '/ 
static  ONEBYTE  sstart  ">  TRUE  ; 

ONEBYTE  'tail; 

PANDL  'inbuff,  'tmpbuff,  •par8e_tool,  *ben; 
int  *partial_msg_flag,  i; 

/*  Deliver  any  input  messages  to  the  GPS  receiver. 

•/ 

sio_putb(  SQPS_PORT,  input- >ptr,  input- >len  ); 
freei  Input- >ptr”’ ); 
free!  input ); 

/*  Create  flag  pointer  and  working  pointers  of  type  PANDL. 

*/ 

partial_msg_flag  ••  mallocf  sizeof (ON EBYTE)  ); 

•partial^msgjlag  -  TRUE;  /♦  Used  in  Parsejnbuff  */ 

inbuff  ••  malloci  sizeofi  PANDL  )  ); 
tmpbuff  -  malloci  sizeofi  PANDL  I  ); 
parsejool  malloci  sizeofi  PANDL  ) ); 
ben  «  malloci  sizeofi  PANDL  )  ); 

/•  Determine  amount  of  Info  in  Rx  buffer  and  Import  to  tmpbuff. 

•/ 

tmpbuff->len  -  (ONEBYTE)  slojqueuei  SGPS.PORT  ); 
tmpbuff- >ptr  -  calloci  tmpbuff- >lon,  sizeof (CJNEBYTE)  ); 
sio_read(  SGPS_PORT,  tmpbuff- >ptr,  tmpbuff- >len  ); 

r  Initialize  the  excess  buffer  that  ia  to  contain  the  partial  msg. 

•/ 

if  (  sstart ) 

{  /•  This  "if"  statement  Is  invoked  '/ 

sxsbuff. len  *  0;  /'  only  the  first  time  the  program  */ 

sstart  »  FALSE;  /'  is  run.  Note  the  scope  of  sstart.*/ 
sxsbuff. ptr  ■■  sxs; 

}: 


/*  Append  this  info  in  tmpbuff  to  excess  msg  from  last  read. 

•/ 

inbuff- >  len  ■  tmpbuff- >  len  +  sxsbuff. len; 
inbuff- > ptr  -  calloci  Inbuff- > len,  sizeofiONEBYTE) ); 
memcpyi  Inbuff- >ptr,  sxsbuff. ptr,  sxsbuff. len  ); 
tail  -  inbuff- > ptr  +  sx.sbuff.len; 
memcpyi  tail,  tmpbuff- > ptr,  tmpbuff- >len  ); 

/*  Free  tmpbuff  memory  and  reset  partial  message  buffer  length  to  zero. 
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*/ 

free!  tmpbuff->ptr  ); 
free(  tmpbuff  ); 
sxsbuff.ptr  -  $xs; 
sxsbuff.len  >  0; 

/*  Initialize  "ben”  for  use  in  returning  all  messages  to  the  calling 
environment. 

•/ 

ben  -  inbuff; 

while  ( TRUE  ) 

{ 

I*  Get  a  full  message  from  inbuff. 

*/ 

parsejool  «  Parsejnbuffi  inbuff,  partial_msg_flag  ); 
if  I  *purtial  msg  flag  ) 

{ 

free!  parsejooi->ptr ); 
break; 


/*  Verify  the  received  message  validity. 

•/ 

if  <  Msg  verification!  parse  tooi,  RX  )  I  -  TRUE  ) 

{ 

free!  par$e.tool->ptr ); 
continue; 

}J 


/*  Convert  the  binary  rnsg  to  GPS  data  structure  for  use  by  cpu. 

*/ 

switch  (parse  tool->ptr(3]) 

{ 

case  'a':  /*  ®@Ba  Position/Status/Data  */ 

/*  insert  conversion  here  */ 
break; 

case  'c';  /*  ©^Bc  xDOP  Table  Status  */ 

/•  Insert  conversion  here  */ 
break; 

case  'f:  /*  (SX^Bf  Ephemeris  Data  */ 

for  ( i  -  5  ;  I  <  EPH  NUM  WORDS  +  5  ;  i+  -3  I 

{ 

hox_ephetiieri8(ll(01  -  parsejool- >ptr(i  +  01; 
hex.ephemerislild )  parsejool- >ptr(i  -t-  1]; 
hex  ephemeris[ll[2]  «  parse  tool->ptr(i  +  2] ; 

) 

/*  Insert  conversion  here  •/ 
break; 

case  'g'l  /*  @@Bg  Satellite  Range  */ 

/•  insert  conversion  here  */ 
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break; 

case  'k':  /*  @@Bk  Position/Status/Data  Cxt 

/*  Insert  conversion  here  */ 
break; 
default: 
break; 


free!  parsejool- >  ptr ); 

};  /•  end  while  */ 

r  Store  any  partial  message  In  the  excess  buffer. 

*/ 

if  ( inbuff"  >len  >  0  1 

{ 

memcpvlsxsbuff.ptr,  lnbuff->ptr,  lnbuff->lan); 
sxsbuff.len  -  inbuff'>len: 

}i 

free!  inbutf  ); 
free!  parse.tool ); 
free!  partiaLmsgJIag  ); 

return  ben; 

}; 


#ondlf  /•  QPSFUN  H  */ 


6 


GPSINIT.H 


/♦  GPEiNIT.H 

This  headet  file  contains  the  Motorola  Proprietary 
Binary  Message  Format  characteristcs  for  GPS 
receivers.  The  first  column  contains  the  second 
letter  of  the  message  ID  of  a  message  prepared  for 
input  to  the  GPS  receiver.  The  fourth  column  contains 
the  second  letter  of  the  message  ID  of  a  message 
received  from  the  GPS  receiver.  For  all  messages  sent 
to  the  GPS  receiver  except  one,  there  is  a  message  reply. 

The  second  column  contains  the  final  length  of  a 
message  prepared  to  be  sent  to  the  GPS  receiver.  The 
third  column  is  the  length  of  a  message  received  from 
the  GPS  receiver  with  that  particular  message  ID. 

Note:  All  message  lengths  can  be  expressed  with 
one  byte  except  one:  The  response  message  Cj  is  294  bytes. 
A  unique  dummy  size  equal  to  1  byte  has  been  used  to 
indicate  the  length. 

*/ 

#ifndef  GPSINIT  h 
#define  ^GPSINIT_H 

#ifndef  _GPSDEFIN_H 

finclude  * c !\borlandc\twite\gpsdef in . h” 

#endif 

0NE3YTE  A_array(301(4l-{{  'a'.  10,  10,  'a'  }, 


{ 

'b' 

,  10, 

10, 

'b'  }, 

{ 

'c', 

,  11. 

11. 

'c'  }. 

{ 

*d' 

,  11, 

11. 

'd'  }, 

{ 

'e', 

,  11, 

11, 

'e'  }. 

{ 

'f', 

12, 

15, 

•r  }. 

{ 

'O' 

,  0, 

8, 

'g'  }. 

{ 

'h' 

,  8, 

8, 

’b'  }, 

{ 

'i', 

9, 

9,  ' 

i'  ), 

{ 

•y. 

8, 

8,  ' 

r  ), 

{ 

'k'. 

,  9, 

9, 

k'  ), 

{ 

'1*. 

9, 

9,  ' 

1'  )- 

{ 

'm 

',  12 

,  12, 

'm'  ) 

{ 

'n' 

,  8, 

8, 

'n'  ), 

{ 

'o' 

,  8, 

25, 

'P’  )- 

{ 

'P' 

,  25, 

25, 

'P'  }. 

{ 

'q' 

,  8, 

8, 

'q'  ). 

{ 

'r', 

8, 

8,  ' 

r'  }, 

{ 

■s' 

,  20, 

20, 

's'  }, 

{ 

't', 

8, 

8,  ' 

t'  }, 

12y 


ONEBYTE 


ONEBYTE 


#endif  /* 


{ 

'u', 

12, 

12.  'u'  }, 

{ 

'v'. 

8, 

8,  'v'  }, 

( 

'w' 

.  B, 

8,  'w'  }, 

{ 

'x', 

9. 

9,  'X'  }, 

{ 

'V'. 

11. 

11,  V  }, 

{ 

'2', 

11, 

11,  '2'  }, 

{ 

'A' 

.  8, 

8,  'A'  }, 

{ 

'B', 

8, 

8,  '8'  }, 

{ 

'C, 

,  9, 

9.  'C'  }, 

{ 

'D' 

,  9. 

9.  'D'  )); 

I  arrav[11l[4 

l-({  'a' 

00 

CD 

(6 

{ 

'b', 

8, 

92,  'b'  }, 

{ 

'c'. 

8, 

82,  'c'  }, 

{ 

'd'. 

8, 

23,  'd'  }, 

{ 

'a'. 

8. 

0,  },/* 

Be*>Cb  •/ 

{ 

'i', 

8. 

80,  r  }, 

{ 

'O'. 

8, 

122,‘g' 

{ 

'h'. 

8. 

0,  },/• 

Bh->Ca  */ 

{ 

T. 

7. 

8,  'j'  }, 

{ 

r. 

80, 

0.  },/• 

Bf->Cc  •/ 

{ 

•k\ 

8. 

CD 

:  arrayll  11141-  {{  'a',  7, 

9,  'a'  }, 

{ 

'b\ 

33, 

9,  'h'  }, 

{ 

'd'. 

27, 

171,'d'  }, 

{ 

62, 

7,  'k'  }, 

{ 

T, 

7, 

7,  T  }, 

{ 

'o'- 

8. 

8,  '0'  }, 

{ 

'i'. 

a. 

0,  M'  ), 

{ 

'j'. 

7, 

1,  T  )./• 

CJ  rasp  msg  294  byte 

{ 

0, 

33,  'b*  },  / 

•  BB'>Cb  */ 

{ 

0. 

62,  V  },  / 

•  Bh->Ce  */ 

{ 

0, 

80.  'c'  )};/ 

*  Bf->Cc  •/ 

GPSINIT_H  */ 
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H.  GPSPORTS.H 


/•  GPSPORTS.H 

lnitialize_master_gpsO 

lnitialize_slave_0P$0 

Open_ports(l 

Cloae.portsO 

Both  initialization  functions  need  to  be  greatly  improved 
to  provide  an  easy  to  understand  and  use  interface  with 
the  user. 


#ifndef  QPSPORTS  H 
#define  .GPSPORTS.H 

#include  <stdio.h> 

#inciude  ''c;\borlandc\twite\head-c.h*' 


Initialize.master.gpsO 

This  function  presents  the  current  configuration 
settings  for  the  master  GPS  receiver.  The  user  Is 
offered  the  opportunity  to  change  any  of  the  parameters. 
In  addition,  the  user  is  requested  to  indicate  which  if 
any  output  files  he  will  want.  For  those  files,  the  default 
file  name  is  presented  and  made  available  for  modification 
or  change. 

Input  is  the  serial  port  number  through  which  the  master 
GPS  communicates. 

Output  is  nothing. 


void  initialize  master  gpsiint  mgps) 

{ 

char  a1(20l  -{'@'/@'/B','aM,?.2,13,10}; 
char  s2(201-{'(g)','@','B','hM,43,13,10); 
8io.flush(mgps,2); 
sio.putb(mgps,s1 ,8); 
sio_putb(mgps,s2,8); 

printfC Master  GPS  receiver  initialized. \n''); 
return; 

); 


/ 
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lnitialize_$lave_gpsO 

This  function  presents  the  current  configuration 
settings  for  the  slave  GPS  receiver.  The  user  is 
offered  the  opportunity  to  change  any  of  the  parameters. 
In  addition,  the  user  is  requested  to  indicate  which  if 
any  output  files  he  will  want.  For  those  files,  the  default 
file  name  is  presented  and  made  available  for  modification 
or  change. 

Input  Is  the  serial  port  number  through  which  the  slave 
GPS  communicates. 

Output  is  nothing. 


V 


void  Initialize  slave  gpslint  sgps) 

{ 

char  8lt201«{‘®','@',’B',‘aM, 22,13, 10}; 

slo_flush(sgp8,2): 

sio_putb(Bgp8,s1 ,8); 

printf(*Slave  GPS  receiver  initlalized.\n”); 
return; 

): 


/ . * . * . . . . 

Open^DortaO 

This  function  initializes  the  two  serial  ports 
on  the  PCD-744  to  be  used  by  the  master  and  slave  GPS 
receivers.  Before  initialization,  it  either  loads  the 
drivers  for  the  standard  and  PCD-744  serial  ports  or 
verifies  that  they  are  loaded.  It  then  configures  the 
Baud  rating,  word  length,  parity,  and  stop  bits  as  well 
as  the  modem  and  flow  control.  All  commands  that  begin 
with  "sio^..."  are  PCD-744  commands  and  may  be  referenced 
in  the  PCLS-B02  PC-ComLIB  Serial  Communication  Programming 
Library. 

The  PCD-744  RS-232  8  port  intelligent  10  card  is 
made  by  ADVANTFCH  in  Sunnyvale,  Ca.  14081  24&'6670. 

Input  is  the  slave  GPS  port  number. 

Output  Is  nothing. 

V 

void  Open  ports!  int  mgps,  int  sgps  ) 

int  BAUD  -12;  /•  9600  (p.  12  PC-ComLIB)  */ 

■It  lOMODE  -  (0x0310x0010x001;  /*  8-N-1  (p.  12)  •/ 

nt  MODMODE  -0x00;  /*  DTR  and  RTS  off  (p.26)*/ 
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int  HWMODE 
int  n; 


=  0x00; 

(p.33r/ 


/*  HW  and  SW  flow  Ctrl  off 


n  -  sio  ioctK  mgps,  BAUD,  lOMODE  I; 
if  (n1  -  0) 

{ 

printfCOpening  comport  %d  iocti  error.\n'',mgps); 
BbortO; 

); 

n  >  sio  Ictril  mgps,  MODMOOE  ); 
if  (n1  -  0) 

{ 

printf(*Opening  comport  %d  ictrl  error.\n'',mgps); 
abort!): 

}; 

n  »  sio  flowctri!  mgps,  HWMODE  ); 
if  (n~!  -  0) 

{ 

printfCOpening  comport  %d  flowctri  error.\n”,mgps); 
abort!); 

}j 

n  >>  sio  open!  mgps  ): 
if  !n1  -  0) 

{ 

printf!"Openlng  comport  %d  error.\n*,mgps); 
abort!); 

}! 


n  >  sio  iocti!  sgps,  BAUD,  lOMODE  ); 
if  !n1  -  0) 

{ 

printf!''Opening  comport  %d  iocti  error.\n*,sgps); 
abort!): 

}; 

n  -  sio  Ictri!  sgps,  MOOMODE  ); 
if  !n1  -  0) 

{ 

printf!"Open'ng  comport  %d  Ictrl  error.\n",8gp5); 
abort!); 

}: 

n  ••  sio  flowctri!  sgps,  HWMODE  ): 
if  !n‘! «  0) 

{ 

printf!''Opening  comport  %d  fiowctrl  error. Nn", sgps); 
abort!); 

}; 

n  -  sio  open!  sgps  ): 
if  !n"l  -  0) 

{ 
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printfrOpening  comport  %d  error. \n",sgps); 
abort!); 

}; 

printf(”Opened  ports  %d  (Master  GPS)  and  %d  (Slave  GPS).\n",  \ 

mgps,  sgps): 


return; 

}; 


/ 


Close_ports() 

This  function  closes  and  resets  all  serial  ports  for 
the  PCD*744  Serial  10  card.  See  the  Open_ports  function 
for  more  information  on  the  commands  used  herein. 

Inputs  are  the  master  and  slave  GPS  port  numbers. 
Output  is  nothing. 


void  Close  ports!  int  mgps,  int  sgps  ) 

{ 

printfl" Closing  all  ports\n''); 
aio_Glose(  sgps ); 
aio^elosel  mgps  ); 
aiojesetO; 
return; 


}; 


#endlf  /•  ,GPSPORTS_H  */ 
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I 


GPSTRUCT.H 


r  GPSTRUCT.H 

•/ 

#lfnddf  GPSTRUCT  H 
fdefine  ’GPSTRUCT.H 

#ifnd*f  .GPSINIT.H 

#include  ''c!\borlandc\twite\gp8lnlt.h'' 

#andif 

typadaf  atruct{  I*  Pointer  AND  Len  */ 

ONEBYTE  *ptrj 
int  len) 

}PANDL; 

t*  Struct  and  global  data  dafinitione 
V 

struct  T  FLOAT  EPHEMERIS  { 

UNSIGNED.ONEBYTE  lode; 

/*  Field  name  abbreviations  were  extracted  from  ICD‘GPS‘200. 

*1 

double  toe)  /*  Time  of  ephemeris.  */ 

double  mO)  /*  Mean  anomaly  at  refernce  time.  *1 

double  delta.n;  /*  Mean  motion  difference  from  computed  value.  *1 

double  e;  ’  I*  Eccentricity.  */ 

double  aqrt.a;  /*  Square  root  of  aemLmajor  axis.  */ 

double  omsga.O;  I*  Longitude  of  ascending  node  of  orbit  plane 
at  weeKly  epoch.  *1 

doublo  iO)  /*  Inclination  angle  at  reference  time.  *1 

double  w;  /*  Argument  of  perigee.  */ 

double  omega.dot;  /*  Rate  of  right  ascension.  */ 

double  l.dot;  /*  Rate  of  inclination  angle.  *1 

double  cue;  /*  Amplitude  of  the  cosine  harmonic  correction 

term  to  the  argument  of  latitude.  */ 

double  cus;  /*  Amplitude  of  the  sine  harmonic  correction 

term  to  the  argument  of  latitude.  */ 

double  crc)  /*  Amplitude  of  the  cosine  harmonic  correction 

term  to  the  orbit  radius.  */ 

double  era;  /*  Amplitude  of  the  sine  harmonic  correction 

term  to  the  orbit  radius.  */ 

double  cic;  /*  Amplitude  of  the  cosine  harmonic  correction 

term  to  the  angle  of  Inclination.  */ 

double  cis;  /*  Amplitude  of  the  sine  harmonic  correction 

term  to  the  angle  of  inclination.  */ 

double  tgd;  /*  Estimated  group  delay  differential.  */ 
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double  toe; 
double  af2; 
double  af1; 
double  afO; 
double  n; 
double  e.eqrt; 

double  wk1 ; 
double  wkO; 
double  aln.w) 
double  ooe.w; 
double  a; 
double  sqrta.ej 

}  float.ephemeria 


/*  Clock  data  reference  time.  */ 

/*  Polvnomial  coefficient  (SV  clock  correction)  */ 
/•  Polynomial  coefficient  (SV  clock  correction)  */ 
/*  Polynomial  coefficient  ISV  clock  correction)  */ 
/*  Corrected  mean  motion.  •/ 

r  Square  root  of  eccentricity‘'2  aubtracted 
from  1.  •/ 

/*  Omoga_dot  •  WE.  */ 

/*  Omega  0  ■  WE  •  toe.  */ 

/•  Sine  wr  */ 

/*  Coiine  w,  •/ 

/*  Semi-maicf  axis.  •/ 

jeonst;  /*  Squarb  root  of  semi*major  axis  * 
eccentricity  *  PCONST.  *! 


struct  T  RECEIVER  CHANNELS  { 
UNSIGNED  ONEBYTE  avid: 
UNSIGNED~ONEBYTE  mode; 
UNSIGNED  ONEBYTE  atrength; 
UNSIONED~ONEBYTE  flags; 


struct  T  GEODETIC  { 

TWOBYTE  degrees ; 
UNSIGNEDJWOBYTE  rninutea  ; 
double  ~  seconds ; 


struct  T.,POS.CHAN_STATUS  { 


UNSIGNED  ONEBYTE 
UNSIGNED  ONEBYTE 
UNSIGNED  TWOBYTE 
UNSIGNED  ONEBYTE 
UNSIGNED_ONEBYTE 
double 

struct  T  GEODETIC 

struct  T.GEODETIC 

double 

double 

double 

double 

double 

UNSIGNED  ONEBYTE 
UNSIGNED  ONEBYTE 
UNSIGNED  ONEBYTE 


month) 

day; 

year; 

hours; 

minutea; 

seconds: 

latitude; 

longitude; 

datum_helght:  /*  meters  */ 
m8l_height;  /*  meters  */ 
velocity:  /*  m/sec  */ 

heading;  /*  degrees  */ 
current_dop: 
dopjype; 
visible  sets; 


UNSIGNED  ONEBYTE  sate  tracked; 

struct  T  RECEIVER  CHANNELS  channeKNUM  CHANNELS]  ; 

UNSIGNED  ONEBYTE  tcvr  status; 


#endlf  /*  GPSTRUCT  H  */ 
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HEAD-C.H 


/*  moxa-c.h  verl.OO  10/11/1991 
*  Definitions  for  MOXA  Seriai  I/O  Controller  ioctrl 

•/ 

/*  BAUD  rate  setting  */ 


#define  B50 

0x00 

#define  B75 

0x01 

#define  B110 

0x02 

#deflne  B134 

0x03 

#define  B150 

0x04 

#define  B300 

0x05 

Idefine  B600 

0x06 

#define  B1200 

0x07 

^define  B1800 

0x08 

Idefine  B2400 

0x09 

^define  B4800 

OxOA 

^define  B7200 

0x08 

#deflnt  B9600 

OxOC 

#define  819200 

OxOD 

#dafine  838400 

OxOE 

#define  857600 

OxOF 

1*  MODE  setting 

*/ 

^define  BIT  5 

0x00 

/*  Word  length  define 

#d«fine  BIT'6 

0x01 

#define  BIT  7 

0x02 

#def!ne  BITJ 

0x03 

#define  STOP  1 

0x00 

/*  Stop  bits  define  *1 

#difit»««  STOPl2 

0x04 

fdefine  P  EVEN 

0x18 

r  Psrlty  define  •/ 

#define  P^ODD 

0x08 

#dflfine  P  SPC 

0x38 

#define  p’mrK 

0x28 

#define  P~NONE 

0x00 

/•  MODEM  CONTROL  setting  */ 
#deflne  C  DTP  0x01 

fdefine  C_RTS  0x02 

/•  MODEM  LINE  STATUS  */ 
#define  S  CTS  0x01 

^define  S  DSR  0x02 

^define  S  Rl  0x04 
fdefine  S_CD  0x08 

/*  function  declaration  •/ 
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long  siojqueuelint); 
long  siojfree(int); 
long  sio.oquauednt); 
long  sio.ofreelint); 


APPENDIX  F:  NP8  QEODETIC  MARKS 


SPANAGEL  GEODETIC  MARKS 

WGS-84 


Marks 


8PANAOXL  HAAL 


2  3  1 

Laval  B 


mmm^^  80*15  DS^ 


spanagal  Hall  runs  SK  to  HE.  Tha  narka  ara  on  tha  8  th 
laval  naar  tha  North-Eaat  oornot  of  tha  building. 

Marks  1  and  2  ara  standard  aurvay  disks «  markad  SPA-l  and 
SPA~2.  SPA-3  is  a  aquara  of  thin  aluminum.  All  marka 
ara  attaohad  to  tha  roof  undar  a  rail. 

Polaa  ara  mountad  on  tha  rkil  abova  tha  marks.  Tha 
haight  to  tha  marks  ara  11. B3  m  (  WGS-84  alipioid  ). 


Mark  Pole 


SPA-1 

36  35.70375 

121  52.47777 

11.83 

14.08 

SPA-2 

36  35.70180 

121  52.48074 

11.83 

14.10 

SPA-3 

36  35.70277 

121  52.47926 

11.83 

14.09 
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Tabla  IX  WGS  84  spanagal  hall  eighth  floor  geodetic  marks. 


Loca- 

tion 

Ztaa 

1 

dagraaa 

dagraaa 

ainutaa 

dagraaa 

nlnutas 

aaoonds 

matars 

(abova 

gaoid) 

Spa-**.! 

;'Lat' 

i"!  -• 

ill.. 

;36.59506253 

36“  35.70375’ 

36“  35’ 
42.2251” 

2321 

Lon 

121.87462956 

“W 

121“ 

52.47777’ 

121“  52’ 
28.6664” 

Ait 

11.83 

Spa’>2 

Lat 

36.59503006 

“N 

36“  35.70180’ 

36“  35’  , 
42.1082” 

2322 

Lon 

121.87467903 

*W 

121“ 

52.48074’ 

121“  52’ 
28.8445” 

Alt 

11.83 

Spa-3 

Lat 

36.59504622 

‘’N 

36“  36.70277’ 

36“  35’ 
42.1664” 

2323 

Lon 

121.87465439 

ow 

121“ 

52.47926’ 

121“  52’ 
28.7558” 

Alt 

11.83 
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